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PREFACE 


This book treats, in a self-contained manner, the problem of the 
deformation of an clastic solid without making the assumption of the 
classical infinitesimal, or linear, theory of elasticity that the deforma- 
tion is so small that, the squares of the strain components are negligible. 
Matrices are used methodically, and an account of those properties 
of matrices which are used in the book is given in the first chapter. 
The importance of non-isotropic media is emphasized, and the fact 
that applied stress can make non-isotropic a medium which is isotropic 
when free from stress is pointed out. The relations between the 
third-order, as well as the second-order, elastic constants for various 
typos of crystalline media arc derived. The agreement between the 
predictions of theory and Bridgman’s experiments on the compressi- 
bility of media under hydrostatic pressures up to 10 5 atmospheres is 
shown. The final chapter applies the theory to the problem of the 
deformation of spherical shells and circular tubes under extreme 
external and internal pressures. 

Wo have, so far, very little experimental knowledge of the third- 
order elastic constants or of the effect of applied stress upon the second- 
order elastic constants. If the mathematical treatment given here 
serves to stimulate the procurement of experimental knowledge of 
those phenomena we shall have attained our aim. 

Part of the material of the book was given in lectures at the Dublin 
Institute for Advanced Studies (summer, 194S), the Carnegie Institute 
of Technology (fall, 1918), and the Eseola Tecnica do Exercito, Rio de 
•lanciro (1919). Wc have thought it well to let the book maintain 
the elementary form of treatment adopted in these lectures and have 
not attempted to provide references to other treatments of the topics 
discussed in thc'boqh.- Owing, doubtless, to technological demands, 
these topics arc currently much discussed, and the ‘interested reader 
may find many references in the latest issues of such abstract journals 
ns Mathematical R crimes and Science Abstracts, 

The inscription that closed the preface to my previous book in the 
Wiley Applied Mathematics Tnfrbdvction .to .Applied Mathe- 
matics, may well ho repented here: • , 

To the Glory of God. Honor of Ireland, 
and 

Solidarity of the Americas. 


I'te it Jarci'-.i 

Art..:' 


F. D. Mt-RNAGIIAX 
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VECTORS AND MATRICES 


The theory of finite deformations of an elastic solid is most easily 
presented and understood by the use of matrices. Very little beyond 
the elements of the theory of matrices is needed, but since many 
engineers who may be interested in the theory of finite deformations 
may not be familiar with matrices this introductory chapter gives a 
self-contained and elementary account of that part of the theory of 
matrices which is necessary for our purpose. 

1. The iuntri\ concept 

The concept of a matrix arises naturally and inevitably when we con- 
sider several different iablc functions of several variables. In ele- 
mentary calculus we learn that, if u = «(x) is a differentiable function 
of a single independent variable x, then the differential du of u is 
obtained by multiplying the differential dx of x by the derivative tt x of 
» with respect to x: 

du — iij dx. 

Thus the derivative it- of it with respect to x may be regarded as a 
magnification factor flint converts, by multiplication, the differential 
dx of X into the differential du of u. When we pass to the considera- 
tion of ttro differentiable functions (u, v) of two independent variables 
(x, it) we have two differentials (du, dv) that arc connected with the two 
differentials (dx, dy) of (lie independent variables by the formulas 

du = u z dr -r u y dy. 

dr ~ r- dx + r v dy. 

Ibuv a, and u 7 are the partial derivatives of u = u(r, y) with respect 
to x and ?/. respectively, and r. and t\. are (he partial derivatives of 
i i' r\s, */) with respect to x and y, respectively. The collection of 
four partial derivatives (»-. u s , r r , r ff ) now plays the role previousb' 
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plajed bj the one derivative and this collection of four denvatives 
ma> be considered as a magnification /actor that converts, by multipli- 
cation, the pair of differentials (dx, dy) into the pair of differentials 
(rfu, di) It e regard the pair of differentials {dx, dy) as a single entity, 
the differential lector of the independent \ ariabtes, and, similarly, we 
regard the pair of differentials {du, dv) as a single entity, the differen- 
tial vector of the dependent variables We denote these differential 
vectors b> hold-face type, thus dx is the differential vector of the 
independent variables and {dx, dy) are the coordinates of dx. In 
s j mbols 

dx =* i(dx, dy) 


(the translation of this stenographic notation being dx is the vector 
whose coordinates are dx and dy) Similarly , 

du = t {du, do) 


(translate this stenographic symbol) It is convenient to write the 
coordinates of dx in a column, i e one abov e the other, thus 


(X) 


We refer to this sjmbol as a 2 X 1 matrix , i e , a matrix of two rows 
and one column The first coordinate dx of the vector dx is the 


element in the first row of the matm 


©■ 


and the second coor- 


dinate dy of the vector dx is the element m the second row of tho 
Similarly, wc write the coordinates of du as the elc- 


(S) 


mentfl of a 2 X 1 matm 


e> 


and then the collection of four 


drnvatives (u„ u # , t*, i v ) maj be regarded ns constituting a mngnifica- 


fdx\ 

c W 


into the 2X1 matm I 


write the four derivatives as the 


turn factor that converts, bj multiplication the 2 X i matrix ( 

[ (t) Wc " n 

elements of a 2 X 2 matrix (ic , a matm of two rows and two col- 
umns), and wc denote this matm bv the svmbol 

(*, V) 

(«, t>) 


(tf, v) m / V , u A 
(X, y) ~ v* *v) 
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Note that the upper letters of the symbol 

the columns, of the matrix 


tell the rows, 


whereas the lower letters tell the columns, of the matrix. Thus 
the elements of the first row are the derivatives of u with respect to 
x and y, and the elements of the second row are the derivatives of v 
with respect to x and y. On the other hand, the elements of the first 
column are the derivatives of u and v with respect to x, and the ele- 
ments of the second column are the derivatives of v and v with respect 

y. This 2X2 matrix is known as the Jacobian matrix 


to 


(*» ?/) 


(after K, G. J. Jacobi [1801-1851], a German mathematician) of the 
pair of functions («, i>) with respect to the pair of independent variables 
(x, y); we denote it by the symbol ./. We now introduce a shorthand 
notation in which we denote any one-column matrix by the symbol 
for the element in its first row; thus we denote the 2 X 1 matrix 

.‘'imply by <!x and the 2X1 matrix simply by du. In 
this notation the pair of equations 

du = Uj dx -1- u u dy, 
dv — !> x dx + v v dy 
appears in the abbreviated form 

du — ./ dx. 


The 2X2 matrix ./ has two 2 X 1 matrices, or column rectors as we 

shall term them, namely, (^‘^j and and two 1 X 2 matrices, 

or row vector s ns we shall call them, namely, (u x , u v ) and (r x , v u ). We 
term the expression u x dx + t/ v dy (which furnishes the first element of 
the 2 X 1 matrix du) the product of the 2 X 1 matrix dx by the first 
row vector of ./; similarly, the expression v x dx + v v dy (which fur- 
nisbet the second element of the 2 X 1 matrix du) is the product of 
the 2 X 1 matrix dx by the second row vector of In this ter- 
minology. then, the product ./ dx of the 2 X 1 matrix dx by the 2X2 
Jacobian matrix J is a 2 X 1 matrix du whose elements are obtained by 
methodically multiplying the 2 X 1 matrix dx by the row vectors of 
J, one after the other. 

Example. The formulas connecting plane polar coordinates (p, <£) 
amt rectangular Cartesian coordinates (x, y) arc x = p cos 4>,y = p sin 6. 

The Jacobian matrix of (x. y) with respect to ( P , t) is 
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cos <t> — p sin 
^sin 4> p cos 


:) 


and the product of the 2X1 matrix dp = hy ia 


dx 


( cos <frdp — p sm ^ d<f> 
sin <f> dp + P cos <f> d<f> t 


:) 


The first element of this matrix is the differential of x and the second 
clement is the differential of y 

Remark "You should note the analogy between the relation du — 
J dx and the basic relation dy = y x dx of elementary differential 
calculus Realize clearlj that it is the 2X2 Jacobian matrix J that 
plajs for two differentiable functions of two independent variables 
the role played by the derivative of one differentiable function of one 
independent \anable The four elements of J have mdnidually no 
fundamental significance It is the 2X2 Jacobian matrix of which 
they are the elements that is the important thing 
There is no difficulty m extending the theorj of the preceding para 
graphs to three differentiable functions of three independent variables 
or generallj to n differentiable functions of n independent variables 
b== 3 4 5 Thus for three independent variables we denote 

/dx\ /du\ 

bj dx the 3 X 1 matrix I dy J and b> du the 3 X 1 matrix I di J 
\dz/ \dw ) 

T1 c Jacobian matrix J = ~ - is the matnx 
(x y z) 

/ «* u v u A 

•' = ['. ".I 

\w z tr„ w,/ 

of which for example the 1X3 matnx ( u, u„ u,) is the first row 
vector J maj be regarded as a magnification factor that converts 
bj multiplication the 3 X 1 matrix dx into the 3 X 1 matrix du 

du <=> J dx 


In this kind of multiplication of dx by J each element of the 3X1 
matnx du is the product of the 3X1 matnx dx bj the corresponding 
row \ector of J for example the element dv in the second row of the 
3X1 matnx du is gn en bj the formula 


dv = t, dx 4- v y dy + i, dz 
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Example. The formulas connecting space polar coordinates 
(r, D, <>) mid rectangular Cartesian coordinates (x, y, z ) are 

x *= r sin 0 cos <t>t 

y — r sin 0 -sin <f>, 

z ~ r cos 0. 


The Jacobian matrix ./ 


------ of (x, ii, z) with respect to (r, 0, £) is 

(r, 0, <f>) 


J - 


fsin 0 cos <t> 
sin 0 sin d 
cos 0 


r cos 0 cos p 
r cos 0 sin 4> 
— r sin 0 


■ r sin 0 sin 
r sin 0 cos 4> 
0 


T),(. elements of t lie 3 X 1 matrix <lx are furnished by the formula 
dx J dr; thus 


(lx ~ sin 0 cos $ dr + r cos 0 cos 4> do — r sin 0 sin <p d(j>, 


dy — sin 0 sin </> dr T r cos 0 sin </> d0 -f- r sin 0 cos <£ d$, 


dz = cos 0 dr — r sin 0 d0. 


It should be clear from the preceding discussion what is meant by 
an m X n matrix A where m and n arc any two positive integers. This 
is a collection of inn numbers arranged in m rows and n columns. We 
denote the tlnnrnt in the rth row and sth column by the symbol 
a,*. Note carefully that the upper label tells the row and the lower 
label tells the column. In the applications of matrix theory which we 
shall have to make neither m nor n will surpass 3, and we can use the 
readily visualized concepts of one-dimensional vectors (i.c., vectors on 
ft line) or of two-dimensional vectors (i.e.. plane vectors) or of three- 
dimensional vectors (i.e., space vectors). For example, if m = 2, 
n **• 3, so that A has two rows and three columns, .4 possesses three 


plane column vectors of whirh the second, for instance, is rO), and 

\ 0 *V 

it abo posse— e*. t wo space row vectors of which the first , for instance, is 
Wi’i os 1 , oj 1 ’- If ft is an « X p matrix (so that B possesses as many 
nnw sis .1 possesses columns) it is possible to multiply any column 
vector of ft by any row vector of .4 since the column vectors of B and 


the row vectors of .1 are all of the same dimension, natnclv, n. We 
ariH' the e various* products of column vectors of 13 by row vectors 
of A n.s tm »?! X p matrix tlm element in the jth row and Art h column 
of C iK’iac the pr.vluet of the /rth column vector of B by the jth row 
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vector of A, 3 «= 1 , m, k *= 1, , p, and we term C the 

product of the n X p matrix B bj the m X n matrix A We indicate 
this product operation as follows 

C - AB 


EXERCISES 


1. Show that if A is polite m multiplication, i e , if AB is the same aa B tor every 
B, A must be square (i e , hav e the same number of columns as rows) with all it9 
diagonal elements 1 and all its other elements zero 
A ole The diagonal elements of a square matrix A arc the elements ai 1 , fli*, 
whose row and column numbers are equal 

Hint If A vs an m X n matrix and BwinXp matrix, AB is an m X p matrix 
For A B to be the same as B w e must have in = n so that A is square 

Remark The matrices that are polite in multiplication arc known aa unit 
matrices Terming the common number of rows and columns of a square matrix 
its dimension, we denote the unit matrix of dimension n by the symbol E R Thus 



2 Form the follow ing products 

<«>G -?)(n> -D(i _?)• 


to 




Answer <«) ( ”). (6) (c) 


A ole The difference between the answers to (a) and (6) shows that when both 
products AB and BA can be formed (show that each product is a square matrix 
when this happens) and when they have the same dimension (show that A and B 
are each square when this happens) the products arc not necessarily the same 
Multiplication of square matrices of the same dimension is not, in general, commutalue 
The results of (e) and (d) show that if we denote by cE n the n X n matrix whoso 
diagonal elements are all — e and whose non-diagonal elements are zero then if A 
is any p X q matrix each of the products ( cE p )A and A ( cE t ) is obtained by multi- 
ply ing each and every element of A by e Wo denote each of these products 
simply by eA 
3 Form the two products 


Ansuer 


(a) (-4, 1,2) 


(o) (-2',), 




(-4,1,2) 
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1. Show tlmt, if D is a diagonal matrix, DA is obtained by multiplying each row- 
vector of A by the corresponding diagonal element of D whereas AD is obtained by 
multiplying each column vector of A by the corresponding diagonal element of D. 
S' oir, A diagonal matrix D is a square matrix of which every non-diagonal element 
i« zero; thus df 1 « 0 if f: j. 

5, Show that if D\ and />; are any two diagonal matrices of the same dimension 
then D:D\ f ar DiDi; in words, multiplication of diagonal matrices is commutative. 

(n Show that if A is any m X n matrix and It is any n X m matrix, so that both 
AH and HA arc square matrices, then the sum of the diagonal elements of AD is 
the same as the sum of the diagonal elements of BA. Hint. The sum of the 
diagonal elements of either AD or DA is the sum of all products obtained by multi- 
plying any element of A by the associated element of D where we understand that 
the element of D which is associated with n/ is b q p (note the interchange of row 
and column ; for example, the element of B which is associated with the clement in. 
th'> second row nnd third column of A lies in the third row and second column of 
D). Note. Wc term the sum of the diagonal elements of any square matrix the 
trace of the matrix, and we denote the trace of A by Tr A. The result of the pres- 
ent exercise appears then ns follows: 

Tr All = Tr BA 


(although All and BA are in general different and, indeed, in general of different 
dimensions). 

7. Verify the result of Exercise 0 for the two matrices of Exercise 3 and for the 

/ 6 '\ 

two nmtriees A «* (aj, a», 03) and D ~ ( Ir ). Wlmt are the dimensions of 

V>7 

AD and DA ? Note that AD is the product of the column vector B by the row 
V'-ctor A . 

B. Show that matrix multiplication is associative; i.e., if .-1 is any hi X n matrix, 
D any n X p matrix, and C any p X q matrix then both the products (AIl)C and 
A(DC) may be formed and they turn out to be the same. Hint. The element in 

n p 

the rth row and zth column of either product is the double sum n af bf c,’. 

»"1 J“1 

A of/-. Wo denote the common product simply by ABC. 

0. Show that matrix multiplication is distributive with respect to addition, 
i.e., that (A 4- B)C. •» .-tr 4- DC. S’otc. Addition of matrices is defined only 
for matrices of the same type, i.e., of the same number of rows and columns, anil 
the rule of addition may be phrased as follows: Respect the roir and column position. 
Thus if A and D are each m \ n matrices the element in the pth row nnd gth col- 
umn of A 4- D is (by definition! 

«? P 4- 

10. Phow that addition of matrices is n commutative operation (what does this 
mean?). 

11. Phtiw that .1 -f A «* 2.1, 

12. Show that Tr (.4 4- D) - Tr A 4- Tr D. 

v 13, * v hftw that the trn-.iparr of A 4- D is the mm of the transposes 0 f A and II. 
•" U - 5 , f A sh - v r * '< n matrix # c term the matrix obtained by interchanging the 
rvtu r.tsd columns of .! the t-mupo.-e of .4 and w.- denote the transpose of .4 bv an 
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attached star A * Thus the present exercise may be written as follows (A -f- B) * 
= A* + B* The row vectora of A * are the column vectors of A, and the column 
x ectors of A* are the row vectors of A 

14 State in words and prove that (A*)* *= A 

15 Form the cofactor matrix of each of the matrices ^ ^ Note 

If t is any square matrix of dimension n we understand by the cofactor matrix of A 
(denoted by co A) the matrix obtained by replacing each element of A by its 
cofactor the cofactor of o p * being the product of the determinant of the (n — 1)- 
dunensional matrix obtained by erasmg the gth row and pth column of A by 
(_!)!*-« 

(J 'I) (-1 

16 If A is a square matrix of dimension n show that A(co A)* *» (det A) £■ 

where det A denotes the determinant of A Note We Bhall not stop for a defini- 
tion and discussion of the elementary concept of the determinant of a square 
matrix Note, however that determinants may be dejined, by induction from 
n — 1 to n by means of the result of this exercise, it being understood that the 
determinant of a one-dimensional matrix is the (one and only) element of the 
matrix Thus to define det A where A is the two-dimensional matrix ^ ^ "C 
have, first, co A =- (see Exercise 16) Hence (co A)* ■» ^ 

and so A(eo A)* ~ (^ d ” be ^ ® = (ad - bc)Ei Hence det A - ad - te 

/ 2 4 5\ 

17. Calculate co A and det A for A = ( 0 - 1 — 1 1 

\-3 -2 0/ 

/ 2 J -3\ 

Anewer co A ** I 10 15 —16 I det A «• 1 

\ l 2 -2/ 

18 Show that if A is a square matrix of dimension n whose determinant is not 
rero the matrix B ■« (det A) -1 (co A)* has the property that AD — DA " E» 
Note We term the matrix D defined in this way the reciprocal of A and we denote 
it by the symbol A~ l A square matrix whose determinant is not zero is termed 
non etngular 

19 Show that if AB ■* E n , where A is a squarp matrix of dimension n, then A 
is non-emgular and B = A~ l Hml A is non-singular Bince det (AB) * det E* 
“ 1 and det ( IB) ~ det A det B On multiplying the relation AD «* E% 
through on the left by A ~ 1 and using the associative property of matrix multiplica- 
tion wo obtain B — A - 1 

20 Show that if BA “ E n , where A is a square matrix of dimension n, then A 
is non-singular and B ■ A“* Note The results of Exercises 19 and 20 may be 
briefly expressed as follows Only non-singular matrices A can satisfy either of the 
relations AD — h „ BA •« £*, and when A is non-singular the associated matrix 
B is unambiguously determinate, being A -1 
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the element of area in the plane 

21. Show that if A nnd B arc non-singular matrices of dimension n so also is AB 

and determine (/IB) - '. Hint. = E n . 

Answer. (AJ3) -1 *=* B~ , A~ 1 . 

Xott, The important result of this exercise may be phrased ns follows: The 
reciprocal of the product of tiro non-singular matrices of dimension n is the product of 
their reciprocals in the reverse order, 

22. Generalize the result of Exercise 21 to the product of any number of non- 
singular matrices of dimension n. 

23. Show that if A is any m X n matrix and B any n X p matrix, so that the 
product AB can be formed, then the product B*A * can be formed and B*A* =* 
(AB)*. 

21. Generalize the result of Exercise 23 to the product of any number of matrices. 
The important result of this exercise may be phrased as follows: The trans- 
parr nf the product of any number of matrices is the product of their transposes in the 
revert' order. 

23. Show that if A is any square matrix then Tr A* = Tr A. Hint. The 
diagonal elements of A* are the same ns the diagonal element a of A. 

2. The matrix element of arc and the element of area in the 
plane 

Let the reel angular Cartesian coordinates (x, y ) of a variable point P 
in a plane be functions of a single independent, variable, or parameter, 
a. Then P traces, as « varies, a curve (which reduces to a point if 
r and y are each constant functions of a). We term the 2X1 matrix 

the matrix element of arc of the curve that is traced by P. The squared 
magnitude of the. vector v(dx, dy) is furnished by the formula 

(ds)~ — (dx) • -f- (dy)' — (dx) * dx 

whore the 1X2 matrix (dx)* = (dx, dy) is the transpose of the 2X1 
matrix dx. We term ds, i.e., the positive square root of (dx)* dx, the 
scalar clement of arc of the curve traced by P. 

Let, now, the rectangular Cartesian coordinates (x, y) of P be 
functions of two independent variables, or parameters, a nnd /3. We 
have now two 2 X 1 matrix elements of arc d a x and d$x where 



and we set up the 2X2 matrix whose column vectors are d a x and 
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The absolute value of the determinant of this matrix is the element of 
area in the plane, it being understood that this determinant is not 

zero, 1 e , that det ^ ^ 0 Since an interchange of the param 

p) 

eters a and 0 changes the sign of det ^ we may choose these 

PJ 

parameters in such an order that det ~ > 0 and we suppose 

this done Denoting the element of area in the plane by dS, we have, 
then, 

dS “ det ^ da dp 
(«, &) 


It frequently happens that x and y are not given directly as functions 
of the independent variables a and 0 but rather as functions of two 
other variables, which are functions of a and 0 For example, in the 
theory of plane deformations a point P a whose coordinates relative to 
any convenient rectangular Cartesian reference frame are (a, b) is 
deformed into a point P x whose coordinates relative to any convenient 
rectangular Cartesian reference frame (not necessarily the same as 
before) are (x, y) We have now two elements of area, an initial ele 
ment of area dS„ defined by 

dS„ = dot j J da dp 
(«, p) 

and a final element of area dS x defined by 

dS x = det da dp 
(« P) 

It follows immediately from the rule of multiplication of determinants 
(in accordance with which the determinant of the product of tv o square 
matrices is the product of the determinants of these matrices) that 

dS x = dot dS a 
(a, b) 

{prove thi3 Hint r ~7 ) The easiest and quich- 

\ (a, 0) (a, o) (a 0) / 

est way of evaluating dS x is to proceed according to the following 
formal procedure Write down dz = x a da + Xb db, dy = y a da + db 
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and multiply dy 1 >y dx where it is understood that the differentials da 
and db obey the following alternating rules of multiplication: 

da da = 0; db da — —da db; db db = 0. 

Wo find that dx dy = (ar „?/ t - x h y a )da db, and on setting da db = dS a 
we see that rfiS' x is dx dy. 

Example 1. Consider the plane shear defined by the formulas 

x = a -f- kb, 

where k is a constant. We have dx — da + k db, dy = db , and so 
d.S'. = dx dy = da db — dS a . 

Thus the final element of area dS r is the same as the initial element of 
area dS a . 

Example 2. Let the independent variables (a, /3) be plane polar 
coordinates (r, 0). From the equations x = r cos 0, y - r sin e we 
read off dx — cos 0 dr — r sin 0 dO, dy = sin 0 dr + r cos 0 do, and on 
evaluating dx dy (where dr and dO obey the alternating rules of multi- 
plication, dr dr — 0, do dr — —dr <10, dO dO - 0) we obtain dxdy — 
r(cos i! 0 ~r sin 2 0) dr dO - r dr dO. In words, the element of area in 
plane polar coordinates is rdrdO. 


EXERCISES 

1. Show that, in tin* plane shear defined by x *= a + kb, p — b, vectors parallel 
to the ci-avi*' are unchanged whereas vectors parallel to the 6-axis arc turned 
through the angle are tan k and are increased in magnitude by the factor (1 -f- 
Hint, dx ■» da -j- >■ dh, dp « db. When db ~ 0, v(dr, dp) — v(da, dh), and, when 
iffl •" n, v(dr, dvj « r{l: db, rib). 

1. Determine the effect of the plane shear of Exercise 1 on a vector that is 
inclined at an angle of r /-. J to tlie c-axL. 

.1. Show that in tlm plane torsion defined by 

x “ r ens {(? k) n a cos k — h sin k 

< k constant, 

>, •' r sin (ft 4- I.) a sin k -f f, cos 

the final element of art a HS t i^ the s ame as the initial element of area dS a . Show 
nho that the magnitude of any vector in the plane is unaffected by the plane 

ViJi 


(1. Matrix demon tv of are and of surface area and the element of 
xojume in space 


l.’f t the rectangular Onrt'Man coordinates (j 
ip. spne** he function*- of a single independent 


?/, r) of a variable point P 
variable, or parameter. «. 
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Then P traces, as a vanes, a curve (which reduces to a point if x, y, 
and z are each constant functions of a) We term the 3 X 1 matnx 

© / x a da\ 

- I y a da 1 
\z„ da J 

the matnx element of arc of the curve that is traced by P The 
squared magnitude of the vector v(dx, dy, dz) is furnished by the 
formula 

(da) 2 - (dx)* + (dy) 2 + (dz) 2 = (dx)* dx 

where the 1X3 matnx (dx) * = (dx, dy, dz) is the transpose of the 
3X1 matnx dx We term ds (i e , the positive Bquare root of (dx)* dx) 
the scalar dement of arc of the curve traced by P 

Let, now, the rectangular Cartesian coordinates (x, y, z) of P be 
functions of two independent vanables, or parameters, a and p We 
have now two 3X1 matnx elements of arc d a x and d f x where 

( x a da\ / Xff dp\ 

Va da ), dfiX = ( yp dp J, 
z„ da) \zp d/3/ 

and we set up the 3X2 matnx whose column vectors are d a x and 
dgx 

( x a da Xu d( 3\ 

Va da yp dp J 

Za da Zp dp) 

If we regard this 3X2 matrix as obtained by erasing the first column 
of a 3 X 3 matnx the first column of the cofactor matnx of this 3X3 

matrix is ( dS* ) w here dS x = det da dp, dSP — det da dp, 

\dS'/ (a, P) (a, P) 

dS * «• det [ X ’ da dp We term this 3X1 matnx the matnx ele- 

(a, P) 

merit of area of the surface traced out by P The vector v(dS x , dS*, dS‘) 
is the vector, or cross, product d„x X d^x, and so v(dS x , dS”, dS') ha3 
the direction of one of the two normal directions at P to the surface 
traced out by P (it being understood that d„x X d,jx is not the zero 
vector) In other words, the column vector of the matrix element of 
area is normal to the surface If dS x denotes the matnx element 
of area the scalar clement of area dS x is the positive square root of 
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(dS*)* dS* where ( dS x )*, the transpose of dS x , is the 1X3 matrix 
(d,S x , dS' J , dS‘). 

It frequently happens that x, y, and z are not given directly as 
functions of the independent variables a and £ but arc given rather as 
functions of three other variables, which are functions of a and p. For 
example, in the theory of space deformations a point P n whose coordi- 
nates relative to any convenient rectangular Cartesian reference frame 
are (a, b, c) is deformed into a point P x whose coordinates relative to 
any convenient rectangular Cartesian reference frame (not necessarily 
the same as before) are (x, y, z). We now have two matrix elements 
of area: an initial matrix element, of area, 

I dot ~~ da dp 

[ («, P) 

; del - -- -- da dp , 

. («, p) 

1 dot 7 -^— da dp 
and a final matrix element of area. 


dS~ = 


Adopting the alternating rule of multiplication of differentials (in 
accordance with which da da ~ 0, dp da = -da dp, dp dp = 0), we 
may write dS 3 in the. more compact form 

fdb t/c\ 
dS a «= I dc da |, 

\da dbj 

and v«e may write dS- in the more compact form 


/ det 



- da dp ^ 


(«, P) 

det da dp 

dot da dp 

(«, P) 





'dy dz\ 
d: dx j. 
K dx dy) 
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Since y a = y a a„ + y b b a + yefia, etc , we have 


det 


(y, z) 

(«, &) " 


(b, c) (tr, 0) (c, a) (a, 0) 


, A f Z ) ^ * ( g > b > 4 

+ det - — — det - — — t etc , 
(o, b ) (a, 0) 


(prove this), and so 

, (y , «) 


dS* = det Sfli' dS“ + det dS b 
(b, c ) (c, o) 


+ det ^ rfs- 

(«. 6 ) 


Jacobian matrix ; 


This and the two similar equations, which furnish dS u and dS* (write 
down these equations), tell us that the 3X1 matrix dS* is the product 
of the 3X1 matrix dS a by the 3X3 matrix co J where J is the 
(z, V, «) 

(o b, c) 

dS * = (co J) dS a 

Taking it for granted that J is non singular (i e , that det J j* 0), we 
may write this relation in the equivalent form 

dS x «= (det JKJ*)-' dS a 


Note that det J = det J * and that (J*)~ 


Vdet J*j 


co J 

(y , *> , 


It 


- det dS a + det dS b -f- 
(6, c) ( c , a) 


is clear from the relations dS x 

det ~r^ dS c , etc , that dS x may be calculated by multiplying out 
(a, b) 

dy dz = (y a da y b db + y c dc)(z a da + z b db + z c dc ), subjecting the 
differentials da, db, dc to the following alternating rules of multiplica- 
tion da da — 0 = db db ~ dede, dedb — —dbdc, dadc = — deda, 
db da — — da db, and finally setting db dc = dS a , dc da — dS b , da db = 
dS c 

Example 1. Determine the matrix element of area of the sphere 
x = r sin $ cos 4>, y = r sin 0 sin <£, z = r cos 0, r constant We have 
dy — r cos 6 sin <f> dO + r sin 0 cos <t> d<f>, dz = — r sin 6 dd, and so 
dS x = dy dz — ~r* sin 2 0 cos 4> dif> d0 — r 2 sin 2 0 cos dd d<t> Simi- 
larly, (show this) dS? — dz dx = r 2 sin 2 0 sin <j> dd d<t> and dS‘ = 
r 2 sin 9 cos 0 dd d 4 >, and thus the 3X1 matrix element of area is 


dS x 


( sin 8 cos <A 
sin 9 sin <f> J 
cos 6 J 


d9 d<f> 
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1 fence (show thin) the scalar element of area dS~ — r 2 sin 6 dO d<f>, and 
the column vector of the matrix element of area 1ms the direction of 
the outward drawn radius. 

Example 2. Determine the matrix elements of area of the surface 
into which the surface of a rectangular box whose faces are parallel to 
the various coordinate planes is deformed by the shear x = a + kc, 
y « //, z — c, it being understood that t is a constant and that the 
point P„:(n, b, c) is deformed into the point Pt‘(x, y, z). 

We have dx — da -f k dc, dy = db, dz — dc, and so 

(dy d:\ (db dc \ / dS a \ 

dS~ — ( dz dx j = ( dc da j — ( dS h J. 

\dx dy] \da db + 1: dc db] \dS c - k dS°] 


Tims the matrix element of area of a face of the box which is parallel 
to the rt-eoordinate plane (for which d,S h = 0 = dS c ) is changed by the 


shear from 



to dS a 



The scalar clement of area is 


magnified by the factor (1 + /r) ! \ The. matrix clement of area of a 
face of the box which is parallel to the 5-coordinatc plane (for which 
dS n (1 - r/,S v ) is unchanged by the shear (prove this), and, similarly, 
the matrix element of area of a face of the box which is parallel to the 
c-coordinatc plane (for which dS a — 0 — dS b ) is unchanged by the 
shear (prove this). 

Let, now, the rectangular Cartesian coordinates (x, y, z) of P be 
functions of three independent variables, or parameters, a, p, and 7. 
We have now three 3 X 1 matrix elements of arc d„x, d$x, and dyX 
where 


(x a dci\ (x}dd\ ( x y dy\ 

d„x - I y„ da J, d.-T = ( y t dp J, d r r. = ( y 7 dy J, 
V„ da/ \zs dp/ \z y dy/ 


and we set tip the 3X3 matrix whose column vectors are d„x, dsx, 
and <hx: 

( x„ dn X; dp x y dy\ 
y a da y. dP it, dy j. 
z„ da dp Zy dy] 


Tim 

if! to 

t.e.. 


absolute value of the determinant of this matrix is the element of 
' 1 •’ in ‘Triee. it being understood that this determinant is not zero, 
. , (x, v. r) 

mat net 0 . Vo can onier the independent variables 
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so that det j- 1 > 0 (show this), thus, on denoting the element of 

(a, ft V) 

volume by dV x , we have 

dV x = det Z \ da dp dy 
(«, ft 7 ) 


It frequency happens that x, y, anti z ate not given ihrectiy as inactions 
of the independent variables a, p, and y but rather as functions of 
three other variables, which are functions of a, ft and 7 Tor example, 
in the theory of the deformation of a three-dimensional medium, a 
point P a whose coordinates relative to any convenient rectangular 
Cartesian reference frame are (a, b, c) is deformed into a point P* 
whose coordinates relative to any convenient rectangular Cartesian 
reference frame (not necessarily the same as before) are (z, y, z) 
We have now two elements of volume an initial element of volume 
dV a , defined by 

dV a — det " da dp dy, 

(“> ft y) 


and a final element of volume dV x , defined by 

dV x — det ^ da dp dy 
(a, ft 7) 

(it being understood that the independent variables a, ft 7 have been 
so ordered that dV„ and dV z , defined in this way, are positive) It 
follows immediately from the rule of multiplication of determinants 
that 

dV « = det ft Z { dV„ - (det J ) iV „ 

(a. b, c) 


where J is the Jacobian matrix of (x, y, z) with respect to (a, b, c) 
(prove this) An easy way to evaluate dV x is to multiply out 
dxdydz =* (x a da + x b db + x c dc){y a da + y b db + y c dc){z a da + 
z b db + z c dc) (it being understood that the differentials da, db, dc obey 
the alternating rules of multiplication da da da = 0 = da da db — 
etc, da dc db = —da db dc, etc ) and to set da db dc = dV a If p a and 
Px denote the initial and final densities of the medium that is being 
deformed the principle of conservation of mass tells us that 

Px dVx = Pa dV a , 

and so the relation dV x — (det J) dV a may be written in the equiva- 
lent form 

p a — px det J 
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EXERCISES 

1. Phow that the final 3 X 1 matrix element of area dS : is connected with the 
initial 3 X 1 matrix element of urea dS a by the formula 

dS- » - (,/*)-’ r/.S°. 

P: 

2. Show that volumes arc jjrescrvcd under the shear x = a + he, y *> h, z « c. 
flint, dr dtj dr «* do dh dc. 

3. Show that volumes arc preserved under the torsion x ® r cos (0 + £c) = 
a eos Jr — 6 «in ;. fl y « r (-in (0 + J.r) *» a sin Arc -f b cos J:c, z «= c, k constant. 

t. Show that the element of volume in cylindrical coordinates (r, 0, z), where 
f ». r cos 0, y » r sin 0, z « ; is r dr d 0 dr. 

3. Show that the element of volume in space polar coordinates (r, 0, <£), where 
- « r sin Dean tfi, y « r sin 0 sin <^>, ; ■=■ r cos 0 is r s sin 0 dr d 0 d <f>. 


t. The canonical form of a linear operator 

The discussion throughout the present section will concern a square 
matrix of dimension 3 but the argument is applicable to square matrices 
of any dimension (repeat the discussion in full detail for a square matrix 

(u v. w) 

of dimension 2). We have seen that the Jacobian matrix J = 

(*, V, s) 

of three differentiable functions (v, v, w) of (x, y, z) may be regarded as 
n magnification factor that converts, by multiplication, the 3X1 
matrix dx into the 3X1 matrix du: 

die — J dx. 


Similarly, any 3X3 matrix may be regarded as a magnification 
factor that converts, by multiplication, a 3 X 1 matrix a into a 3 X 1 
matrix b: 

b — X fa. 

us now consider the space vectors a and b whose coordinates 
relative to any convenient rectangular Cartesian reference frame are 
the elements of a and b, respectively: 

a “ r(o 1 , a", a 3 ), b = v(b\ fr, 6 3 ). 

l.f-t the coordinates of a and b relative to a second rectangular Car-* 
w-Mun reference frame (which may be obtained from the first coordinate 
reference frame by a rigid rotation) be furnished by the 3 X 1 matrices 
a' and //, respectively. (A 'ole. A prime is often used to indicate the 
transpose of a matrix but we have adopted a star to indicate the trans- 
pose and this leaves the prime free for the present purpose; a ' is not 
the transpose of a; in fact, n’ is, like a, a 3 X 1 matrix, whereas the 
nn t -jio-v< n* of o b a i X 3 matrix.) Let the coordinates of unit 
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vectors along the axes of the new reference frame relative to the original 
reference frame be furnished by the three 3X1 matrices Ui, « 2 , and 
U 3 , and denote bj R the 3X3 matrix whose column vectors are 
U 2 , and U 3 The fact that the three unit vectors along the axes of 
the new reference frame are not only of unit magnitude but are also 
mutually perpendicular is conveniently expressed by the matrix 
equation 

R*R - E z 

(prove the validity of this equation) The fact that det R* =* det R 
assures us (why?) that det R = ±1, and the fact that the second 
reference frame may be obtained from the first by a rigid rotation 
restricts us to the relation det R = 1 (why?) ( Hint det R is a 
continuous function of the elements of R, hence it can assume only 
one or the other of the values ±1 Since it is originally, when R = 
E z , +1, it must be permanently +1 ) Thus R is a 3 X 3 matrix 
whose determinant is +1 and whose reciprocal is its transpose We 
term any such 3X3 matrix a three dimensional rotation matrix If 

M 

c - ( c 2 } is any 3X1 matrix we have Rc = c l ui + c 2 u 2 + c 3 «j 

w 

This equation tells us that the coordinates m the original reference 
frame of the vector, w hose coordinates m the new reference frame are 
furnished by the elements of c, are provided by the elements of Rc 
On setting c ~ R*a and using the fact that RR* = Ez, we have the 
following result 

The coordinates in the original reference frame of the vector, whose 
coordinates in the new reference frame are furnished by the element of 
R*a, are provided by the elements of a In symbols 

a' = R*a 

or, equivalently (why?), 

a = Ra’ 

Similarly, b = Rb’, and so the relation b = Ma yields RV — MTta' 
or, equivalently, 

V = MW 

where M' = R*MR We say that the elements of M furnish the 
coordinates ra the original reference frame of a linear vector operator 
M and that the coordinates of this same linear vector operator ra the 
new reference frame are furnished by the elements of M' = R*MR 
(just as the coordinates of the vector a in the original reference frame 
are furnished by the elements of the 3X1 matrix a whereas the coor- 
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and b — t>(6 1 , i» 2 , b 3 ) are said to be orthogonal (or perpendicular) when 
b*a — 0 (Show that the orthogonality of two complex vectors is a 
commutative relation, 1 e , that b*a = 0 if, and only if, a*b = 0 
Hint b*a = (a*b)*) A 3 X 3 complex matrix each of whose column 
vectors is of unit magnitude, these unit vectors being mutually per- 
pendicular, is termed a unitary matrix, and we denote a unitary matrix 
by the symbol U 


EXERCISES 


1 Show that a 3 X 3 matrix U is unitary if, and only if, V*U = Et Note 
The result of this exercise tells us that a 3 X 3 matrix is unitary if and only if, its 
reciprocal is its star, i e , the conjugate of ita transpose 

2 Show that a 3 X 3 matrix U is unitary if, and only if, UU* =** 

3 Show that every rotation matrix R is unitary 

4. Show that if U js unitary det U is a complex number of unit modulus i e , 
that det V — e {> , e real Hint det U* is the complex conjugate of det U 

5 Show that if each of the column vectors of a unitary matrix is multiplied by a 
complex number of unit modulus (not necessarily the same for the various column 
vectors) the resulting matrix is still unitary 

6 Show how to construct a unitary matrix whose first column vector is an 
arbitrary given unit vector Hint The 6econd column vector Ui is any unit 
vector perpendicular to the given unit vector ui, and the third column vector is a 
unit vector whose coordinates are proportional to the conjugates of the coordinates 
of (Hi X u 2 ) 

7. Show how to construct a untmodular unitary matrix, i e , a unitary matrix 
whose determinant is unity, whose first column vector is an arbitrary given unit 
vector Hint Combine Exercises 6 and 5 

8 Show that the general two-dimensional unitary matrix is of the form 



where B is real and fia + Sb « 1 

9 Show that the general ummodular two-dimensional unitary matrix 
form 


(: 1) 


where 6a + Sb = 1 


of the 


10 Show that the product of two unitary matrices is unitary Hint 


{U\Ut)*UiUs = Ut*Ui*Ui(Jt =» Ui*U t = E t 

11. Verify that the product of two matrices of the type where 6a + 

Sb “ 1, is again of this type 

Ummodular unitary matrices U play m the complex field the role 
played by rotation matrices R m the real field Thus, if a is any 3X1 
complex matrix whose elements furnish the coordinates of a three- 
dimensional vector a, the coordinates of this same vector in a second 
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complex rectangular Cartesian reference frame are furnished by the 
matrix o' *= U*a where the column vectors of U furnish the coordi- 
nates, relative to the first coordinate reference frame, of unit vectors 
along the axes of the second coordinate reference frame. If M is a 
linear vector operator that converts, by multiplication, the complex 
vector a into a complex vector b, 

b = Ma, 

we have the two equations 

b = Ma, V = MW 

where the elements of the two 3X3 matrices M, M' furnish the 
coordinates of the linear vector operator M in the first and second com- 
plex rectangular Cartesian reference frames, respectively. The con- 
nection between the two different presentations of the same linear 
vector operator M is furnished by cither of the two equivalent formulas 

iff « U*M'U, M' = UMU* 


(prove the equivalence of these two formulas). 

Wc arc now prepared to attack the problem of determining a particu- 
lar rectangular Cartesian reference frame (in general, complex) in 
which the presentation M' of a linear vector operator M (whose pre- 
sentation in a given rectangular Cartesian reference frame, usually 
real, is a given 3X3 matrix M) is particularly simple. Let X be a 
runt of the cubic equation det. (M — \E 3 ) = 0. Then there exists a 
non-zero 3X1 matrix a that satisfies the equation jlfa = Xa (why?). 
If X happens to be real and if the elements of M are real the elements of 
a may he taken to be real, but if X is a non-real complex number the 
elements of a will not all be real when the elements of i\f are real (for 
if the elements of a were real the elements of Ma would be real and we 
could not have Ma *= Xa since X is, by hypothesis, non-real). In any 
event let n, be a 3 X 1 matrix, whose column vector uj is of unit 
magnitude, which is proportional to a (i.e., «i = l:a where k is some 
complex number). Construct a 3 X 3 unimodular unitary matrix 
f - 1 wIhko first column vector is u\. Then, on denoting bv cj the 3X1 

/A 


matrix 


0 

, 0 , 


I, we have the two relations 


«= Aw 1, Wj - U1C1 


(show this), and these imply the relation MU,ci - XDVi or, equiva- 
lently. 


IV-VfV, - Xfi. 
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In other words, the first column of U\*MU\ is [ 0 1 Thus the 

w 

complex rectangular Cartesian reference frame determined by Ui is 
such that the presentation, in it of the linear vector operator M. has 
zeros in the second and third rows of the first column We may write 
this presentation of M in the form 


(:'»') 


where N is a two-dimensional matrix and the stars denote numbers 
(complex) about whose \alues we can say nothing further Applying 
to N the argument just presented concerning If, we construct a two- 
dimensional ummodular unitary matnx V which is such that 


The matrix 172 = 


/I 0 OX 

= 10 „ J is a three-dimensic 

\0 / 

unitary matnx (prove this) which is such that 

’ - )""G : 3 

1 it follows that 

U 2 *Ui *MUiUz = [0 fi * ) 

\0 0 v] 


three-dimensional ummodular 


(prove this), and it follows that 


Ui*Ui*MUiUz = 


The matrix U 1 U 2 is a three-dimensional ummodialar unitary matnx 
that we denote simply by U We have, then, the following funda 
men tall} important result 

If M is any three-dimensional square matnx, there exists a three 
dimensional ummodular unitary matnx U which is such that U*MU 
has zeros below the diagonal We term the matnx M' = U*MU, 
which has zeros below the diagonal, the canonical form of the linear 
vector operator w hose coordinates, in the original reference frame, are 
furnished hy the elements of If 
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If .1/ is Midi that M * is the panic as .V we term it. symmetric if its 
element « are real and Hermit inn if its elements are complex. It is 
clear that M' - U*MV satisfies the same relation (.IP)* = M' 
(prove this) ami so .1/' is diagonal with real diagonal elements. If 
the elements of ,1/ are real it follows (prove this) that the elements of 
P may lie taken to he real (so that U is a rotation matrix). We have, 
then, the following important special ease of our principal result: 

If M i*. any three-dimensional symmetric matrix there exists a rota- 
tion matrix It such that Ii*MIt is diagonal with real diagonal elements: 

fin i 0 0 \ 

It* MR = 0 wi 2 0 • 

\0 0 vi J 

The three real numbers (?»i, »i 2 , in 3 ) are known as the characteristic 
numbers of the linear vector operator M whose coordinates in the 
original reference frame are furnished by the elements of the real sym- 
metric matrix M . 
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real complex number x*Mx is never negative (no matter what is the 
complex vector x) we say that the linear vector operator M, whose 
coordinates are furnished by the elements of Af, is non-negative It is 
clear that M is non negative if, and only if, its characteristic numbers 
are non-negative (Prove this Hint = mix 1 x l -f 

mzi^x 2 ' + m, 3 X z x v ) We now construct the 3X3 matrix 



where mi 44 denotes the unambiguously determinate non negative 
square root of mi, etc , and it is clear that the square of this matrix is 
the canonical form M' of M We denote it, therefore, by (A/') 44 , 
and it is clear that (A/') 44 commutes with AT, 1 e , that = 

A/ 7 (A/') 44 , since both M' and (A/') H are diagonal matrices On setting 

Jf M = U(M'> M l/* 

(so that the elements of A/ 44 are the coordinates m our onginal refer- 
ence frame of the linear vector operator M 44 whose coordinates, m the 
reference frame in which the coordinates of M are furnished by the 
elements of M', are furnished by the elements of (AT ) M ), it is clear 
that M 44 is a non negative linear vector operator (why?) and that the 
square of M 44 is M (prove this) It is also clear that M H commutes 
with M, 1 e , that M m M = MM 44 (prove this) 

If A is any 3X3 matrix it is clear that the 3X3 matrix M defined 
by M = A* A furnishes the coordinates of a non negative linear vector 
operator In fact, M* = M (prove this) and x*Mx = x*A*Ax = 
(Ax)* Ax is never negative, being the squared magnitude of the vector 
whose coordinates are furnished by the elements of the 3X1 matrix 
Ax If the matrix A is non singular, M is not only non-negative but 
also positive, 1 e , x*Mx is positive save when x is the zero 3X1 
matrix (m which event x*Mx is evidently zero) In this case M 
is non singular (why?), and this implies that A/* 4 is non singular 
(why?) Denoting the reciprocal of A/ 44 by A/ -44 , we can easily see 
that the matrix AAf -44 is unitary In fact, (A A/ -44 ) *(AM~ V> ) = 
A/ -44 A*AA/- h = = (A/ _44 ) 2 A/ = M~*M = E% On 

denoting the unitary matrix AM~ 44 by U, we have A — UM 44 We 
obtain, in this w ay, the following fundamental result 

Every non singular 3X3 matrix A whose elements are complex 
numbers may be written in the form 


A = UM l 
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v/liorc U is n unitary 3X3 matrix and is a positive (and, hence, 
Ilermitinn) 3 X 3 matrix. 

If the elements of A are real. M = A* A may bo transformed into its 
canonical form by a rotation matrix (why?), and so the elements of 
M~ Xi arc real (why?). Hence the elements of U = AM~ Xi are real. 
If dot A > 0 it follows (since det M~ v ' > 0) tiiat det U > 0 and so 
V is a rotation matrix. Wc have, then, the following particular 
instance of our principal result: 

If A is a non-singular 3 X 3 matrix, with real elements, whose 
determinant is positive, A may be written in the form 

A = RM ]i 

where It is a rotation matrix and M — A* A; A / 54 is a symmetric 
matrix wlio^e characteristic numbers wi !i , 7H» 54 , 3 ^ are all positive 

(mi, ?ii}, vh being flic positive characteristic numbers of M = A *.4). 
This result is of fundamental importance in the theory of elasticity. 

6. Conclusion 

If one has mastered fairly well the contents of this introductory 
chapter lie is in possession of all the matrix theory required in this 
treatment, of the finite deformations of an elastic solid. For an intro- 
ductory study at a level suited for senior high school or beginning 
college students the reader is referred to the first five chapters of 
Analytic Geometry, 1 where analytic geometry is treated from the very 
beginning from the point of view of vectors and matrices. More 
advanced matrix theory is presented in the first two chapters of 
Introduction to Applied Mathematics 3 and a treatment of those parts 
of the theory of matrices which are of importance in group theory is 
given in Theory of Group Representations. 3 Try to become familiar 
with matrices and to use them regularly as a helpful tool. Do not 
fall into the error of regarding them as a complicated device invented 
by mathematicians to make the theory of elasticity harder than it 
actually is. 
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2. Write down the most general 3X3 matrix that commutes with the diagonal 
matrix 


where dt ^ d i 
Answer A = 


(dt 0 0\ 

D - { 0 dr 0 I 

\ 0 0 dt) 


we can say nothing 


0 S 5 


where the stars indicate elements about whose values 


(ar, y, z) 

3 Calculate the matrix M = ,4 M where A is the Jacobian matrix-- — r — of the 

(a, b, c) 


shear i «= a + Are, y ** b z >=■ c (A: constant) 

4 Repeat Exercise 3 for the torsion x = r cos (8 -f Arc) « o cos Are — b sin Arc, 
y => r sin (0 + kc) — asm Are + 6 cos Ac, z = c {k constant) 

5 Determine the axis and angle of the rotation whose matrix is 



Answer The line from (0, 0, 0) to (1, 1, 1), 
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1. The strain matrix 


Wc consider two states of a three-dimensional deformable medium: 

(1) The initial or unstrained state. Wc denote the rectangular 
Cartesian coordinates of a typical particle of the medium, when in 
the unstrained state, with respect to any convenient reference frame, 
by (a, b, r). 

(2) The final or strained or deformed state. We denote the rectangu- 
lar Cartesian coordinates of the particle, whose coordinates when the 
medium was in the initial unstrained state were (a, b, c), by (x, y, z). 
The reference frame with respect to which (x, y, z) are measured need 
not be the same (although it is usually convenient to have it the same) 
as the reference frame with respect to which (a, b, c) are measured. 

If we consider a collection of particles occupying, when the medium is 
unstrained, a certain volume l'„, a, b, and c will each be functions of 


three independent variables (a, 13, y) which are such that det 


(a, b, c ) 


(a, jS, 7 ) 

0. Wo suppose the independent variables, or parameters, «, f, and 

> so ordered that det •— — > 0, and wc denote bv dV a the product 
(o, a. y) 


det 


(«. b, c) 

(«, 8, y) 


da d(i dy. 


dV'a is the initial element of volume, (i.e., the 


element of volume of the medium when it is in its unstrained state), 
and we may, by adopting the alternating rule of multiplication (what is 
this?) of the differentials da, d,3, dy, write it in the form da db dc : 


dr„ — dadbdc t= det 


(o. b. c) 
(or, d. y) 


da d8 dy. 


1 iir- coordinates (r, y, ;) of the final position of the particle whose 
initial coordinates are (a, b, c) are assumed to be differentiable func- 
tions of (a, h. c), and it is assumed further that the determinant of the 
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Jacobian matrix 7 1 \ is positive Then the final element of volume 

(a, 6 , c) 

is given by the formula 

dV x = dx dy dz = det ^ dF„ = (det J) dV, 


(a, b, c ) 


where J denotes the Jacobian matrix 


fo V, z) 


If p a and p x denote, 


(a, b, c) 

respectively, the initial and final mass densities, the principle of con- 
servation of mass assures us that 

Pz dV z - Pa dV a 


and so the compression ratio — is furnished by the formula 


— — * 
p a det J 

Let us now consider a collection of particles that he initially on a 
curve C 0 Then a, b, and c are functions of a single independent 
variable or parameter a, and the matrix element of arc of C„ is the 
3X1 matrix 



The scalar element of arc ds a of C n is the positive square root of (da)*da 
(ds a ) 2 = {da)* da 

The particles of the deformable medium that originally lay on the 
curve C„ he, after the deformation, on a curve C x , and the matrix 
element of arc of C r is the 3X1 matrix 

/dx\ 

dx = I dy J = J da 
\dzj 

The scalar element of arc ds x of C x is the positive square root of 
{dx)* dx 

(ds x ) 2 = (dx)*dx « ( da)*J*Jda 


If da* is the same as ds a f ot every curve C a (1 e , for an arbitrary 3X1 
matrix da) all curve lengths are unchanged by the "deformation” so 
that, in particular straight lines are “deformed” into straight lines 
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fv.-hy?). This moans that r, y, and z are linear functions of a, b, and 
c (v,"hy?) and so J is a constant matrix. Since this constant matrix 
satisfies the relation 

(da ) *./ *./ da = (da)* da 

for arbitrary choices of the 3 X 1 matrix da we must have J*J = 
E z (prove this). Since det .1 > 0 it follows that ./ is a rotation matrix. 
Thus the transformation from the initial coordinates (a, b. c) of any 
particle of the medium to the final coordinates ( x , y, z) of the same 
((article is merely that involved in the transformation of coordinates 
from one rectangular Cartesian reference frame to another. In other 
words, the change from the initial to the final position of the medium is 
a mere rigid displacement and no deformation or change of shape of 
tlie medium (or of any portion of the medium) is involved. When 
./ is not a rotation matrix we use the difference between J*J and E 3 as 
iv measure of the deformation of the medium in the neighborhood of 
the point (a, b, c). We set 


UJ*J -E a) = 1J 

and we say that the 3 X 3 symmetric matrix y measures the deforma- 
tion of the elastic medium at the point (a, b, c). This matrix y is 
known as the strain matrix, and it satisfies the relation 

(ds T )~ — ( ds „ )" — 2(da)*y da 

(prove this). The use of the factor h in the definition of y or, equiva- 
lently, of the factor 2 in the relation just written, is due to the fact 

that this relation involves the difference of the squares of ds x and ds„. 
If the strain or deformation is small enough so that we can replace 

dr, -!-(/,*„ by 2ds„, (ds,)" — (ds a ) z becomes 2 (ds x — ds a )(ds a ) (why 

cannot wc replace <h x - ds„ by zero?), and we. have the approximate 
relation 



The column vector of the 3 X 1 matrix ~ is a unit vector; if this 
t ds a 

3 X 1 matrix is denoted by », the quadratic form u*yu in the elements 


of v measures, approximately, the relative magnifications — 

ds a 

m the neighborhood of the point (o, b, r). It is dear from the relation 
{<?*,)• - (as fir !!(</<])*•; da that the quadratic form u*yu measures 

. 1 - (dsfir 

fN . lr .|y „ __ (prove this). 
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The moat direct way of calculating rj for any given deformation is 

first to set up J = \ and then to calculate the symmetric matrix 

b, c) 

J*J The diagonal elements of J*J are the squared magnitudes of 
the column actors of J, and each non-diagonal element of J*J is a 
scalar product of two of these column vectors After J*J has been 
evaluated, subtract Ez from it, i e , dimmish each of the diagonal ele- 
ments of J*J bj unity one half the resulting matrix is the desired 
strain matrix 

Example Calculate ij for the shear z = a + kc, y = b, z = c, 
where k is a constant Here 


and so 


/l o k\ /I 0 k \ 

J- 0 1 0 J*J- 0 1 0 ) 

\0 0 1/ \l 0 1+ /by 



If k is so small that its square may be neglected we say that the shear is 
infinitesimal thus the strain matrix for an infinitesimal shear parallel 
to the a axis in the ac-plane equals 



EXERCISES 

1 Calculate the strain matrix 17 for the uniform compression or dilatation * — 
a + ka y »* b -f- A 6 z *• c +kc k a constant Note This is a compression if 
A < 0 and a dilatation if k > 0 

Answer ij ■ (1+ jA^Fj 

2 What is the strain matrix for an infinitesimal uniform compression or dilata 
tion? 

Answer ij * AFj 

3 Calculate the strain matrix >j for the torsion x — r cos (6 + Ac) — a cob kc — 
b sin kc y — r sin ( 9 + Ac) — a sin Ac + 6 cos kc * — e A a constant What 
does 1 ? reduce to when the torsion is infinitesimal? 


( 0 

0 

-Ab \ 

/ 0 

0 

— Ab\ 

Answer i ( 0 

0 

ka ) { 

M 0 

0 

Aa \ 

\-*b 

Aa 

»*(■’ + 6*)/ 

\-Ab 

Aa 

0 / 


4 Calculate the strain matrix 7 for the simple tension and cross-sectional con 
traction of a rod defined by X " a — aka y — b — akb z ■« e + kc where k and a 
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nrc constant?. Note. The constant <r is known tvs Poisson's ratio (after S. D. 
Poison 1 1781— 1 S10], a French mathematician). 

/-rifc +$***» 0 0 \ 

A nutter, v m [ 0 — tr/: + • 5 <r ! /: J 0 r 

V o 0 k+ik-J 

5 , What doe 1 * th»' strain matrix of Exercise -1 reduce to when the strain is infini- 
l primal, i.e., wiien the square of f: may be neglected? 
t-ck 0 0\ 

Artixter. v " ( 0 ~~°k 0 )‘ 

V 0 Ok/ 

2. The behavior of the strain matrix under transformations of 
the initial and final coordinate reference frames 

A rotation of the initial coordinate reference frame, i.e., of the rectangu- 
lar Cartesian reference frame with respect to which the coordinates of a 
typical particle of the clastic medium are (a, b, c), is described by the 
formula 

a — > a' — R*a 


where the coordinates with respect to the old reference frame of unit 
vectors along the axes of the new reference frame are furnished by the 
column vectors of R. We understand by the notation a— > a’ ~ R*a 
that the new initial coordinates of the particle, whose original initial 
coordinates were furnished by the elements of the 3 X 1 matrix a, 
are furnished by the elements of the 3 X 1 matrix R*a. It follows 

(n, b, c) 


from the relation o' = R*a that a — Ra ' and, hence, that 


R. Since 


(x, y, z) (r, y, z) (a, b, c) 


(o', 6', c') 


we have the following result: 


O («, b, c) (o', h', c/) 

Under the transformation a— * o' = R*a of the initial rectangular 
Cartesian reference frame, J —* J' — JR, the final coordinate reference 
frame being unchanged. 


EXERCISES 

1. Show that under the transformation x—» x' <= R’x of the final rectangular 
C'artcrian frame J — * J’ >» R *J. 

2. Show that under the simultaneous transformations a — ► o' «=> R , n, x — ♦ x’ = 
R‘r ct both the initial and the final rectangular Cartesian reference frames J — * 

r - irJR. 


Since M JV and since the transformation J — * J 1 — JR implies 
the transformation J* -> J* f = R*J* (why?), it follows that M — M' 
" Since the strain matrix r? = i(-V - E s ), this fact yields 

the following important result: 
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Under the transformation a—*a' = R*a of the initial rectangular 
Cartesian reference frame n — ► rf <= R*^R 
In other words, the strain matrix jj is sensitive to the initial rectangu- 
lar Cartesian reference frame There exists (at each point of the med 
lum) at least one particular initial reference frame in which ij appears 
in canonical ( = diagonal) form We term the axes of any such par- 
ticular initial reference frame -principal axes of strain and we write the 
strain matrix ij when referred to a reference frame whose axes are 
principal axes of strain, as follows 



Note carefully that, unless ij is a constant matrix, the principal axes 
of strain will vary from point to point of the medium Note further 
that if £ is a rotation matrix that transforms v into its canonical 

/vi 0 o\ At o o\ 

form, so that R*rjR = f 0 0 J we have tjR — R I 0 ij 2 0 1 

\0 0 W \0 0 t), / 

Hence the column vectors of R are characteristic vectors of g (why*) 
The easiest way to determine R is to determine the characteristic 
vectors of r> It frequently happens that one of these is relatively 
easy to determine (the determination of the others being comparatively 
difficult) The best thing to do m this event is to construct a pre 
lirmnary rotation matrix Ri whose first column vector haa the direction 
of the easily determined characteristic vector of »j Then will 

be of the form 



and our problem becomes the comparatively simple one of & plane 
strain whose matrix is indicated by * * 

In contrast with its behavior with respect to a change of the initial 
rectangular Cartesian reference frame the strain matrix v is insensitive 
to a change of the final rectangular Cartesian reference frame In 
fact, under the transformation x— *x' = R*x, J —* J' = R*J, J*—* 
J*' = J*R, M -> M' = J*RR*J = M Hence v -» v’ = v (why?) 
Thus the elements of the strain matrix ij are invariants under trans 
formation of the final rectangular Cartesian reference frame, to 
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emphasize tliif* point we denote these elements by the symbols 

r tn')t CtC." 


> 

Von 

Vo b 

Voe\ 

Via 

Vhb 

%<:] 

flat 

Vcb 

Vre/ 


Note Hmt the first of the two subscripts tells the row and the second the 
column; thus y hr is the element in the second row and third column of y. 

Itanart: /. We use y„ n , Vai, etc., instead of y a a , vi a , etc., since v is 
not the presentation of a linear vector operator but rather the matrix 
of rocflioicnts of a quadratic form u*yu where w is the 3X1 matrix 

element of arc ~~ of C a whose column vector is of unit magnitude. 
( 1 % 

The fact that the strain matrix y is symmetric is expressed by the 
relations 

Vcb 2=1 Vbc* Van ~ Vcai Via Vab* 

Rrmnrl: 2. The i nsemitivencss of y to a change in the final rectangu- 
lar Cartesian reference frame may be expressed by the statement that 
the strain matrix is unaffected if the given deformation is followed by a 
rigid rotation of the medium; thus both ,7 and I{*J yield the same ij. 
The fcnxilirciMS of y to a change in the initial rectangular Cartesian 
reference frame may he expressed by the statement that the strain 
matrix is affected if the given deformation is preceded by a rigid rota- 
tion of the medium. Thus .7 and JR yield different strain matrices. 


EXERCISE 

.1. Show that the strain matrix y transforms under the simultaneous transforma- 
tion a • o' m fl'n, x— • x' «* /f*x of the initial and final rectangular Cartesian 
refert nrr frames according to the formula ij — ♦ ij' *= Hint, ij is insensitive 

to the transformation x' ■* /?*x. 


Example 1: Determine the principal tixcs of the plane shear x = 
a / 7), y h, k constant. Here .7 — ^ and so y — l- 

hetting It r~ ^ (so that the new initial reference frame 

obtained by rotating the old initial reference frame through the 
angle a), wc have 


2rj{ / f: t ’ in « k cos ft \ 

V‘ eos o -f- };- sin a —7; sin « -f- Jr cos aj 

„ ‘•hi - n "t* l" sin* ft J; cos 2<rr -f- id; z sin 2a \ 
V* 2a -'r 7 . Jr sin 2a —k sin 2a -f Jr cos" a ) 


2 R'rji 
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In order that R*qR be diagonal we must have cot 2a =* — -or, 


2 


equivalently, sin 2 a = (o lying in the interval - < « < -i 

if h > 0) On eliminating a, we find that ??j 


I 2 

(4 + kY + T + 


4(4 + &)* " “ (4+1’)“ + T - i(4T jV* ^ Intllem,ta '- 
tesimal theory of plane shear (in which l 2 is neglected) the angle a 


independent of A. To this degree of approximation tji = 


k -fc 

" 5 ”" _ 2 

The next approximation to ni and (in which ive neglect A. 3 and 

The 

2 4 2 4 

-k 


k k * 

higher powers of A; but not A/ 2 ) is i?i *=-+—, r 


formula cot 2a = 


- which furnishes the exact angle of rotation is 


so simple that there is no point m approximating a, if w e write a = 
j + 0 we have tan 20 = ^ and so the exact angle of rotation may be 


1 


Example 2 Determine the principal axes of the torsion x - 
r cos (0 + Ac) y — r sin (0 + Ac), z = c, k constant Here x * 
a cos kc — b sin kc y — a sin kc -+- b cos kc z — c, and so 


/ COS kc 

— sin kc 


/ o 

0 

I sm kc 

cos kc 

kx ) 

v = 0 

0 

\ 0 

0 

i/ 

\-kb 

ka 


— A&\ 

fca } 

fcV/ 


Since det ij = 0, one of the characteristic numbers of v is 0, and it is 
clear that v(a, b, 0) is an associated characteristic vector Setting 
up the preliminary rotation matrix 




PRINCIPAL AXES OF A SIMPLE TORSION 


3o 


/0 0 0 
2 i - ( 0 0 At 

VO At A-V 


and we have only to deal with the plane 


shear (in the i/c-plane) whose strain matrix is Thus 

if we start with an initial reference frame whose first and second axes 
have the directions of plane polar coordinate lines at the point (a, b) 
in the (a, 6)-planc the principal axes of the torsion are obtained by 
rotating these axes around the radial direction through the angle 

The canonical form of the torsion strain matrix 

At 


t 1 At 
4 + 2 !:m ~2 


is furnished by the formulas m = 0, V 2 


k*r* 

4(4 4- A-V) H ' 1,3 


(4 + *V) W + 4 + 


At 


(4 4- A-V) 




+ 


, O O 

k~r~ 


l; 3 r 3 


4(4 4~ A-r“) 


54- 


If we are 


dealing with a cylinder of length l which is twisted through a quarter 

r r 


turn at c - / we have Id — - and so At - . , 

4 4 l 


Thus the approxima- 


tion furnished by the infinitesimal theory will not be valid at the sur- 
face of the cylinder (for a torsion of this magnitude) unless the ratio of 
the radius of the cylinder to its length is so small that its square may 
be neglected. 

If we use the same reference frame for the initial and final positions 
of the medium wo may introduce the displacement vector v{U, V, IF) 
whose coordinates are defined bv the formulas 


(" «= x — o, 


V 


IF 


It follows that J = /?a 4- K where 7v is the Jacobian matrix 
mnl so M « ,/V ** E- 4 - K* + K 4- K*K. Hence 


(U,V, W) 

(a, b, c ) 


v « &(.V - /v.i) - i(K* 4- IC) + iK*K. 

In the infinitesimal theory of strain, K i* regarded as an infinitesimal 
matrix, i.e,, as a matrix whose elements are so small that their squares 
• and products are negligible. To this degree of approximation, then, 

( U. 4(F n 4 - U h ) m\ + W a )\ 

V, ~ i(/c» 4- K) » i(F« + U',) V b 4(1F, 4- V e ) . 

WA'c 4- IF 0 ) 4(lFfc 4- ID JF r / 

ro obtain the exact expressions for the elements of ;; we must add 
lK*h to h{K % 4- K ) : remembering that the elements of K*K are 
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the squared magnitudes and the scalar products of the column vectors 
of K, we obtain 

<i„ = u. + j(tV + IV + »V), 

<m - V„ + £(£V + IV + IIV), 

*. - ir, + i(EV + IV + iiV), 

me = i(in + v.) + i mu, + v h v, + ir k ir«), 

- iw. + »'.) + HUM, + r,r. + ir.irj, 

<u - i(v. + u t ) + ta >.u t + r.n + ir.ir,) 

3. The invariants of the strain matrix 

We have seen that the strain matrix v is sensitive to a rotation of the 
initial rectangular Cartesian reference frame, under the transformation 
o— » o' « R*a, 17 — ► 17' = R*vR It follows that, no matter what the 
number X is, tj' — \E j = R*( 17 — \Ez)R (prove this) Hence det (n' 
— \E 3 ) = det (ij - XEj) (why?) If we denote by 1 1, J 2 , and lj, 
respectivelj, the trace of v, the trace of co n and the determinant of v, 
w e have 

det (17 — XE3) = Ij — Ii\ + IiX — X 3 

(prove this) If //, TV, I3 are, respectively, Tr ij', Tr co ij', and det ij' 
it follows that, no matter what the number X is, 

/ 3 ' - /A + /A 2 -X s = I 3 - / 2 X + /,X 2 - X 3 
Hence (why?) 

h’ = h U = /*, U = h 

Thus, although tj is sensitive to a rotation of the initial rectangular 
Cartesian reference frame, the three functions 1 1, / 2 , and I3 of the 
elements of 17 (or, as we shall simply say, of »j) are insensitive to any 
such rotation We term them the three invariants of 17 They are 
furnished by the formulas I\ = tj aa + 17&& + Vce, 1° words, 1 1 is the 
sum of the three diagonal elements of 17, 

t P766 We 1 P7<e Veai t r7oa *7ab 
i 2 = + + 1 

F?c6 Vcc\ |17oe T7<xa| *7661 

in words, J 2 is the sura of the three tw o-rowed diagonal minors of r\, 


I3 — det ij, 
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in words, Js is the determinant of y. On writing y in its canonical 
(ts diagonal) form, we see that 

l j «: V ) 4- ijs 4- m, h - mm 4- mm 4- iu»?2, h - mmm- 

If y is mi infinitesimal matrix /,, h, and h are, respectively, infini- 
tesimals of (at least) the first, second, and third orders. 

It is easy to express the compression ratio — m terms of the mvari- 

Pa 

ants 7j, /», and 7a of jj, In fact, p„ = p x (det J) and dot. ( E 3 4 2rj) = 
del M » dot ./*,/ = (det J) 1 and so 

^ = dot (E3 4 2ij) = 1 + 2/j 4 47; + 87a. 

Hence 


— — (1 + 27] + 47; + 87.n) 

P« 

The approximation to this exact formula which is furnished by the 
infinitesimal theory is 

- » (1 + 27 1)“* 4 « 1 - 

Pa 


To this degree of approximation, I \ = — — is the relative decrease 

Pa 

in density or, equivalently (in view of the principle of conservation of 
mass), the relative increase in (local) volume. 

EXERCISES 

J. Show that for the uniform dilatation (or compression) x <=> (1 4 k)a, y ■= 
(1 4 k)b, : "■ (1 4 l)c, k constant, ij «= (k 4 ■jArJEj, and calculate I\, /«, and Iz. 
•lnrtrer, /, » 3 {/; 4 U‘), /. - 3(/: 4 Ik')'-, h ■=■ (k 4 ;t 5 ) s . 

2. Verify that, for the uniform dilatation (or compression) of Exercise 1, /; — 
If i s 1 I j »- uV/i s . 

3. Show that if the relations of Exercise 2 are satisfied the strain matrix is n 
'"■“hir matrix, i.e., a multiple of Ej. /tint. The cubic equation that determines 
the charneteriMir numbers of *is (X -i/j)* -0. 

t. Same that two deformations having the same strain matrix have their 
■Jacobian matrices J t and J * connected by the relation J- «* HJ j where R is a rotn- 
itin matrix. Hint, J\ ~ J. » and so J : ~ EJj where R « 

**• hhow that the compression ratio for the uniform dilatation (or compression) 
r.reret*i* j h furnished by the formula - »{1 4 
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6 Show that for the uniform dilatation (or compression) of Exercise 1 the strain 

1 /p.V M 

matrix is the product of Ei bj —1 — 1 — 1 Hut det M — (1 +2r) J 

where ij = cEj 

7 Calculate the invariants of the shear x «* a + Ac y ** b 2 ■=» c and deter 
mine the relations beta een them 

Answer I\ ■» ■jA' /1 *=■ —■}<:* /j = 0 1 2 *» — yli J* =» 0 

8 Repeat Exercise 7 for the torsion x = r cos (9 + kc) = a cos Ac — 6 sin kc 
y «* r s n (0 + Ac) = a s n kc -f b cos kc 2 ■= c and verify that this is at any 
point {a b c) a shear whose coefF cient is kr = A (a* + 6 1 )* 4 

lnsirer I\ = |t ! r* 1 2 = — ■ jA*r* 1 3 = 0 

4 The compatibility relations between the elements of the 
strain matrix and their derivatives 

Suggestion The mathematical computations involved in this 
section are quite formidable for a reader who is not well trained in 
tensor anal} sis One may safely omit this section in a first reading 
at some later time when the reader has leisure and wants to under 
stand what is behind the compatibility relations of the infinitesimal 
theory of elasticity he may return and work through the computations 
of this section 

The basic relation (dx)* dx - ( da)*M da tells us that the quadratic 
form with coefficient matrix M in the elements of the 3 X 1 matrix da 
appears when written as a quadratic form in the elements of the 
3X1 matrix dx as a quadratic form with the coefficient matrix E s 
whose elements are constants In the terminology of metric differen 
tial geometry this constancy of the elements of the coefficient matrix 
of the quadratic form (dx) * dx is expressed as follows 

The curvature tensor of the three-dimensional metric space the 
coordinates of w hose metrical tensor m the (a b c) system of coordi 
nates are furnished by the elements of M is the zero tensor 

This fact furnishes six relations involving the elements of M and 
their first- and second order dematives with respect to (o b c) or, 
equivalently since 1/ = E 3 -f 2ij six relations involving the elements 
of ij and their first and second-order derivatives with respect to 
(a b c ) These are the (exact) compatibility relations between the 
elements of jj and their first and second derivatives with respect to 
(o 6 c) If ij is an infinitesimal matrix and we neglect in the (exact) 
compatibility relations all products of an element of y or of any 
derivative of an element of ij by an element of rj or any derivative 
of an element of 17 we obtain the compatibility relations of the infini 
tesimal theory of elasticity 

A knowledge of tensor analysis will not be assumed here and the 
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discussion will simply furnish that portion (and no more) of tensor 
analysis which is necessary for the purpose. Following the notation of 
tensor analysis we use (a 1 , o', a“) instead of (a, b, c ) and (x 1 , x*, x 3 ) 
instead of (x, y, z). The double occurrence (once below and once 
above) of a Greek label in a term indicates summation, with respect 
to that label, over the range 1, 2, 3. For example, if x r is a typical 
element of a 3 X 1 matrix x and if y t is a typical element of a 1 X 3 
matrix y, the matrix product yx is denoted by the symbol y a x° (or, 
equivalently, yttx 3 , j/ 7 x’, etc.). In this notation the fundamental 
relation (tlx) * dx - ( da)*M da may be written as 

77i a$a a it a B <t — 1, if s = r, 

== 0, if s r. 

Here is a typical clement (the clement in the pth row and gth 
column) of M and the comma in the symbols a° >r , a 3 ,, denotes differen- 
tiation; a p . r is the derivative of a p with respect to x r . The occurrence 
of the two Greek letters a and /? in the symbol 7 n o3 a n , r o fl ( ., indicates a 
double summation (one summation on a and one on p) so that m„sa" r a fl lt 
is the sum of nine terms (write down five of these terms). We now 
differentiate the. relation = constant with respect to any 

one of the x’s, say x'. In doing this differentiation we remember that 
the elements of M arc given as functions of the o’ s so that the deriva- 
tive of with respect to x' is ?n P7i7 a -r ,i where »j p7 , r denotes the deriva- 
tive of m pt) with respect to n r . We obtain 

w«).-,a a , r a a ,,n y , t + ui aS a a trt a 9 t , + in a) a n %r a B , f t - 0 

where a p tTi , for example, denotes the second derivative of a v with 
respect to x r and x' (so that « p , r( = a p Jr [why?]). We solve this 
system of equations for the second derivatives a p <rt as follows: First 
interchange r and t and rewrite the triple summation m 0 9 ,- ) a a ja B t ,a y tT in 
the equivalent form ?;i 7fltn o n >r n i, ,,o' r . e (why arc these two triple summa- 
tions the same?). We obtain 

”b3.»o <, .rO i! . s a , 1 , + jn o3 o" tr/ o s ,„ + 7u nS a n it a e ,„ — 0. 

Interchanging s and l in our original relation and rewriting similarly 
our triple summation, we obtain 

-F tn„ 3 <i n . r ,o s , ( + m aS a a tT a s t , s = 0. 

On subtracting our original relation from the sum of the two relations 
obtained in this way and denoting by r ps: , the combination F(m rr ,, 



40 THE SPECIFICATION Or STRAIN 

m 9 r p — m p/ j r ) of the first derivatives of the elements of M, we obtain 

r„a r a s , # a T i + m a3 a“ „a B t = 0 

(in denying which we have used the fact that the two double sums 
m a8 a a t a e „ and m a3 a" r ,a B , are the same (why?) (Hint The matnx 
Af is symmetric ) (Note that r = T pq T ) This relation can be 
written in the equivalent form 

(r*a ,o" T o 8 * + m ay a a T ,)a y t — 0 

On multiplying this relation by x l p (the denvativ e of x‘ with respect 
to o p ) and summing with respect to t, we obtain, since a 3 <x T p is zero 
unless j = p, in which event it is unity (whj 7), 

r„a p a a ,o? , + m pa a a r , = 0 

Finally, on multiplying this relation by m 3P and summing with 
respect to p, where m JP is the element in the jth row and pth column of 
M'' 1 , we obtain 

m’T.ii ,a a T a? , + a 3 „ = 0 

We now differentiate the relation T af p a a r a 8 , + m pt a a r , - 0 with 
respect to x‘, interchange s and t in the result, and subtract the two 
equations obtained in this way Since a k = a* t , and a*„» — 
a k T i , v.e obtain 

(r a /J p) 7 a “ r a? ,o 7 t - (r aff p ) r a“ ,a B ,a y , + T a8j> a a rt a B , 

— T a fi p a a r ,a f i + m p f y a B r ,a y t — m P s y a B ri®* . = 0 
We simplify this relation by first writing the triple summation 

(r.j p) y a° rd? t a y , 

in the equivalent form (IV, P ) 8 a“ r a B t a y t (why are these two triple 
summations the same?) and then eliminating the second derivatives 
of the o’s with respect to the x’b by means of the relation a 3 „ = 
— w/Taj ,a a , On the appropnate rearrangement of the summa 
tion labels every term has a “ ,a s t a y t as a factor, and then if we multi- 
ply through by x T q x‘ 3 x l * and sum on r, s, and t there remain only the 
terms for which a — q, p — j and y = k We obtain in this way the 
following relation 

(r,i p) k ~ (Fgt p) j — m rt V 9] pTqi » 4* »F,* p 

— m* f m p€ i 4* vi €, m p , }V 9 k # “ 0 

From the definition of the symbols r P9 , it readilj follows that m pq r = 
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y rft -f r and, on using this to eliminate the first derivatives with 
respect to the <i's of the elements of the matrix M, we obtain 

R CP .jh *= (r,jy ;r ),jt — (r,i ; p).y + tn' f r P ; ;e r Q k;f — m e/ r P i ;c r q};f — o. 

These are the (exact) compatibility relations connecting the elements 
of M and their first and second derivatives with respect to the o’ s. 
The terms involving the second derivatives of the elements of M are 
furnished by (I\,y ; ;,),* - (IV-.a)./ an( * arc ~ m jq.pk ~ m vi-.ii 

-r It follows readily that the curvature tensor R qv -j> ; is alter- 

nating in p and q and that, it is symmetric in the pairs qp and jk. 
'J’herc are, then, only six distinct compatibility relations, namely, 

/?23;23 = 0, R 31;3J — 0, Rl2,l2 ~ 0, 

R" 3;31 0, R% 1-12 = 0. /?12;23 = 0. 

When the strain matrix is treated os an infinitesimal matrix, the 
symbols r., 7;r are infinitesimals of at least the first order and so the 
term'' involving the products of F’s in the compatibility relations are 
of at least the second order. Neglecting second and higher order 
infinitesimals the compatibility relations reduce to 

^pk.jq t m qk,) P ~ 0 

or, equivalently, since M differs from 2 tj by a constant matrix, 

’hr. >7 1 *b ’hfr.jr ~ V)i}.pk ~b npk.jq- 

There are three relations of the type 

’tlJ.tl + T?n,23 - »?J2.3I + ^31. 12 

and three of the type 

V22 . 33 + 1)33,22 ~ 2 , )23,23. 

In our original notation, where we used (a, b, c) instead of (a 1 , a 2 , a 3 ), 

the six compatibility relations (of the infinitesimal theory of elasticity)’ 
are 



-f 

Vmijtf 

~ Vab.ea 

+ 

Vca.nht 


-r 

T] *>, f 

~ Vic, til 

+ 

Vat,, bet 

Vtib.tr 

4* 

Vrr,t>** 

“ Vca,lc 

+ 

Vic, cat 

*lt * ,r<* 

-r- 

Vcr,r 


*» 



-i- 


— O 





'c-t.V- 4 " 
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THE SPECIFICATION Or STRAIN 


EXERCISES 

1. Verify the (infinitesimal theory) compatibility relations by setting 
rjaa ** U a , V be “ + l'e), etc 

2 Show that the (infinitesimal theory) compatibility relations are satisfied by 
every homogeneous strain, l e , by every strain matrix whose elements are constants 
Is this statement valid for the exact compatibility relations? 

Answer Yes 

3 Show that the (infinitesimal theory) compatibility relations are satisfied by 
any strain matrix whose elements are linear functions of (a, b, c ) Is this state 
raent true for the exact compatibility relations? 

Answer No 



3 


THE CONNECTION BETWEEN STRESS 

AND STRAIN 


1. The matrix 

Wo consider any thiec-dimonsional portion V x of our deformable 
medium (when in the final or deformed state) and we denote the sur- 
face of Y x by .S-. Since our deformable medium is supposed to be in 
equilibrium, when in the deformed state, F r is in equilibrium under 
the action of various forces. Those forces are of two types: 

(1) Mass or body forces (such as the weight of the material in 
F r ). p. being the density at a typical point P x :(x, y, z) of V x , the 
tnnss of the element of volume dY x is p x dV x and we denote the force 
per unit mn«s acting on this element of mass by F ; the coordinates of F 
are furnished by the elements of a 3 X 1 matrix 
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the elements of the 3X1 matrix 



and the force acting on V x across the element of area dS z has its coor 
dinates furnished by the elements of the 3 X 1 matrix fdS z where 
dS x = {(dS*)*dS*} M is the scalar element of area of S z We term 
the vector whose coordinates are furnished by the elements of the 
3X1 matrix / the stress on the matrix element of area dS z Thus 
each coordinate of the stress \ector has the dimensions of a force 
divided by an area Stress is force per unit area 

We now make the basic assumption that these two systems of forces 
which maintain that portion of the deformable medium which occupies 
V x in equilibrium would maintain it in equilibrium if it were rtgidified 
This implies that the two systems of forces (mass forces and surface 
forces) which act on the deformable material occupying V x must 
satisfy, when taken together, the conditions imposed upon a system of 
forces which maintains m equilibrium a rigid body These conditions 
are most conveniently expressed as follows 

The virtual uork of all the forces of the system in any virtual rigid dis 
placement is zero 

A virtual rigid displacement 13 a particular case of the more general 
concept of a virtual deformation In order to define the latter we must 
imagine that the coordinates (x, y, z) of P z are functions not only of the 
coordinates (a, b, c) of P a but also of an accessory independent variable, 
or parameter, which we shall denote b> 9 The four variables a, b, c 
and 6 are independent, and if / is any differentiable function of these 
four vanables we write its differential in the form df + Sf where 

df = /„ da + f b db + f e dc, 5f = f,de 

Thus df is calculated under the assumption that 0 is held constant 
whereas 8/ is calculated under the assumption that a, b, and c are held 
constant A virtual deformation is defined when the 3X1 matrix 



is given as a function of x, y, and z, i e , when hx is furnished at each 
point P x of the deformed position of the medium When the 3X1 
matrix hx is such that 5(ds z ) 2 - Owe say that the virtual deformation 
is a t irtual rigid displacement 



VIRTUAL DEFORMATIONS 


4o 


Since the variables a, b, c and 0 are independent the order in which 
second-order derivatives with respect to these variables are calculated 
in immaterial (these second-order derivatives being assumed con- 
tinuous). Thus, for example, x, at — i,« s (where x, a $ denotes the 
derivative with respect to 0 of x„, and x,»„ denotes the derivative with 
respect to a of a>) and, after we multiply through bj- dO, this relation 
appears in the form 5(x„) = (5x) n . If .1/ is any matrix, of which the 
element in the #th row and 71th column is 7n p ? , we denote by 5il f the 
matrix of which the element in the qt\\ row and pth column is 5in p Q . If 
we denote by (6x) 2 the 3X3 matrix 


( (5x) z 

(5x)y 

(ox).’ 

(5x) z = ( (by). 

(5t/)y 


\(o=)r 

(5r) y 

(5-).-, 


(ox, by, 62) 
(x, ?/, 2) 


the relation 5(x„) — (5x)„ and the similar relations obtained by replac- 
ing a by b or c and x by y or z yield 


5.7 = (5x) z ,7 


(if 17 

where ./ is the Jacobian matrix — {Prove this. Hint. (6x)„ = 

(a, b, c) 

(5-r).x n -f ( oT) v y n + (5x) t - n .] It is clear that. 6.7* = ( 8J ')* (prove 
this), and so 


U* = .7‘{(SxLl*. 


Since M = ,/*,/ it follows that 


5.1f « (5./*), 7 + j*u 

- J'UCox),}* + (5x) z |J. 

l\e denote by D the symmetric matrix 


I) i ; |(5.r) r { % -f (5x) z j 

/ (or),. £{(ty) r + (ox) v { 

I v I (*//).- 4 (5x) v \ (or/) v 

\i » < H -?* («-)*■ i A { (52). + (S, ; ) f ( 

and then we have the relation 


$K*r)« + (5r) z |\ 
^•|(52)j, -f- (oy).\ j, 

(5=).- / 


5.1/ - 2J*DJ, 

or, equivalently, since .If ^ E^ + 2y 


J*DJ. 
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Since ( ds x ) 2 = (da)*M da, we have {remember that a is independent 
of 0 so that 5 (da) — 0] 

&(ds x ) 2 - (da)*SMda = 2 (da)*J*DJ da 
» 2(di)*D dr 

For the virtual deformation to be a virtual rigid displacement we must 
have S(ds x ) 2 = 0 for arbitrary choice of the 3X1 matrix dx Hence 
we have the following criterion for a virtual rigid displacement 
A virtual deformation Sx is a virtual rigid displacement if, and only if, 
the symmetric matrix D = £{[(5x)x]* + (Sx) z j is the zero matrix 
It follows that the system of forces, both mass and surface, acting on 
V x must be such that their virtual work, m any virtual deformation for 
which D is the 2 cro matrix, is zero The simplest virtual rigid dis- 
placements are the virtual translations for which the elements of the 
3X1 matrix ix are constant functions of (x, y, z) Setting, m par- 
ticular 5x - 1, by = 0, Sz = 0, we find that the virtual work of the 

mass forces is p x F dV x whereas the virtual work of the surface 

Jv. 

forces is J f dS x Hence 


L>- Fdv ‘ + L 


f dS x 


(which merelj expresses the fact that the resolved parts in the direction 
of the positive x axis of all the forces, both mass and surface, acting on 
V x add up to zero) On letting the volume of V x — > 0 we see that 

lim (V x — * 0) J f dS x *= 0, and on applying this result to a flat cyhn 

dncal volume whose height tends to zero we see that / changes sign 
when the direction of dS 1 is reversed Applying our result to the 
tetrahedron whose faces are the coordinate planes and a plane whose 
coefficient vector is v(u l u 2 , u 3 ), we find that 

/ - fxit 1 + f v u 2 + fzU 3 

where f x is the x component of the stress on an element of area whose 
normal has the direction of the positive x-axis, etc Writing down the 
corresponding equations for g and h, we see that the 3X1 matrix f 
whose elements furnish the coordinates of the stress across thp matrix 
element of area dS x which has the direction of the unit vector v(u\u 2 ,u 3 ) 
is given by the formula 

/ = Tu 
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the other hand, the stress matrix has nothing to do with the initial 
position of the medium, it is determined by the final position of the 
medium Thus it does not make sense to ask Is the stress matrix 
sensitive to a rotation of the initial rectangular Cartesian reference 
frame ? It is logical to inquire, however Is the stress matrix sensitive 
to a rotation of the final rectangular Cartesian reference frame? 


EXERCISES 

3 Show that if the stress on an element of area is always normal to the element 
the stress matrix is scalar le a multiple of 27j Note When T — —pE* 
p > 0 the stress on an element of area is always normal to the element and is a 
pressure rather than a tension This is the situation in hydrostatics and we term 
for this reason a scalar stress a hydrostatic stress A hydrostatic stress is a pressure 
when T is a negative multiple of Ej and a tension when T is a positive multiple 
of Et 

4 Show that m a hydrostatic stress the magnitude of the stress on an element of 
area is independent of the direction of the element of area 

5 Show that the resolved part in the direction of the element of area dS* of 
the stress on dS* is u*Tu where u is the 3X1 matrix whose elements are the coor 
dinates of the unit vector which furnishes the direction of dS 1 


2 The conditions of equilibrium and the virtual work in any 

virtual deformation 

The virtual work of the mass forces acting on any portion V x of the 
deformable medium in any virtual deformation is furnished by the 

i olume integral J p x (Sx)*F dV x and the virtual work of the surface 
forces acting on V x in this virtual deformation is furnished by the 
surface integral J* (Sx)*fdS x Since f = Tu and since (dS z )u = 
dS x the virtual work of the surface forces may be written m the form 
J (5x)*T dS x If, now, £ = (£, y, f) «s any 1X3 matnx the surface 

integral 1 fdS* = (£<*S X + ydS" + t dS') = (£dydz + 

ydzdx + f dx dy) may be w ntten, by virtue of Green s theorem, as 
the volume integral L (fc + 1 J» + f.) dv, = J r (div* {*) 1 17, where 

ne understand bj the divergence (m the i-coordmate system) of any 

3 X m matrix the I X m matrix obtained by taking the divergence of 
each of its column vectors Thus 
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EQUATIONS OF EQUILIBRIUM 


EXERCISE 

1. Show tli when the tire-'* in a hydrostatic procure the equations of equi- 
librium are — (^rad r p)‘ -{- pJ' - Onnd that the equations of motion are - (grad* 
/.)* -- PfT *-> where prnrl.- j< is the 1 X 3 matrix 7b>, ?>.•)• 


On using the fuel t hat T* - 7’, i.e., that the stress matrix is sym- 
metric, in the expression for the virtual work, we obtain (show this) 
the following fundamental result: 

Thi i irlual irorl: of till t he forces acting on any portion V. of our 
thformahh mtdium in any virtual deformation is obtained by integrating 
the trace of the product TD over F-; ^ 


( VI 

Virtual mark = f Tr(TD) dV x . 
Jr, 


C 

L 




We have seen that or, » ./*//./ and so TD = T{J*)- l SnJ~ l . Since 
7'r .•!/>' ~ ’/'r /LI it follows (show this) that 

1 

Tr{TD) » TriJ-'JX.rr'tn). 


Hence mir expression for the. virtual work of all the forces acting on any 
portion l . of the deformable medium in any virtual deformation may 
be written in the equivalent form 
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/ 6 l \ 

and b = ( I are two real 3X1 matrices, the 
V>7 

scalar product of the column vector of b by the column vector of a is 
the 1X1 matrix a*b This scalar product may be written as Tr a*b 
or, equivalently, as Tr b*a If A and B are two real 3X3 matrices, 
either of the two equal numbers Tr A*B, Tr B* A may be termed, by 
analogy, the scalar product of the two 3 X 3 matrices A and B Thus 
the virtual work, in any virtual deformation, of all the forces acting on 
any portion V* of our deformable medium is found by integrating the 
scalar product of J~ l T(J*)~ l and Srj over V x It is clear that this 
scalar product is (as the name implies) insensitive to a rotation of the 
final rectangular Cartesian reference frame In fact, under the rotation 
z~*x' = R*x, R*J, T—> T’ *= R*TR and so J~ x T{J*)~ l 

-» J~ X RR*TRR*(J*)~ X = J~ X T(J*)~ X thus J~ x T{J*r x is insensi 
tive to a rotation of the final rectangular Cartesian reference frame and 
we know that >j, and hence 6ij, is insensitive to any such rotation 
In the infinitesimal theory of elasticity the (approximate) result con 
cerning virtual work may be phrased as follows The virtual nork of 
all the forces acting on any portion V x of the deformable medium m 
any virtual displacement is found by integrating the scalar product of 
T and Srj over V x Since T is sensitive to a rotation of the final rec 
tangular Cartesian reference frame, whereas Stj is not, it is clear that 
the scalar product Tr(TSi)) of T and Stj depends on this final reference 
frame, this makes it evident that the statement of the infinitesimal 
theory of elasticity concerning virtual work cannot be an exact state 
ment, for the virtual work cannot depend on the accidental choice of 
the final rectangular Cartesian reference frame 

3. The virtual work when the stress is hydrostatic 
The expression for the virtual work of all the forces acting on any 
portion V x of the deformable medium, m any virtual deformation, 
takes a particularly simple form when the stress is hydrostatic, l e , 
when T = —pEz In this case J~ X T(J*)~ X = —pJ~ x (J*)~ x = 
— p(iljT) 1 and so the virtual work is found by integrating the scalar 
product of —p{M)~ x and St/ over V x Since M = Ei -f- 2ij &M = 
2St? and so 

Virtual work = — $ J p Tr(M~ x SM) dV x 

If, now, A is any 3X3 matrix, S(det A) is the sum of three determi- 
nants (each obtained by differentiating one of the column vectors of A) 



THE ELASTIC ENERGY OF A DEFORMABLE MEDIUM 53 


jirij 1 

a 2 J 

as) 

a, 1 

Soj 1 

oa 1 ; 

K 

a ~l 

003 1 

15a i* 

o 

«3 : j + 


oa o 2 

< 13 “) + jar 

ay 

Sa 3 : 

|5ai s 

a : 3 

of! 

ot 3 

oa 2 z 

«3 3 ! 

oi 3 

a : 3 

Sa 3 3 


/prove this), and this result may be written in the equivalent form 
6(dct A) = Tr( co A)*6A. 

If A is non-singular we obtain, on dividing through by det A, 
Tr(A~ l oA) = (det A) _1 5(dct A) 

and ko the virtual work may, when the stress is hydrostatic, be written 
in the form 

Virtual work = — f p(det Af)“ 1 5(dct M) dV x . 

~ J v. 

Since dot M » (det J)~ we have 

Virtual work = — f 7 r(det ./) -1 <5(del .7) dV x , 

(IV 1 

and since det J = -r~> so that 5(det«7) = - o(dV x ), we obtain, 

ft d\ a 


finally, 


Virtual work 


Sr. 




If p is a constant function of (r, y, z), i.c., if the stress is a constant 
hydrostatic, pressure (or tension), this equation simplifies further. On 

writing V~ ~ f dV x « f r/F a , we obtain 5F, - j —~~~dV a 

Jy. Jr.dV„ Jy. dV a 

n j (Note that, in the integration over V„, the limits of 

J y, 

integration are independent of the accessory variable, or parameter, 
0. whereas in the integration over V x the limits of integration vary, in 
general, with 0.) Hence, when the stress is a constant hydrostatic 
pres.-ure, or tension, we have 

Virtual work = — p j o(dV x ) — — poF~. 

J y. 


The elastic energy and the relation between stress and strnin 

We ft*-'Uinc that the work done by all the forces acting on any portion 
V f of our deformable medium in any deformation is stored up in V x ns 
rauftV energy or energy of deformation (none of it being dissipated in 
hwtO The energy of deformation may well depend not only on the 
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actual state of strain of the medium but also on the previous history 
of the medium, in other words, the energj of deformation may be 
distributed throughout V x with a mass density ^ 

Energy of deformation = L Pi’/' dl r i 

where ^ is a function not only of n but also of the previous values of 
17 (from 0 to 17) Such hysteresis effects are matters of common expen 
encc but we shall ignore them here We make the fundamental 
assumption that the mass density ^ of the energy of deformation (1 e , 
the energj of deformation per unit mass) is a function of the strain 
matrix rj so that $ is unambiguously determinate when tj is given 

\fr ~ tf-(lj) 

Remark 1 It is more convenient to deal with the cnergy-per-umt 
mass rather than the energj -per unit-volume since the mass of anj 
clement of volume dV x of the medium remains constant in anj defor- 
mation whereas dV s changes, in general 

Remark 8 Since 17 is a sjmmetric matrix, ^(jj) is unambiguouslj 
determinate when the six following elements of 17 are given — ri aa , 
Vbb, Veet Vbei Vca, Vab — and we can write 

= vK’faa, Vbbr *7ee> Vbct Vcoi lja&) 

However, the gain obtained in reducing the number of variables of 
which ^ is a function from nine to six is overbalanced by the loss of 
sjmmetrj involved in this reduction We shall agree to write any 
function of the symmetric matrix 17 symmi tncally where we under- 
stand that f(ij) is written symmetrically if /(A) = /(A*) where A is 
any 3 X 3 mntnx If /fa) is any function of 17 it may be written sym- 
metrically by setting /fa) = / = pfa), for it is obvious that 

g(A) = g(A*) where A is any 3 X 3 matrix (prove this) The mam 
advantage of writing a function /(17) of 17 symmetrically is that we do 
not lo^e symmetry when we take the gradient of / with respect to 17 
Just as the gradient of any function of (x, y, z) is the J X 3 matrix 
obtained bv taking the three partial derivatives of the function with 
respect to x, y and z, so the gradient of any function f{A) of a 3 X 3 
matrix i isthe 3 X 3 matrix obtained by taking the partial derivatives 
of /(A) with respect to the various elements of A Thus the clement 
in the pth row nnd 9th column of the gradient of /(A) with respect to A 
13 the partial derivative of /(A) with respect to In calculating 
these partial denv fltiv cs we pay no attention to the particular structure 
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of A ; we regard the nine elements of A as independent whether they 

actually are independent or not. Thus, in particular, if -4* = A we 

regard a/" as constant when we are differentiating with respect to 

a « even though we know that a q p - a p v . We denote the gradient of 

Of , df /df\* . . 

/(-I) by the symbol —• It is clear that — = f — J and so we have 

the following important result: 

If J( A) is a function , trriV/en symmetrically, of the symmetric matrix 
Of 

A then — is, like A, a symmetric matrix. 

Example 1. The first invariant /j = y aa -F yib + Vcc is a function, 

d J 

written svrnmetricallv, of the symmetric matrix y. — is the sym- 
' “ dy 

metric matrix Zfj: 


3«t 


E 3 . 


Example 2. The second invariant 7; - (ijMtj ce — y irVcb) + 

— y fn y nr ) *k — VoiVbn) is a function, written symmetrically, 

< 7 / « 

of the symmetric matrix jj: — is the symmetric matrix I iE 3 — y' 

f)y 


dy 


hEz 


Example 3. The third invariant h — del y is a function, (sup- 

Dl 

puddly written symmetrically), of the symmetric matrix y. — - is the 

dy 

symmetric matrix eo y: 

dh 
dy 


cn y. 


I he principle of conservation of energy may be expressed as follows: 
Die virtual work of all the forces acting on any volume F, of 


our 
PaMF*. 


deformable medium, in any virtual deformation. = 6 j 

r r J v, 

f)n writing p.f dV. in the equivalent form Pa&dV„, we hav 

’ J^. <v ”’ *■ 5 -* r ’ n.oV'dl'fl *>= | p.of <!V t , and so the virtual work 
U.e, the integral over F. of Tr{J- % T{J *) ~'cy,) is the same as the 
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integral over V x of p z fy Since this must hold for an arbitrary volume 
V x of our deformable medium, we have (why?) 

TriJ-'TV*)- 1 ^) = p z 5 f 


Since yp is (by hypothesis) a function (written symmetrically) of ij, we 


have (by the very definition of a differential) = Tr 
this), and so 




(show 


7Vy-'3V»)-‘Jl) - P.Tr J,) 

Since this relation must hold for an arbitrary (symmetric) matrix 
we have (why ? ) 


or, equivalently, 


j-'tw)-' 




On setting <t> - p a p, so that <f> is the energy of deformation per unit 
initial volume (why?), we obtain, finally, 



This is the fundamental relation of elasticity theory connecting stress 
and strain Ij 4> is known as a function of 17 and if J is given, we can 
calculate the stress matrix T from this relation (17 being determined by 
means of the formula 17 = — E 3 ) and the compression ratio 

being determined by means of the formula — = In the 

p a det J 

infinitesimal theory of elasticity the difference between J and E s is 
neglected, (so that the difference between the compression ratio — and 

Pa 

d<f> 

1 is neglected) and so T = — In w’ords, 

dr] 

To the degree of approximation afforded by the infinitesimal theory 
of elasticity, stress is the gradient of energy density with respect to 
strain, that is, the stress matrix is the gradient, with respect to the 
strain matnx, of the energy of deformation per unit initial volume 
Be sure to understand that this is not an exact statement, it is only 
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an approximation to the exact result, which is furnished by the formula 

\Pa/ (>V 

Since J*J M we hare J - where R is a rotation matrix; 
hence the fonnula furnishing T may be written as follows: 


T 


~ RM k — M h R\ 

Pa dg 


On comparing this with the relation 

JJ* = RMR*, 

we obtain the following result: 

In Dial particular final rectangular Cartesian reference frame in which 
thr coordinates of JJ* arc furnished by the elements of M the coordinates 


of T are furnished by the elements of 


M 


» ,.i S 

dr; 


if, in particular, the initial rectangular Cartesian reference frame is 
so ohm-on that M is in diagonal form, we see that: 

In the final rectangular Cartesian reference frame that is furnished by 
the principal axes of JJ* (he coordinates of T arc furnished by the ele- 
ments of 




3 6 

”'i r~ 


(mpihj)' 

: d(() 

dr] ah 


(ntsmi) M 

36 

dtp,* 

3<f* 

1)1 2 

‘ drj(,b 


(ntaWa) 54 

(»l3W,) S> 

<!<*> 

(nhm 2 )- 

t d<t> 

36 




dr]r n 

driti 

dl]cc 
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of unit magnitude and directed an ay from the medium) Since 
dS z u = dS s we have 

T dS * « dSJ 

Since T — J~J* and 

dS* = (del J) W) 'dS‘ = ^ (./*) <fS“ 

(See Chapter 1 Section 3) this relation may be written m the form 
dd> 

J -2- dS a - dS z f 
dr/ 

W e denote by T a the 3X3 matrix J — and obsen, e that T u dS a — 
drj 

T dS x , in words, the result of operating on the final matrix element of 
area dS x uith the stress matrix T is the same as the result of operating on 
the initial matrix clement of area dS a mth the matrix T a Since, in a 
deformation problem, it is the initial matrix element of area dS a which 
is given (the determination of the final matrix element of area dS x 
being part of the problem ) it is mth the matn\ T a and not with the 
stress matrix T, that we must operate when we wish to satisfy the 
boundary conditions These conditions are 

T a dS a = dSJ T a ~ J — 

dri 

If, for example, part of the surface of our deformable medium is free 
from applied force / is the zero matrix over this part of S x and so 

T a dS a — 0 Since J is non singular it follows that — dS a — 0 and 

dr} 

dif> 

so dS a is a characteristic vector of — corresponding to the character 
drj 

istic number zero Note In some problems the stress matrix T 
(rather than the stress f) is furnished over S x Thus the medium may 
be subjected to a constant hydrostatic pressure Then the boundary 
conditions simply state that T must coincide, over S X) with the given 
stress matrix For example, an important problem that we shall later 
examine thoroughly is the following Determine the deformation of a 
circular cylindrical tube under the action of internal and external 
pressures p, and p e Assuming the displacement of each particle of 
the cylinder to be radial and a function of the distance of the particle 
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from the axis of {ho cylinder, the deformed tube is again a circular 
cylindrical tube. The vector T„dS a must have the direction of the 
radius and, if r. nnd r e are, respectively, the initial internal and externnl 
radii of the tube, the magnitude of Tads' 1 must be p,-dS x when r = 

n and p, dS. when r « r r . From the relation dS~ — ^ dS a , 


we obtain 


dS T 




1 (dS a ) *./ “ 1 ( J *) - 1 dS a } }6 
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1. The energy of deformation of an isotropic elastic medium 

We have seen that the strain matrix 17 is sensitive to a rotation of the 
initial rectangular Cartesian reference frame Under the rotation 
a — > a' •» R*a, ij — * vf = R*rjR This change of the strain matnx does 
not mean that the rotation of the axes affects the strain or deformation, 
the strain is the same as before but its coordinates (1 e , the elements of 
the strain matnx) with respect to the new axes aTe, in general, different 
from its coordinates with respect to the old axes The energy of 
deformation, per unit initial volume, must be insensitive to the rota- 
tion a a' = R*a of the initial rectangular Cartesian reference 
frame This fact implies, since ij is sensitive, in general, to such a 
rotation, that the form of the function which furnishes the energy 
per unit initial volume, must be, in general, sensitive to a rotation of 
the initial rectangular Cartesian reference frame In the new refer- 
ence frame the energy of deformation, per unit initial volume, will be 
furnished by a new function <J>' of the new strain matnx ■>/« and the 
fact that the energy of deformation, per unit initial volume, is a scalar 
(1 e , that it is insensitive to any rotation a — * o' = R*a of the initial 
rectangular Cartesian reference frame) assures us that 

4>(v) “ 4>'W) = <f>'(R*i}R) 

If it should happen that the form of the function 4> is insensitive to a 
given rotation a — > a' — R*a of the initial rectangular Cartesian refer- 
ence frame, we say that the medium is elastically insensitive to the 
rotation R Thus the medium is elastically insensitive to a given 
rotation R if, and only if, <j>'(R*tiR ) = <p(R*ijR) for an arbitrary sym- 
metric matnx rj In view of the relation 4>(v) = we can 

rephrase this definition of elastic symmetry (with respect to a given 
rotation of the initial rectangular Cartesian reference frame) as follows 
The deformable medium is elastically insensitive to the rotation a — » 
60 
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a' ~ It* a of the initial rectangular Cartesian reference frame if, and only 

if, 

£ 0 ?) = £(#*>?/?), 

thin relation being an identity in the ( symmetric ) matrix y. 

Instead of regarding t lie elements of y' - R*yR as furnishing the 
coordinates of the strain with respect to (he new axes we may consider 
the new strain who«e coordinates with respect to the original axes are 
the same as the coordinates of the old strain with respect to the new 
axes (so that the coordinates of the new strain with respect to the 
original axes are furnished by the elements of jj' = R*yR). Then the 
mlntion <f>(>7) “ $(R*yR) may Ik: interpreted ns indicating that the 
energies of deformation, per unit volume, of these two (in general, 
different) strains arc the same. The new strain is that which would 
he produced if the original strain were preceded by the rigid rotation 
a a ' - R*a of the medium. Note that this is an actual rotation of 
the medium and not merely a rotation of the initial rectangular 
Cartesian reference frame. 

A deformable medium is said to be elastically isotropic (or simply 
isotropic) if it is elastically insensitive to every rotation of the initial 
rectangular Cartesian reference frame. Thus, 

.■1 deformable medium is isotropic if, and only if, 
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30 d<^ 30 

isotropic medium r— — — Ej + — (I \Es — 7 ) + rr* co ij, hence 

a ij oi 1 di 2 d/ 3 

when the initial rectangular Cartesian reference frame is so chosen that 

« (and hence \I and l/ 51 ) are diagonal so also is It follows that 

dtj 

— commutes with AT M so that 1/ M ~ 1/ Ji = Hence, 

3>j dij 9 17 

For (in tsotroptc medium the coordinates of the stress tensor T arc 
furnished in the final rectangular Cartesian reference frame in uhtch tie 
coordinates of JJ* are furnished by the elements of M = J*J, by the 
( Px\ 30 

elements of the matrix 1 — 13/ — 


In the final rectangular Cartesian reference frame that is fur 
nished by the principal axes of JJ* T is diagonal (it being understood 
that the axes of the initial rectangular Cartesian reference frame are 
the principal axes of AT = J*J) 



30 30 . 30 30 , . . 

where = ~r 4- (172 + nz) ~rr + v*vz irr etc Thus the principal 
3iji dli 3/* 3 / 3 

axes of T are the principal axes of JJ* In an arbitrary final rectangu 
lar Cartesian reference frame wc hive 



where J => Note The principal axes of ij arc the principal 

axes of J*J whereas the pnncipal axes of T are for an isotropic 
medium the principal axes oi JJ* Thus the principal axes of T are 
not even for an isotropic medium in general the same ns the pnncipal 
ixes of v If howexer wo follow the deformation b> the rigid rota 
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/i 3 — 2m III 2 + n/ 3 We shall not carry the development of 

<(>, as a power senes in 1 1 , J 2 , Is, farther (so that we agree to neglect, 
in the development of <f>, infinitesimals of higher order than the third 
in the elements of the strain matnx tj) We have, then, 

<f> = «£0 + + $2 + <f>3, 


4>i = all, 


d<t> l 
di) 


= Bs, 


fa - ^ Ii* ~ ^ “ \IiE 3 -{- 2fiij; 

& OIJ 

*= - *-±p! V - 2m/,/, + n/„ a -ii 

O OTJ 

= (II i l — 2 mI 2 )E 3 + 2mlir) -f n co 17 , 


— = aE 3 + (\IiEj -f- + (II i 2 — 2m/ 2)^3 d - 2m/ iij d* n co 17 

*> 

Hence T - /1<«E, + (X/,E, + 2m) + (Hi 1 - 2 mI,)E, + 

2 ml\i) + « co When the strain is zero, 17 is the zero matrix 

and J is a rotation matrix On denoting by T 0 the initial stress 

matrix, 1 c , the stress matrix when 17 = 0, it follows (since — = 1 

Pa 

w ben 7 \ = 0 [w hy?J) that 

To =* aE 3 


In other words, the initial stress is either a hydrostatic pressure or 0 
hydrostatic tension The only assumption about the deformable 
medium which we have made is that the medium is isotropic and this 
has been enough to force the conclusion that the initial stress must be a 
scalar stress, i e , a hj drostatic pressure or tension Hence, 

No elastic medium that is initially in a state of stress which is not 
scalar, te , of the nature of a hydrostatic pressure or tension, can be 
isotropic 

Assuming that the medium is isotropic, we denote the initial hj dro- 
static pressure by p 0 , then a = —p 0 and the stress matrix T, when tlie 
strain matnx is 17, is furnished by the formula 

r = (j)/{ — poE, + (X/,E, + 2m) + (//,= - 2 m/,)E, 

d- 2m/i»7 + n co rt\J* 
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'I’liia formula involves the six constant* po, X, /i, I, in, n (the significance 
of t lie t»>rm constant being that they are independent of the strain 
matrix r?). Of these six constants the first is the initial pressure and 
so its dependence on the initial pressure is obvious. It is natural to 
expect, then, that the remaining five, namely, X, g, l, m, and n, depend 
also on p 6 . These five constants are known as the clastic constants of 
the isotropic medium. In the infinitesimal theory, which neglects, in 
the formula for the stress matrix, infinitesimals of the second order in 
the elements of 17 , the three elastic constants (l, in, n) do not appear, 
and T is furnished by the formula 

T - (-) ./{ -Vo Ei 4 (X/iEn 4 2p» 7 )}./* 

V>«/ 


where the compression ratio — — 1 — / j . The two elastic constants 

Pa 

X and n which appear in this formula are known as the clastic constants 
of l Mint (after Gabriel LamO [1795-1870], a French mathematician). 

Warning. Do not be misled by the term constant into overlooking 
the fact that X and n depend on the initial hydrostatic pressure. 

Writing J - RM”, we have, to the degree of approximation con- 
templated by the infinitesimal theory, 

irril - (1 - /,)(£, 4 217 ) | -roEl + (X /,£ 3 4 2 / 117)1 

- —polis + (X 4- Po)I]Ez 4 2(/i — po)n- 

When the initial stress is zero, this equation reduces to 

11 * TR — Xq/ j/y 3 4 2fiov 

"here (X 0 , go) are the values of (X, /j) that correspond to po = 0, i.c., to 
zero initial pressure. On comparing this with the previous relation, 
which may be written in the form 

R*(7 4 pt>Es)R ~ (X 4 Po)/iE 3 4 2(a — po)n, 

it is clear that the numbers X n and no that appear in the formula fur- 
nishing the stress, when the medium is free from stress in the initial 
state from which the strain is measured, are the values of X 4 p 0 
and n ~ po. respectively, when p 0 is zero. If, for example, we neglect 
the dependence of X and a on p 3 the effect of a non-zero initial hydro- 
Malic pressure is in so far as the increment T 4 pc,Ez in the stress 
matrix E concerned, to increase X by p 0 and to decrease a by p 0 . 

In the ^second-order approximation, in which we neglect, in the 
formula for 1 , third-order (but keep second-order) infinitesimals in 
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the elements of the strain matrix i;, the formula for R*TR is 


R*TR = - (E* 4- 2ij){ -p 0 E, + (X7 t E s 4- 2m) 

Pa 

4* (E i 2 — 2mI 2 )Es 4* 2m/\tj 4" ti co tj| 

where now ^ - (1 + 27, + 41 ,)-" = 1 - 7i + (| 7, J - 27.) On 

collecting terms, we obtain 72*772 ** — p 0 Ej + (X + po)/iEj 4* 
2 (p — po)»l + 2(X + po — p + m)Ii + [(l — X — %po)h* — 2 (m — 
Po)/ 2 1E s + 4/iij s 4- n co v When the medium is initially free from 
stress this reduces to R*TR = (Xo7iE3 4* 2/i 0 ij) + 2(X 0 — + m 0 )7 1 

4* { (io — Xo)/i" ~ 2 mo/ 2 jE 3 4* + no co ij where (fo m<> no) are 

the values of (l, m, n) that correspond to po = 0 

In comparisons of theory with experiment the compression ratio 

— is often furnished by experiment and it is unnecessary to use the 

Pa 

approximation (1 — 7i in the linear thcorj and 1 — h + (f/i s — 
2 / 2 ) in the second order approximation) If we do not use these 
approximations we obtain the following results 

Linear Theory 

R'TR = (E, + 2,)i-p„E, + (X7,E a + 2m) | 

- (^) (“PoEj + X7.C, + 2b - p,),j 


When po = 0 this reduces to 


R*TR 


(X0/1E3 + 2vov) 


Second order Approximation 
R’TIl = (£) (£, + 2,)| -v,E, + (X7.E, + 2m) 

+ (l/i* — 2mIi)Ez + 2m7i»j 4* n co »j) 

— HPoE 3 + \I 1 E 3 4- 2(^ — p 0 )v 
4- (lh 2 - 2mU)Es 4- 2 (m 4- \)hv 4* « co ? + W\ 
When po = 0 this reduces to 

R*TR *=* lXo7iEa 4* 2/ioij 4~ (fo7x s — 2mnIz)Ez 

4" 2 (mo 4" Xo)7ii; 4- «o co ij 4" 4/ioij 2 } 
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3. The stress in nn isotropic medium when in a state of scalar 
slrnin 

We shall write our scalar strain in the form 


7) = — CiFs 


(the reason for the negative sign being that most of the experiments 
with winch wc wish to check the theory deal with compression rather 
than dilatation. If we wrote jj — cE 3 the coefficient c would be 
positive in dilatation and negative in compression; when wc write 
tj *s —c/43 the coefficient c is positive in compression and negative in 

dilatation'). Since the medium is isotropic, — is a scalar matrix (i.c., 

drj 

n multiple of E j) when n is a scalar matrix (prove this) and it follows 
(why?) that R*TR is a scalar matrix. Hence (why?) T is a scalar 
matrix (being the same as R*TR), and on setting 


T 


•?E s 


(so that p is the pressure in the medium) wc obtain the following rela- 
tions between r and p in the linear theory and the second-order approxi- 
mation, respectively: 


I. i m nr Theory. 
written ns follows 


Since (1 — 2c) 



the formula for p may be 


p » 



Ipo + (3X 4- 2p)cj. 


When wo um» the approximation 1 4 3c to — this equation becomes 

Pa 
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in other words, the compressibility is the derivative with respect to p 
of the logarithm of — Terming the reciprocal of the compressibility 


the incompressibility, we have the following result 
To the degree of approximation furnished by the infinitesimal (or 
linear ) theory of elasticity the (local) incompressibility t $ given by the 
formula 

Local incompressibility => X + -f p + £po 


When the initial pressure is zero the local incompressibility is \q + $po 
Note carefully that the (local) incompressibility 13 sensitive to the 
value of the initial pressure If we neglect the dependence of X and u 
on po (so that we replace X and u by Xo and ua, respectively), the (local) 
incompressibility is the linear function Xo + §po + -J-po of po If w r e 
make the less drastic assumption that X and u are linear functions of po, 
the (local) incompressibility is a linear function X 0 + •§-/■< o + kpo of 
po (the coefficient k of p 0 being a constant that must be determined by 
experiment) Assuming that the strain and po are constant over the 
homogeneous medium, we have, then, 

(dp ° g v)p~p, Xo + ■f/to + kp 0 
Since po is arbitrary, this equation is equivalent to the relation 

d_. 1 = 1 

dp ° g V X 0 + f-Mo + kp 

On integrating this relation, we obtain the following equation of slate 
for the isotropic deformable medium 

i a. L - (V*\* 

+ M+im” \r) 

or, equivalently, 

p = ((£)*_,} 


where Fo is the value of V which corresponds to p = 0 The value 
of k which corresponds to the (drastic) assumption that X and u are 
independent of po is ■£, and the correspondmg equation of state is 

V - (3X + 2„) {(£)" - 1 j 
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The Second -ord rr Approximation. Here 

p = (^j" |p 0 + (3X + 2n)c - (9/ -f n)c 2 }, 
and (bis reduces, when po — 0, to 

p — {(3Xo + 2p 0 )c — (9?o + »o)e 2 |- 

Since (-^) - (1 — 2c) -54 = 1 + c + ^ c 2 + • • • , these results 
VW 2 

run}' be written in the equivalent forms 

p ~ po + (3X -{- 2/r -f- po)c + (3X -f- 2p — 9f — « -f f-po)c 2 , 

p — (3Xo + 2po)c + (3Xo + 2po — 9?o no)c". 

Limiting ourselves to the case in which p 0 and c are constant over the 

medium, we have — — ~ = 1 + and so (1 — 2c) 3i = 1 -f —— 
pz Vo Vo E 0 

2 


TI 1 AV . 

Hence c — — - — — V 


•ii. ( AT V 2AV 1 f A 'V 

■ 1 - 2r = (' + n) " 1 f in - 5 {rj + 

— ) + * • • and so the formulas 
o/ 

ows: 

A P ~ 4 r ~ p + ~ paj 


for p may ho written as follows: 

i \ AV 

Po 


+5( x+ l , -“-i* + W(^) 

4 '-Y. 


v - - (x. + 1 »,) ~ + i (x. + 1,„ - a, - 1 (0 

ifiop y E denoted, for the moment, by £, it follows that 


dp 

A 


l> A , 2 1 \ d J p 2 ,2 4 

s - K + 5" + 3 t ) rfF ■ x + 3“ - 2( - 5" + 3 *• 


(hi ditierentiatiiiE the first of these relation 1 - with respect to £ and sub- 
Mitutinc the result in the second, we obtain 
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and if this value is dn ided bj the expression for ^ it follows that 


U'+b)’ 


f(X ~h ■§ it) — 21 — yn -f* 

* + f p. + ip 


For values of p which are negligible in comparison with Xo + ipo a 
good approximation to the fraction on the right is 

2 [ 9I 0 + n 0 l 


3Xo + 2/ioj 


Thus the k that appears in the equation of state 
(Xo + -§mo) (/^o\* ,1 


'ffl 


of the linear theory is (approximately) 1 — 


2(9fp + rip ) 


The fact 


3(3X„ + 2p 0 ) 

that Xo + fao maj be determined by experiment (as the value of 

37 when p = 0) furnishes a method of determining by experiment 
d? 

the combination 9l 0 + n 0 f 

Remark 1 We have termed the equation 

(Xo + |mo) l/Vo\ k A 


w 


the equation of state of the linear theory Actually, the equation of 
state of the linear theory is the equation 

p = (3Xo + 2 ^ ci) e, 

which we obtain by setting p 0 — 0 in the equation 


p = 


jpo + (3X + 2 fi)e) 


If we use the approximation — = 1 -f- 3e this linear approximation 

Pa 

appears m the form 

p — (3Xp + 2fto)e ~ ^X 0 + - poj — — ^X{ 
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This equation is in disagreement with experimental results unless — 

is infinitc'imal. and we shall not slop to discuss it. The equation that 
we prefer to call the equation of state of the linear theory is really the 
equation of state of an integrated linear theory. We use the linear 

AE 

theory to express Ap as a multiple of — not merely at p = 0 but also 

at nn arbitrary value of p; the multiple X -J- 4 h> Is nn unknown 

function of p, and we approximate it by a linear function of p (the 
coefficient of p in this linear function being determined by experiment). 
The final relation between p and V’ is obtained by integrating the rela- 
tion between dp and dV obtained in this way. 

llrmark £'. The formula p — 1;) 0 4- (3X 4- 2 n)c — (9/ 4 

n)e 2 } of the second-order approximation may be checked against the 
5 ( ! f r 

relation p -- — -- where ‘l* = V n tf> is the energy of deformation of the 

V * T 

entire medium. In fact, on setting 7i - —3c, / 2 = 3c 2 , h — — c 3 in 
the formula 

X 4- 2m / — 

$ « -pel t 4 1\~ - 2pl» 4- / 1 3 - 2 tnlih + nl 3 , 


we obtain 
and so 


d> ~ Spue 4* * (3X 4 2/i)c* — (9/ 4- n)c" 

4 - V„ 1 3;v 4- ?(3X 4- 2a)r 2 - (9/ 4 n)c z \. 


/I’A 2 

_i ji ~(AX’\ 

I ~r ) we have e 

\1 J 

- 1 \' J I 

O'l 


Hence - — - - 
x 

(-) i* + < 


i /J'A“ i w* 

- S \r.) k , ;™” s - Si-, ” {rj <p° + < 3X + 2 ") c - 

(M + n)r 4 | <-(^\ ij>« + (:» + ‘.W - (!>( + Ti)r-|. 

t. Cotnpnri*on or theory with experiment in the cn^e of extreme 
hxdro-tnttc pressure 

The experiments of Bridgman on the compressibilities of various media 
up to the extreme procure of 10 s atnm-pheres f = 10 s grams per square 
e* ntimeter) are available to u-t the various relations connecting 
p-v-ome with volume which art* fumidied bv the various theories we 
h'.ve considered {the simple linear theoty, the integrated linear 
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theory and the second order approximation) Bridgman 1 gives the 
V 

relative volume — at pressures up to 10 5 atmospheres (1 atmosphere 


being I kilogram per square centimeter) The values for the metal 
sodium are as follows (the unit of pressure being taken as 10 4 atmos 
pheres = 10 7 grams per square centimeter) 


V " 


25 

3 

4 

5 

6 

7 

8 

9 

10 

To ' 

X To/ 

-0789 

0770 

0 737 

0 708 

0 683 

0 661 

0 641 

0 623 

0 606 

AV 
~ To 

V 

" 1 ~ To 

-0211 

0230 

0263 

0 292 

0 317 

0 339 

0 359 

0 377 

0394 


It is immediately clear that the formula furnished by the (simple) 
linear theory is in disagreement with the experimental results This 
formula is 


V 


-(*+!-)£ 


and if we use the experimentally determined value of — at p = 2 5 

To 

to determine the value of the (single) constant Xo + Imo appearing in 
the formula we find that 


Xo + !po - 11 85 

The following table gives the comparison between the experimentally 

V V 

determined values of and the values of — calculated from the 
To To 

formula 


p - 


V_ 

To 


i + ^ = 1 r_. 

To Xo + § Mo 


25 3 4 5 6 


- 0 0844p 
7 8 9 


10 


Proc Am Acad Arts Set Vol 76 pp 55-70 (1948) 
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When p *=> 4 


the value of 


JP 

Fo 


given by the formula is too small by 


about. 10 per cent. The experimental values of — have an accuracy 

l 0 

of 2 or 3 per cent, and so the simple linear theory is not valid up to a 
pressure of *1.10' atmospheres (less than twice the pressure used in 
determining the constant Xo -f- -ilaci of the formula). Y\ e shall not con- 
sider further, therefore, the formula furnished by the simple linear 
theory and shall pass to a study of the formula provided by the inte- 
grated linear theory. 

There arc two constants, namely, k and Xo + %no in the formula 


_ (Xo + jpy _ 5 1 


of the integrated linear theory. If pi and p» are two pressures to which 
correspond the volumes Fi and P«, respectively, we have 



and so p z (jrj ~ P> ~ 0>2 - P\) = 0. The value of k 

may lie determined from this ccpmtion (by the method of trial and 

y y 

error) by using the. experimental values of — - and ~ : k having been 

r i V 2 

determined, the. value of Xo + ®/io may he obtained from the equation 


a . 2 kp« 

Xo + g po - 77777" 


In order to see how the integrated linear theory works both backwards 
and forwards, i.e., when we extrapolate to lower as well as to higher 
values of p, we select the values pi = 5, p ; - G in order to determine 
the two constant*., k and X f; -f gpe. of the formula. Our calculation 
will s»rve to determine the value of the integrated linear theorv as a 
mean* of pmHdion; thus we may imagine that Bridgman had suc- 

c*>Hlr-d in measuring r only up to p = G (i.e., up to a pressure of 
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60,000 atmosphcroa) and we wish to predict the values of up to 

Vo 

p =» 10 (i e , up to a pressure of 10 5 atmospheres) If our predicted 
values check well with the experimentally determined values we may 
use these \alues to predict the (as yet experimentally undetermined) 
V 

values of — up to p = 20 (l e , up to a pressure of 200,000 atmospheres) 

' o 

The equatioh that serves to determine K is 



Vi Vt 

where — =» 0 708 — = 0 683 On denoting the left-hand side of 
► o V 0 

this equation by /(A.) we find that (check these calculations) /( 3) ■= 
0 214, /( 4) - -0 0965 /(3 8) « -0 006, /(3 79) - 0 001 In view 
of the errors of experiment the value 3 79 for k is close enough Using 
this value of k we, have 



22 74 
3 242 ' 


7 014 


The formula (of the integrated linear theory) that furnishes — is, 

' o 

accordingly, 




+ 1 851) 


V 

The following table gives the comparison between the values of y~ 
computed from this formula and the experimentally measured value 



V = 2534 

5 6 

7 8 9 10 

V i 

— (measured) - 0 789 0 770 0 737 
V a 1 

0 703 0 683! 

0 661 0 641 0 623 0 606 

V 

— (calculated) =» 0 798 0 776 0 73S 

0 708 0 683 

0 662 0 643 0 627 0 613 


The greatest divergence between the calculated and the measured 
values is at p = 2 5, where the calculated value is in excess of the 
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measured value by less than 1.5 per cent. Since the experimental 
value.' do not claim an accuracy of more than 2 or 3 per cent we may 
<j ay that the values furnished by the two-constant formula of the 
integrated linear theory check the experimentally measured values 
over the entire range of pressures from p — 25000 atmospheres to 
p -■ 100,000 atmospheres. 

Tliis check between theory and experiment gives us confidence to 
V 

predict the values of — - that will be measured when the experimental 
>o 

V 

technique is sufiicicnfly developed that measurements of — up to 

p s- 20 (i.e., up to a pressure of 200,000 atmospheres) can he made. 

We me the experimentally measured values of ~ at p = 9 and at 

V 0 

7) «- 10 to redetermine the values of 1: and ,\ 0 + fiMo- The equation 
that serves to determine 1: is now 


10 



= 0 


T, F„ 

where ~ .023, — - .000. Wc find that (perform the calculation 
1 0 > 0 

or verify t ho result) I: - 2.80, Xo -F fjao - S.9GS. Note. The diiTer- 
ence between these values for /„• and X 0 + aMo and those obtained 
previously (when was 5 and p 2 was 0) merely indicates the differ- 
ence between the linear approximation to X -j- in the interval 
0 $ p 10 and the linear approximation to X + 571 in the interval 

5 £ p ^0. The formula determining “—is, now, 

> 0 

- (0.31 9)(;> + 3.135) 



and the comparison between the experimentally determined values and 
the computed values over the range p - 2.5. 3, • • * , 10, together 
'ahh the predicted values over the range 7) ~ 11, • • • , 20, is furnished 
hv the table on p. 70. 


Tim- the theory predicts that a pressure of 200,000 atmospheres will 
compro" sodium to 0.-197 of its original volume. At a pressure of 

r 

S 00.000 atmospheres the formula yields — - = 0.-100, and if wc have 

* 0 

confidence enough in the formula to use it at p ~ 100 we find that 

V 

the corn-spending value of f- 0.295. In other words, a pressure 

* f 
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of 10 4 atmospheres will compress sodium to less than three-tenths its 


original volume 


The agreement betw een the values of — computed 


from the formula and the experimentally measured values, is within 3 
per cent over the entire range p = 2 5, 3, , 10 of the experiments 

The greatest discrepancy is at the pressure 2 5, where the calculated 


value of — is 0 815 and the measured value 0 789 This measured 

y o 

value is somewhat uncertain as it was obtained by an adjustment 
between two series of expenmen ts with different apparatus Over 
the range p = 4 to p — 10, for which direct readings from one appara- 
tus were obtained, the agreement between the values given by the 
formulas and the experimentally measured values is within 2 per cent 
(the discrepancies for all the values save the first being less than 1 per 
cent) 


V = 25 

3 

4 

5 

6 

7 

8 

9 

10 

y 









— (measured) 0 789 
Vo 

0 770 

0 737 

0 708 

0 683 

0 661 

0 641 

0623 

0 606 

~r (calculated) 0 815 

0 782 

0 750 

0 715 

0 688 

0 664 

0 642 

0 623 

0 606 









p » 11 

12 

13 

14 

15 16 

17 

18 

19 

20 


V 

— (calculated) 0 591 0 577 0 564 0 552 0 541 0 531 0 522 0 513 0 505 0 497 


Remark A noteworthy feature of the formula connecting p and 

V 

— which is furnished by the integrated linear theory is the unsym- 

V o 

metncal way in which it treats pressure (p > 0) and hydrostatic ten- 
sion (p < 0) Thus it is clear from the formula 

p _fia±lai {(&)*_!} 

that, as V — * 0, p — * « and, as V—*», p — » — ^2-lhJE^ In 

other words, the medium cannot, according to the integrated linear 
theory, support, without rupture, a hydrostatic tension of amount 

although it can support an arbitrarily large hydrostatic 
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pressure. From the formula in which the constants (for sodium) were 

determined by using the experimentally measured values of — at 

p, 5 and p : - 0 it would appear that the medium cannot support a 
hydrostatic tension of 1.851 X 10 4 atmospheres, and when the con- 
stants are calculated from the experimentally measured values of 

V 

-7- at ;>i - 0 and p2 - 10 this rupturing hydrostatic tension is raised 
Vo 

to 3.135 X 10 1 atmospheres. 

We now turn to the formula furnished by the second-order approxi- 
mation : 

p k { (3Xo + 2^o )c — (9Io 4* nolc"! 

«=-- (yj ' a ~ 2 *o, & — ~($h 4- no). 

/ F\ n 

Here 1 — 2c - ( J * a °d 80 ^ lG data on sodium may be presented 
ns follows: 

V - 2.5 3 -1 5 0 7 8 9 10 


V 

l”p 


0.780 0.770 0.737 0.70S 0.0S3 0.001 0.0-11 0.023 0.G00 


0.0731 0.0700 0.0020 0.102S 0.1122 0.1200 0.1283 0.1353 0.1-120 




310 2.750 3.013 4.450 5.2S5 0.09S 0.80S 7.087 8.462 


re have 


/ P V 4 

Denoting, for a moment, by g the product of p by f — J . we 

the following two equations to determine the two constants « and $ of 
our formula: 

aci 4- /Jor « q u 
or* 4* d r ; S ~ 72, 

> j and t* being the values of <• which correspond to any two values p\ 
and p*, rvv-]>cetively, of p. The formulas that serve to determine « and 
d are. tvreordinglv. 






< l:r 1 ~ 7iC; _ 
r if Aft ~ d) 
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In order to test the \alue of the second-order approximation for the 
purposes of prediction we again suppose that we are in possession of 
experimental measurements up to p = 6 only and w e use the formula to 
V 

predict the values of — corresponding top = 7, 8, 9, and 10 Since a 
V o 

and p are very sensitive to small variations in ci and e 2 if «2 — «i is 
small, we use the values pi = 2 5 and p 2 = 6 m determining a and 
P We find (check this calculation) that 

a = 2 619, P = 39G 3 


The results of the comparison between theory and experiment is 
furnished b\ the follow ing table 


p (observed) = 

25 

3 4 5 

1 6 

7 8 9 10 

Q *=■ P 1 y ) (observed) ** 

2 310 

2 750 3 613 4 450 

5 284 

6 098 6 898 7 687 8 462 

q (calculated) » 

2 309 

2 739 3 595 4 457 

5 282 

6 080 6 859 7 608 8 362 

p (calculated) = 

2 499 

2 9813 980 5 001 

5 90S 

C 980 7 956 8 908 9 882 


Note For reasons of convenience we have calculated the values of 
p that would correspond (when we use the formula of the second 


order approximation) to the experimentally measured values of — 

Vo 

rather than, as with the formula of the integrated linear theory, the 

V 

values of — that would correspond to the experimentally controlled 

V o 

values of p The greatest divergence between the calculated value of 
p and the experimentally controlled value of p occurs when the 
experimentally controlled value is 10 and the calculated value is less 
than this by less than 1 2 per cent Since the experiments do not claim 
an accuracy of more tiian 2 per cent, we may say that the formula of 
the second order approximation yields correctly the connection 
V 

between p and — up to p = 10, i e , up to a pressure almost twice 
' o 

as great as the larger of the two pressures that were used in determining 
the constants of the formula 

If we determine the constants a and p by means of the expen- 
mental measurements at pi = 2 5, pi = 10, instead of by means of the 
experimental measurements at pi = 2 5, p 2 = 6, we obtain the 
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following values: 

a - 1.902, 0 = *100.2. 

The comparison between the experimentally measured values and the 
values obtained from the formula 

p(f) - 1.902c + 400.2c 2 


is furnished by the following table: 


rvrtl) 

*2.r, 

3 4 5 G 7 8 9 

10 

fvY 

q - ;i l 1 (observed ) 

2.310 

2.750 3.G13 4.450 5.284 G.09S G.S9S 7.G87 

8.402 

q f<-:dcul:lt(sl) 

2.310 

2.74G 3.G15 4.489 5.327 G.13S G.931 7.094 

n 

p (calculated) 

2.500 

2.990 4.002 5.037 0.019 7.040 S.03S 9.008 

B 


'l’lnis the calculated values of p differ from the observed values of p 

by less than 1 per cent (i.e., by less than the experimental margin of 

error) over the entire range from p - 1 to p — 10. 

To find the value of p which would correspond to a given value, say, 

V /FY 5 

0.197 of ~~ wo proceed as follows: From ( — - ) — 0.G275 we obtain 

* o ^ \l'o/ 

r •- 0.1S02 and so p - M.*M, p — IS. 23 (as compared with the 

value p - 20 furnished by the extended linear theory). 

AT/-". The values 2.019 and 390.3 for a and /S. respectively, deter- 
mined from the experimental observations at pi ~ 2.5 and p 2 ~ 6 
differ from the values 1.902 and *100.2, respectively, determined from 
the measurements at- pi ~ 2.5 and p« — 10. This merely indicates 
that the assumption that q is a quadratic function of c over the range 
n V < 10 is only an approximation to the truth; the parabolic 
approximation over the range 2.5 ^ p 0 differs slightly from the 
pimholtc approximation over the range 2.5 Sj p ^ 10. The fact that 
3 ]- '-o large compared with a indicates the importance of considering 
the third-order elastic constants. Whereas X 0 -f is around 0.7, 9/ 
*:* t; is negative and around — -100. 
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In order to test the value of the second order approximation for the 
purposes of prediction we again suppose that we are m possession of 
experimental measurements up to p — 6 only and w e use the formula to 
V 

predict the values of — corresponding top — 7, 8, 9, and 10 Since a 
Vo 

and 0 are verj sensitive to small variations m ej and e 2 if «2 — ei is 
small, we use the values pi = 2 5 and pi = 0 in determining or and 
0 We find (check this calculation) that 


a = 2 619, 0 = 39G 3 


The results of the comparison between theory and experiment is 
furnished by the following table 


p(obser\ed) ** 

25 

3 4 5 

C 

7 8 9 10 

(v\* 

<? - P l pr I (observed) ~ 

2310 

2 750 3 613 4 45C 

5 284 

6 098 6 898 7 687 8 462 

9(calculnted) = 

2 309 

2 739 3 595 4 457 

5 282 

6 080 6 859 7 608 8 362 

p (calculated) = 

2 499 

2 981 3 980 5 001 

5 998 

6 980 7 956 8 908 9 882 


Note For reasons of convenience vve have calculated the values of 
p that would correspond (when we use the formula of the second 

V 

order approximation) to the experimentally measured values of — 

rather than, as with the formula of the integrated linear theory, the 
V 

values of — that would correspond to the expen mentally controlled 

M 0 

values of p The greatest divergence between the calculated value of 
p and the expenmentally controlled value of p occurs when the 
experimentally controlled value is 10 and the calculated value is less 
than this by less than 1 2 per cent Since the experiments do not claim 
an accuracy of more than 2 per cent we may say that the formula of 
the second order approximation yields correctly the connection 
V 

between p and — up to p = 10, l e , up to a pressure almost twice 
'0 

as great as the larger of the two pressures that were used in determining 
the constants of the formula 

If we determine the constants a and 0 by means of the expen 
mental measurements at pt = 2 5 p 2 = 10, instead of by means of the 
experimental measurements at p 2 = 2 5, p% = 6, we obtain the 
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following values: 

ft = 1.902, 0 = 400.2. 

The comparison between the experimentally measured values and the 
values obtained from the formula 


p(~rj = 1.902c + 400.2c 2 
E furnished by the following table: 



2.5 

3 4 5 0 7 S 9 

10 

/i*V* 




'1 ” ;> tyr ) (oh-erve<l i 

2.310 

2.750 3.013 4.15G 5.2S1 0.09S 6.89S 7.087 

8.462 

‘i friltulnU'd) 

2.310 

2.740 3.015 4.489 5.327 0.13S 0.931 7.091 

8.462 

p (f'llolllsltf tl) 

2.500 

2.000 *1.002 5 037 6.019 7.016 S.038 9.00S 

10.00 


Thu« tltt* calculated values of p differ from the observed values of p 
by less than 1 per cent (i.e., by less than the experimental margin of 
error') over the entire range from p = 1 to p = 10. 

To find the value of p which would correspond to a given value say 
V /y\» 

0.407 of ™ uo proceed as follows: From ( — j = 0.G275 we obtain 

r • > (USG2 and so p = 14.44, p - 1S.2.3 (as compared with the 


value - 20 furnished by the extended linear theory). 

Xo!l ■ Tll, ‘ vftluw 2.010 and 300.3 for « and respectively, deter- 
mined from the experimental observations at p, = 2.5 andV = g 
differ from the value-' 1.002 and 400.2, respectively, determined from 
the measurements at />, - 2.5 and p 8 - 10. This merely indicates 
that the a-umptton that q is a quadratic function of e over the ranee 
0 p h 10 E only an approximation to the truth; the parabolic 
appear, mat,,,, oyer the range 2.5 < p ^ G differs slightly from the 

mahnn 0VCr ^ range 2.5 -< p .< 10. The fact tint 
compared with a indicates the importance of considering 

tn- t lard-order ola-tie constants Whereas X -x. ? ; , „ ng 

, • ’ " ncrcas Aq -f -jpo is around 0 7 0 / 

— f. <• negative and around —400 ' ’ y 
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1 . Definition of a non-isotropic elastic medium 
Certain materials (for example, wood) that are important m con 
struction are not isotropic i e , they are sensitive to at least some rota 
tions of the medium (before application of the strain) We have seen 
also that even if a medium is isotropic when it is initially free from 
stress or when the initial stress is a hydrostatic pressure or tension, it 
cannot be isotropic when the initial stress is any stress other than such 
a hydrostatic pressure or tension It is, therefore, highly important, 
from the scientific point of view, to obtain the relation connecting 
stress and strain for non-tsotropic media Unfortunately, however, 
this relation is much more complicated than for isotropic media, even 
m the infinitesimal theory many more elastic constants than the two, 
X and a that suffice for isotropic media are needed In the second 
order approximation the number of elastic constants becomes, usually, 
so large that the theory is too complicated for practical applications, 
but it is possible to construct an integrated linear theory that is useful in 
practice Many deformable media have a crystalline structure so that, 
although they are not isotropic, 1 e , elastically insensitive to all pre- 
lotations of the medium the} are elastically insensitive to certain such 
prerotations Thus the relation 

<t>(R*riR) = <f>(y) 

will be now valid for certain given rotation matrices R but not for 
all rotation matrices When <t> is analyzed into a sum 

= 4>o + 4>i $2 + 4>s + 

of terms of different degrees in the elements of ij (<pj being a homo- 
geneous function of degree j in the elements of v, ] ** 0, 1, 2, ), 

the relation 4>{R*yR) = <#>(77) implies the senes of relations 

4 j(R*vR) = Mv), 3 - 0, I, 2 

(Prove this Hint The relation <*>(R*yR) = is an identity in 
the elements of 17 It remains true, then, if 17 is replaced by ky where k 
80 
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i? an arbitrary parameter; when this substitution is made 7) is 
replaced by ?).) The relation - $j(v) imposes certain 

conditions on the coefficients of the homogeneous function, £,-(??), of 
degree j, of the elements of 77 (save in the case j — 0 in wliich we are not 
interested [whv?j) and we proceed to find the nature of these condi- 
tions, for various given rotation matrices /?, when j - I or j — 2 or 
j 3, 

EXERCISES 


1. Show that if a deformable medium is elastically insensitive to n given rotntion 
R it is elastirnlly insensitive to every integral power, positive or negative, of /?. 
Hint. Iteplnee rj liy Il*r}Jt anti by RyU* in the relation $(R*rjR) ■= 0(17). 

2. Show that if the 1111x111110 is elastically insensitive to a given rotation R it is 
rltiMtrnlly insensitive to the reflexion —it. Sole. A reflexion is an orthogonal 
Irntuforitmtinn of rectangular Cartesian coordinates the determinant of whose 
matrix is —1 {instead of +1 as for a rotation). 

3. Show that if II is a rotation through x around the n-nxis, —7? is a reflexion in 
the n- plane. 

-i. Show that if the medium is elastically insensitive to two given rotations 
7’iand R jit is elastically insensitive to their product R\Rz. 

5. Show that if the medium is elastically insensitive to a rotntion through x 
around the maxis and to a rotation through x around the 6-nxis it is elastically 
itvnudtive to a rotation through x nrotmd the r-axis. Hint. Any one of these 
rotation" is the product of the othrr two. 


Show that if the medium is elastically insensitive to n rotation through 


2 


around the o-nxi- nml to a rotation through ^ around the 6-axis it is elastically 


ij>*rn*Utvo to a rotation through ^ around the r-nxi.«. Hint. 


Denoting these 


rotation*, by }{■>, and R 3 . respectively, It* «■=* R i/?;6’i *. .Vo te. Certain eubie 
crystals jkivo >s the elastic symmetry of this exercise. 

!'!u*w that if a medium possesses the elastic symmetry described in Kxerci“e 


. .... 2x 

® >* i" rlaatieMly insensitive to n rotation through •— around each of the four 
diagonals of the rube whose* vertices are the points (±1, ±1, ±1). Hint. The 


product of a rotation through ~ around the 6-axis by a rotation through ~ around 
tin- r-nxi" a rotation through ~~ around one of the diagonals of the cube in 

o 


2. Hit form of £t(rj) for various rotations R 

is n linear function of the elements of r,. When written sym- 
metrically (what does this moan?) it appears in the form 


*i(“> - Tr(C t v) 
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■where C\ is a symmetric 3X3 matrix The elements of C\ are the 
values, when 77 = 0, of R*TR for 

R*TR « Af M ~ M li 

and when i) — 0, — = 1, M = E 3 , and ~ — C\ Thus, 

Pa 017 

The coefficients of the linear function <fn(v) of the elements of y, when 
this linear function is mitten symmetrically, are the coordinates 0/ the 
initial stress in the reference frame in which JJ* ~ J*J On denoting 
the initial stress matrix by To ue hate 

R'ToR = Ci, R - JM~* 

We non denote the variables (ji aa , mb, n <*, Vbc, Vca, Vab ) by (171, m 
V3t V*, Vs, Vi) and write ^i(tj) in the form 

4>l(v) ™ Cl»?l + C 2172 + O3VZ + C4?74 *f CbV S + Cfi77« 

Thus the symmetnc 3X3 matrix C\ is 

( ci fa* 

c 2 ic* I 

ic 5 £c« C 3 ) 

Under the rotation a— » a’ = R*a, 17 — * 77' = R*yR, and since aa*-* 
aV* = R*aa*R it follows that the elements of 77 transform, under the 
rotation o — * a' = R*a, like the elements of aa* This fact implies 
(show this) that (*71, tj 2 , 173, v* vs, ns) transform like (a 2 , b 2 , c 2 , be, ca,ab) 
Let us consider the case in which R is a rotation through an angle 8 
around the a-axis, since the a axis is arbitrary, this takes care of every 
rotation There is no lack of generality in taking 0 < 0 < v, since R 
may be replaced by R * and we shall do so We have 

/I ° 0 \ 

R = I 0 cos 6 — sin 6 1 
\0 sin 0 cos 6/ 

so that a' = a, 6' = (cos 6)b + (sin 0)c, c' = -(sin 6)b + (cos 6)c 
Hence 77 — ♦ 17' where 

VI - Vu Vi' - (cos 2 0)772 + 2(cos 8 sin 0)174 + (sm 2 8)vi, 

773' *= (sin 2 0)772 ~ 2 (cos 0 sm 0)774 + (cos 2 0)773, 

774' «= -(cos 0 sm 0)772 -f (cos 2 0 — sm 2 0)774 + (cos 0 sin 0)i?3, 
vs' = (cos 0)175 — (sm 0)774, ij 6 ' = (em 0)775 + (cos 0)vi 
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Sinro r,* -f ra ~ vz + *73 it is convenient to introduce, ns replacing 
the variables (nu * * * , vc), the variables (fi, £?, £3, £4, £5, £c) where 

Ii " vu tz = tivz 4- 7};), £s = - *)n), 

£4 = *74. £.'> - *75. £c — *76. 

Then Tjt ~ U. *72 - £: + £n, *73 = £2 ~ £ 3 . *74 = $4. *?r, « $5, Vo ~ to 

t»nd, under the rotation n a' — Ji*n, £ — £' where 

si* si. £; r ~ tz, tz ~ (COS 2fl)$:s + (sin 20)t-\, 

£<' “(sin 20) £3 + (cos 20)£.,, £ f ,' = (cos 0)f 5 - (sin 0)£ fi , 

£c' = (sin 0)£ 5 + (cos 0)£o- 

These relations suggest the introduction of the (complex) variables 
(ft. fs, fr„ f 4, fs, ft) defined by the equations 

fi C! fit f; ~ fit S' 3 - 4(^3 — i£t), 

it ~ {.(£?. H- I $ 4 ), f r, “ 4(fr, — Ur,). ft = 4(f.i 4- *£(•>)• 


Of the six variables (fi, • • • , ft), fi and fn are real, and fj and fi are 
conjugate complex numbers as arc also f& and f c . It, is dear that 

Ii ~ ft, £; ~ f;, £3 ~ fs 4 s' 4 , £4 - >'(lz ~ ft), £5 = fs + fe, £c = 
if f j - ft) and it follows (show this) that, under the rotation a — * a' — 
ll*n, f -- f' whore 

ft' fi, fs' - fs, fs' = c 2i %, U' = c" !<5 f4, 

fs' = c-^fs, f r / = e«f«. 

Tim relations connecting the original (real) variables (iji - • • 17c) with 
the new (complex) variables (fi • • • f c ) are 

“« £1 ft; fi - fi - *71 : 

r i* 4 £3 ~ f: 4 - f j 4 f 4 ; fs *= £2 = id* 7 : 4 - *73); 


it - £3 - f: - fa ~ IV, f 3 = “ (£s - i£i) ~ “ (*72 *73) - ~ >74; 


I 


tt - 


Ufa - u); 


f< - o (£3 4* t’£t) = - (*7-> — *73) + - *7i : 


9. • 2 

fs = 4(£t - i£«) = 4(*7i — i**7e); 

•t *'• %s *= »(iE st); ft = 4(£t 4" j£g^ r - 4(*?;. 4 **?d. 

V,e nenote the coefficients of when this expression is written as a 

unction of the variables (f,. ■ • • , f f ). i,y (,/ t , • • • , ,? c ) and so 


** C7 i» “ h 4- fc; 
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•Alfa) = ClVl + C2V2 + C3173 + C4IJ4 + C5JJ5 -f- C6VB — d if 1 -f* dfo -f- dsfs -f 

di an< ^ d 2 are real, and d 3 and ^4 are conjugate 
complex numbers as are also d 3 and d$ The necessary and sufficient 
condition that the medium should be elastically insensitive to the rota- 
tion R, as far as #1 is concerned, is that the relation 

difi' 4 - 4 * dfiiV = difi 4 - + d 6 fe 

should be an identity in the six variables (f 1, , f 6 ) Upon sub- 
stituting for (fi', , fg') their expressions in terms of (fi, , 

fe)) we obtain the following relations involving (dj, , rf 6 ) and 
the angle 0 of the rotation R 

e ix, d 3 = d 3 , c^'di = d 5 

Note In w ritmg these equations \\ e have put dow n only one of each 
pair of conjugate complex equations, thus the equation e txl d 3 = d 3 
implies the equation e~ 2xl (U *= d 4, and the equation e -,9 d fi «= d 6 implies 
the equation e x, dt - d 6 

If 6 is not the quotient of 2 ir by 2 , 1 e , if 9 7^ ir, neither of the num- 
bers e 2xi , e~~* 9 is unity and we have ^3 = 0 d 3 = 0 The medium is 
then elastically insensitive, as far as <(n is concerned, to a rotation 
through any angle around the a axis (why?) and 0 i(jj) is a linear 
combination of fi = vi and f 2 = %{vi + 17s) We have, then, the 
following result 

If a deformable medium is elastically insensitue to a rotation through 
an angle 6 7* tt around the a-axis it is, as far as <f>i is concerned, elasti- 
cally insensitive to a rotation through any angle around the a-axis and 
^1(17) is a linear combination of vi and ij 2 + 173 

0l(»») = ClTJl + C 2 (lj 2 + vs) 

Thus the medium cannot be elastically insensitive to such a rotation 
unless the initial stress matrix is, in the reference frame in which JJ* = 
J*J, diagonal 

EXERCISES 

1 What 13 the necessary and sufficient condition that the medium be elastically 
insensitive as far as £1 is concerned, to a rotation through any angle around the 
5-axis? around the c-axis? Hint When the (a b c)-axes are subjected to ft 
cycbc permutation the two sets (123) and (456) are subjected to the same cyclic 
permutation 

2 Show that if the medium is elastically insensitive as far as is concerned to 

a rotation through any angle around the a axis and to a rotation through any 
angle around the b axis it is as far as 4>i is concerned, isotropic Hint is °f 

the form c(i» + ij 2 + ijj) *= c h Hole The phrase ' a3 far as <tn 19 concerned’ 
may be omitted since every rotation may be factored into the product of three 
rotations, two of which are about the a axis the third being about the b axis 
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If a ~ r, c ,it ~ 1 and d 3 is arbitrary (instead of being zero). £1(7?) 
K now a linear combination of fi» f;, fj, and £4 or, equivalently, of 
*jl. V:, Vs, VE- 

4 >l(v) — ClV\ *f C2V2 + C3»?3 + c -*V-\- 

The initial stress matrix (in the reference frame in which ././* = J *,F) 
is of the form 

/ci 0 0 \ 

J 0 C; ict 1. 

\0 Ac, r.3 / 

The necessary and sufficient condition that the medium should be 
elastically insensitive, ns far as <£] is concerned, to n rotation through - 
around the o-n\is or, equivalently (why?), to a reflexion in the n-plane, 
is that ^i(t)) — Cjt;i -f* C 2V2 + taija + 04774. 

EXERCISES 

3. What is the neefs-.ary anil sufficient condition that tl>e medium l>e elastically 
ia-'-n- iti\ e, as far ns is concerned, to a reflexion in the f>-plnne? in the oplnne7 
-J. Show (hat if the medium is elastically insensitive, as far as $ 1 is concerned, 
to a relation through ntt angle 0 <■/ r around the o-avis it is elastically insensitive, 
ft* fnras is concerned, to a reflexion in each of the coordinate planes. 

3. The form of <>;(ij) for various rotations R 

~<i>:{v) is a homogeneous quadratic function of the elements of 77, and 
so it may he written in the form 


r ~ (ci]7)V cr,r,7c") + 2{<'iiViV2 cr.cvsvo)- 

Tlie medium is elastically insensitive, ns far as <f> 2 is concerned, to a 
given rotation R if, and only if, the relation 
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independent, in this e\ ent we say that the medium possesses fewer 
than twenty-one second order elastic constants An isotropic medium 
possesses only two second-order elastic constants, the degree of com- 
plexity introduced into the theory of elasticity by lack of isotropy may 
be grasped when we compare this with the twenty -one second-order 
elastic constants of a medium that is completely non isotropic (as 
far as is concerned) 

We now consider the case m which R is a rotation through an angle 
0 around the a-axis, and we introduce the complex variables (f 1( , 

f#) of the preceding section We denote the coefficients of 2<fo(’?), 
when written as a function of (fi, , £#), by d vq 

= diifi* + + d««r« 2 + 2(di 2 fif 2 -f + dssfife) 

Here du and d 22 are real, and djj and d<4 are conjugate complex num- 
bers as are also da and d 66, dj 2 , dj 4 , and djg arc real, dn and du arc 
conjugate complex numbers, etc The necessary and sufficient condi- 
tion that the medium be elastically insensitive, as far as is con- 
cerned, to the rotation R is that the relation 

du(fi') 2 + + d«(fV) 2 + 2(du$\ti + + dsafsT#') 

= dnfi 3 4- + dee?# 2 + 2(d l2 f!f 2 + + <Wsfe) 

should be an identity in the six (complex) \ anables (f lf , tt) 
This furnishes the following relations involung the coefficients d p9 and 
the angle 9 of the rotation R 

e* l, dzz - dzz, e‘ lt d\z = du, e“**d 2 s — djs, 

c _ 2 ‘ , d S s = da, — du, e lt dj 5 = djs, 

e 2 ,, d 2 3 = d 2 j, e Sl, d36 = dj# 

Note In WTiting these equations we ha\o put down only one of each 
pair of conjugate complex equations, thus the equation e*'*dzz = dzz 
implies the equation e~‘ , *d44 = d«i etc 

If 9 is not the quotient of 2ir by 2, 3, or 4 none of the multipliers 
e <lt , c ±: ' t , c ±lt , e ut is unity and so all the (complex) numbers dss, da, 
diz, dis, </« 3 , d 2 5 » dss, dst are zero In this event the medium is elastic- 
ally msensitne to aery rotation around the a axis and it possesses at 
most 21 — 2(8) = 5 second-order clastic constants, # 2 (>j) being a 
linear combination of fi 2 , f 2 2 , fjf 2 , fjf 4 , and fjfg or, cqimalently, of 
iti 2 i (*J 2 + >»)*, vi(vt + Vi), (.Vi — Vi) 2 + 4 »I 4 2 , and ij 5 a 4- Vt‘ Since 
»?i(’72 + is) *» /|i» - T 7 i*r hvi may be substituted for rn(t/t + vt) Jn 
this linear combination, and since ij 2 -f- »?j «= I t — nu / i a may be sub- 
stituted for (173 + 173) 2 Since (ij 2 — >jj) 2 + 4 iji 2 *= (172 + »») 2 ~ 
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— m~) f nina ~ Vi 2 tnay be substituted for (r lZ — gf) 2 + ‘hf 2 , 
and, finally, since / 2 ~ (t?^ — V4 2 ) + vi(gs + m) ~ (*7s“ 4* tfc 2 )) /s 
m.nv be substituted for n 2 + vc 2 - We have, then, the following result : 
// u deformable medium is elastically insensitive (as far as $2 is con- 

2tt‘ t 

court!) to a rotation through an angle other than or - around the 

o Z 

rt-oxjV, 1 / f« elastically insensitive (as far as $2 is concerned ) to evert/ 
rotation around the a -nr is. Such a medium possesses not more than 
foe (instead of hrtnt’j-onr ) second-order clastic constants, and <t>z(y) is a 
linear combination of (he following five functions: 

1 1 2 - I;, ’ll 2 , h i’ll, ~ V4 2 - 
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m> V2, vz, >J4 (involving ten coefficients) and a linear combination of m*, vt~, and 
v m 

7. What is the form of <t>t(v) if the medium is elastically insensitive to a rotation 
through x around the 6 axis (or, equivalently, to a reflexion in the 6 plane)? to a 
rotation through x around the c axis (or, equivalent!} , to a reflexion m the c plane) ? 

8 Show that if the medium is elastically insensitive to a rotation through r 
around the a axis and to a rotation through x around the b axis it possesses not 
more than nine second-order elastic constants, $j(jj) being the sum of a general 
polynomial of the second degree in the three variables iji, 175, jjj (involving six 
coefficients) and a linear combination of vt 1 , i?s 2 , and ije 2 

When 6 = the multiplier e 3,t is unity and so d 3 e is arbitrary, 

instead of being zero The medium possesses (at most) 5 + 2 = 7 
second order elastic constants, and <£2(17) is a linear combination of the 
five functions 1 1 2 , /j, »ji 2 , Ivtu 172173 ~ 174 2 and, in addition, the two 
functions (172 — 173)175 + 2174176, (772 — 173)176 — 2*74175 (which are, respec 
tively, the real and imaginary parts of 8f s f g ) 

EXERCISE 

9. Show that if a medium which possesses the elastic symmetry described in 
Exercise 8 is elastically insensitive to a rotation through — around the line joining 

the points (0, 0, 0) and (1, 1, 1) it possesses only three second-order elastic constants 
and that for such a medium *2(17) is a linear combination of Ii a , It and (ujijj + 
ijjiji + ijiija) Hint 4 >t(.v) is invariant under the transformation ij-* ij' = R*vR 
where 



thus is invariant under the transformation ijj' = yi vt' = i|a> vs' m VU 

tji' = nr,, vt = Vt Vt = Vt Note Certain cubic crystals possess the elastic 
symmetry of this exercise 

When 6 = e M is unity and so d 3 3 is arbitrary, instead of being 

zero The medium possesses (at most) 5 + 2 = 7 second order elastic 
constants, and ^2(17) is a linear combination of the five functions I i 2 , 
1 2, i?i 2 » 1 117 1, 17217s — 17 4 2 and, m addition, the two functions (172 — 173)* — 
4 i? 4 2 , (132 — 173)94 (which are, respectively, convenient multiples of the 
real and imaginary parts of fa 2 ) It is easy to see that ./217s and 174 2 are 



they may be substituted for 172173 — 174 2 and (172 — 173) 2 — 4 174 s Thus 
we have the follow mg result 
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The medium is elastically insensitive (as far as is concerned ) to a 

rotation through ~ around the a-axis if, and only if, is a linear com- 
*■* 

bination of the following seven functions: 

1 1 2 , Is, yi~, I mu vim, ’jr. (v: — nz)iu- 
EXERCISES 
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fa ~ i{(Clim 5 + ) + 3(CH2>J1 2 >J2 -f- ) 

+ 6(ci2S»?Xlj2T73 + )} 

There are fifty-six coefficients c vqr , p ^ q ^ r, m all six of the type 
cm, thirty of the type Cm, twenty of the type c 12 a Thus if fa(>j) is 
sensitive to every three-dimensional rotation R the medium possesses 
(as far as fa is concerned) fifty-six elastic constants, we term these 
third-order elastic constants A medium that is completely non iso- 
tropic as far as fa and fa are concerned possesses, then, 21 + 56 = 
77 clastic constants of which twenty-one are second-order constants 
(i e , coefficients of fa) and fifty-six are third-order constants (i e , 
coefficients of fa) The six coefficients of fa (or first-order elastic 
constants) are the coordinates (in the reference frame in which JJ* - 
J*J) of the initial stress We proceed to investigate what reduction is 
caused in the number of the third-order elastic constants when the 
medium is elastically insensitive to one or more given rotations 
Introducing the complex variables (ft, , te) of the preceding 
sections, we denote the coefficients of 3fa(y), when this is written as a 
function of the variables (ft, , f 6 ), by d pqr , p ^ q ^ r = 1, 
, 6 

3^3(17) = dmfi 6 + + 3 (dn 2 fi 2 f 2 + ) 

+ 6(dt*ifiJ - 2fj + ) 

Of the fifty-six coefficients d pqr eight are real (why?) — dm, dm dm, 
di22, di34, d 1 5 6 , dm, dm — and the rest are pairs of conjugate complex 
numbers (why?) 


(da 33, 

dm), 

(d55S, dggg), 

(dn3, dm), 

(djis, dug), 

(d 2 23, 

dm)t 

(d 2 25i d 2 2 g), 

(dm, dm), 

(d 2 33, d 2 44), 

(d334, 

dm), 

(dsss dm), 

(ds 38, d445), 

(digs, digg), 

(d: 55, 

due), 

(dm, dm), 

(d45j, dm), 

(dsss, dfi6«), 

(dl23, 

dm), 

(dl25, dl2g), 

(dl35, dj 4g) , 

(dl36, dj45), 

(da 3 5, 

d 2 4 b) , 

(d236, d245)p 

(ds45, d 3 4g), 

(d35g d4se) 


The necessary and sufficient condition that the medium be elasticallj 
insensitive, as far as 4 > 3 is concerned, to the rotation R through an angle 
0 around the a-axis is that the following relation should be an identity 
in the six complex variables (fi, , fg) 

dmfi' 8 + + 3(dix2fi' 2 f 2 y + ) -f G(tfj23f iWfV + ) 

= diufi 3 + + 3(d lls fi J rj + ) + C(d ls3 niSh + > 
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EXERCISES 

1 Verify that each of the ten functions written above is insensitive to a rotation, 
through an arbitrary angle 6 about the a axis Hint ,j = f i, 12 + u “ 2 f 2 
(12 — ,3)* + 4,4 s = 16 f 3 f« i?6 4 + vt 1 = 4fjfe (,2 — ,3) (,5 s — ,g s ) — 4,417610 ** 
real part <7/ 1 8?tfs a , and in — nltste + nils 2 — ne 2 ) — imaginary part ol 
— Sfafs* 

2 Show that ,21a — H 2 is insensitive to a rotation through an arbitrary angle 8 
around the a axis Hint 4 (,1,3 — 14 *) = (12 + ,j) s — |(,i — 13) 2 + 4,4 2 j 

3 Show that ,216* + , 3 i & 2 + 2,4,5,( is insensitive to a rotation through an 
arbitrary angle 0 around the a axis Hint 2 (,2,8 s + ,j,g 2 + 2,4, 6 , e ) = 
(12 + ,j)(is 2 + io s ) — l(,i — 13) (1 6 * — n 2 ) ~ 4 , 4 , 8 , j | 

4 Show that /» = ,1(12,2 — ,4 s ) + (12,0 s + 13,5 s + 2,4isio) — (,2 + u)(,5 S 
+ 16 1 ) 


It follows from the result of Exercise 2 that ,213 — *?i 2 may be used 
as a substitute for (172 — 173) 2 + 4,4 2 and from the results of Exercises 
3 and 4 that 7 3 may be used as a substitute for (,2 — ij 3 ) (, 6 2 — 179 2 ) — 
47*1751-9 /i(ij2»j3 — n* 2 ) may be used as a substitute for (172 + vz)(vm 

— Vi 2 ) (why*), and Ii(i75 2 + 7 s 2 ) may be used as a substitute for 
(V2 + I7s)(i75 2 + VC 2 ) (why?) Since IJ2 - Iifav* ~ V 2 ) ~ h(vi 2 
+ vc 2 ) + Iivibn + vz) = I 1(7273 - v 2 ) ~ h(v s 2 + vc 2 ) + 7i 2 (72 + 
vz) + 171(172 + i7s) 2 , it follows that 7 i 7 2 may be used as a substitute for 
71(72 + 1J3) 2 and, since I'm 2 = 71 3 + 7i 2 (72 + 173), that 7m 2 may be 
used as a substitute for ,i 2 (i72 + 173) Finally, since 7i 3 = ,i 3 + 
3 »fi 2 0?2 + 173) 4 - 3 i?i(i ?2 + 173) 2 + (172 + vz)*, Ii* may be used as a sub- 
stitute for (172 -j- 173) 3 Hence we have the following result 

A deformable medium is elastically insensitive, as far as ^3(17) is con- 
cerned, to an arbitrary rotation around the a axis if, and only if, <t> 3(17) is 
a linear combination of the following ten functions of (171, , 17a) 

7l 3 , /i/ 2 , I3, Vl*> 7m 2 , 17l(l72»73 — 174 2 ), 7i(l72,3 — »74 2 )j 

Vl(vc 2 + VC 2 ), 7 j(i 75 2 + VC 2 ), VbVcivi ~ 7 a) + 174 ( 175 * — »76 ) 


EXERCISE 

5 What is the form of £3(1) if the medium is elastically insensitive to every 
rotation about the b-axis? about the c axis? 


When the angle 8 of the rotation around the a-axis to which the 
medium is elastically insensitive is the quotient of 2 tt by 2, 3, 4, 5, or 6, 
some of the multipliers e M , , e~ 2lt are unity and so some of the 
twenty three coefficients d 33 3, , d«6 are arbitrary complex num- 

bers (instead of having to be zero) The medium may possess, then, 
more than the ten third order elastic constants to which it would be 


limited if 6 were different from t, or ^ The results for these 

3 2 5 3 

various values are as follows (check these statements) 
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elastic constants, and <£3(1?) is a linear combination of the ten functions 
already written down and of 10 new functions that may be taken to 
be the real and imaginary parts of fa 3 , fs 3 , f if sfe, and (or 

any convenient ten linearly independent combinations of these) 

When 6 — -> the multipliers e 4 '® and e~ us each are unity, and so the 

three complex numbers dm, d 2 33, and d 4 &5 are arbitrary (instead of 
being zero) The medium possesses (at most) 10 -f- 2(3) *= 16 third- 
order elastic constants, and $3(17) is a linear combination of the ten 
functions of (171, , v e) that we have already written down and of 

six new functions that may be taken to be real and imaginary parts of 
f if 3 2 , f 2 f 3 2 , f 4f 5 2 However, the form of <fo(j?) is most easily derived as 
follows Since the medium is elastically insensitive to a rotation 
through it around the a axis (why 7), 3^3(17) is of the form 

3^3(1?) " cm»ji 3 + c 2 22V2 3 4- 0333173* 4- 0444174* 

4- 3(cii2i?i 2 i72 + 4* C344173174 2 ) 4- 6(0123171172*13 4- 

4* 0234*72*13114) 4- 3(Ci55l7l 4- C2S5*?2 4- C355173 4- C455*14)>15 2 
4- 6(ci56*7l 4- 0258*12 4" 0356*73 + C458*74)*l5l78 

4- 3(ci66*ll 4- 0266*12 + C3 6 6»13 4" C466*l4)*?8 2 

Under the rotation through ^ around the a axis, namely, 

a' = a, b' ~ c, c' = — b, 

v—nj' where i?i' = ijt, *12' = >13, Vs' - *72, v*' *= — 174, *7s' = — * 18 , 
*16' = *75, and since the relation <f> 3(17') = <f> 3(17) must be an identity m 
the six vanables (171, , rj 6 ) we obtain the following relations 

C114 = Cl56 — C234 = C444 *= 0, C112 — C113, C122 — C133, 

Cl24 = C1S4, Cl 5 5 = 0166, C222 = C333, C22S = C 2 33, 

0224 = C334, 0244 = C344, C2S5 = C366, C 25G — —0356, 

0268 ~ O355, C455 = OlBt 

Thus 3^3(17) is of the form 

3^3(17) = Clll»li 3 4- C222(*72 S 4- *13 3 ) 4- 3{cii2?7i 2 (l72 4" 173) 

4- c 122*7l(*32 2 4* *73 2 ) 4- C22 3*72*1 3(*72 4* *13) 4" C224*74(*l2 2 — »?3 2 ) 

4- 0144*11*74* 4* C 2 44*74 2 (*?2 4" *73)} 4* 6 1 Cl 23*Jl'72*l 3 
4* Cl24’7l*7 4 (’72 — *13) 4- C 2 56*15*76(112 ~ *73) 4- 0 4 56*l4’75’l6 } 

4" 3{ci56*7l(*J5 2 4* *16 2 ) 4- 0255(172*75* 4* *l3*?6 2 ) 

4* C266(»>2’76 2 4" 11317S 2 ) 4" 0455174(175* — 1J6 )) 
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1. The modified stress matrix 

We have seen that when a deformable medium is in a state of equi- 
librium in a strained position the stress matrix T is connected with the 
applied force per unit mass matrix F by means of the relation 
(div* T)* -f p z F — 0 

and with the applied force-per-umt-area matrix / by means of the 
relation 

TdS* = dS x f 

( dS z being the matrix element of area, and dS z the scalar element of 
area, of the surface of the medium when deformed) Since the surface 
of the medium when deformed is not given (what is given being the 
surface of the medium when undeformed), we introduce the 3 X 3 
matrix T a , which is such that 

T a dS a = TdS z 


We term this 3X3 matrix T a the modified stress matrix Since T = 


(V)/- 

\pj dr\ 


J* and since dS x * 


) (J*) 1 dS a (see Exercise 1, p 17), 


_ T d<f> 
T a = J-r 
di] 


Note Although the stress matrix T is symmetric the modified stress 
matrix T a is not, in general, symmetric T a wall be symmetric if J 
is symmetric and if the medium is isotropic, for, then, J commutes 

with j? = i(/ 2 — Ez ) and, hence, with — 

dri 

The relation that connects the modified stress matrix T a with the 
applied force-per-umt-mass matnx (replacing the relation (div* T)* + 
p x F = 0) is readily found In fact, the virtual work of all the forces 
96 
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ndiiiR on nny portion of the medium in any virtual deformation is 
furnished hv the formula 

Virtual work = f p.(S.r)*/-’ dV x 4- j (6r)*f dS z 

*/ i * v 

- £ Pn {bx)*FdV a + (bx)*T a dS a 

= f \ Pn (ox)*F + div a (T n *5r)\ dV a 
Jr. 

M 

where, if p « I 7 1 i« any 3 X 1 matrix whose elements arc differ- 
entiable functions of (a, b, c), we understand by div„ p the sum p a + 
q K 4- r f . Similarly, if P is any 3 X m matrix we understand by div„ P 
the 1 X m matrix whose elements are obtained by applying the 
operation div„ to the various columns of P. 

This virtual work must be zero for the virtual translations (for which 
or is a constant 3 X 1 matrix) for every portion V„ of the medium, and 
so 

(div„ r n *y 4- PnF = 0. 

On substituting for T„ its value J so that T a * — ( — ) J* we 

dy \dy/ 

obtain the basic equations of equilibrium of the medium in a form 
suitable for employing the initial coordinates (o, b, c ), rather than the 
final coordinates {r, y, :), of a typical particle of the medium as inde- 
pendent variables: 




4- n,F « 0. 


When the medium is isotropic (so that — is symmetric^ and the 

\ ^ J 

applied mass force is either absent or negligible, these equations 
simplify to 


or, equivalently, 


dn-„^.r 

fir; 


- 0 , 


div rt --- 0. 


•' r ' f >■ In the mtimtesjinnl theory of elasticity, — replaced bv 
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X/i + 2/xij (the initial stress being zero), and the product of any ele- 
3 <> 

meat of — by the derivative of any element of J* frith respect to a, b 

dij 

or c is neglected The equation div a J*^j = 0 reduces, then, to 

^div a J* = 0 or, equivalently (since J* is non-singular), to 
div» (X/ x + 2firj) = 0 


If we assume the medium to be homogeneous, so that X and p are 
constant functions of (a, b, c), this equation yields 

X(/i). + 2n div 0 >i = 0 

Note 2 In the derivation of the basic equation 
(div B T a *)* + Pa F = 0 

we have assumed that the final reference frame is the same for every 
point of the medium (so that the virtual translations are characterized 
by the constancy of the 3X1 matrix 8x ) It sometimes happens (for 
example, when we wish to use space polar or cylindrical coordinates) 
that the directions of the axes of the final reference frame vary from 
point to point of the medium Let Ri be the rotation matrix that 
transforms any convenient fixed reference frame to the (moving) final 
reference frame If <5£ is the 3 X 1 matrix whose elements furnish 
the coordinates in the fixed reference frame of the vector whose 
coordinates in the moving final reference frame are furmshed by the 
elements of the 3X1 matnx 8x, we have 8x = Ri*8£, and the virtual 
translations are characterized by the constancy of the 3X1 matnx 8( 
The basic equation is now (show this) 

{div a (Ta*Rl*)l* + PaR\*F = 0 
When the applied mass or body forces are zero this reduces to 
div B (T a *Ri*) = 0 

When lv e are using orthogonal curvilinear coordinates it is convenient 
to choose the fixed Cartesian reference frame to be that determined by 
the orthogonal curvilinear coordinates at the initial point (a, 6, c) 
Then Ri is the rotation matrix that rotates this reference frame into 
the reference frame determined by the orthogonal curvilinear coordi- 
nates at the final point (x, y, z) In performing the operation div c on 
the 3X3 matnx T a *R 1 *, we must use the form of the divergence 
operator which is appropriate to the orthogonal curvilinear coordinates 
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we arc using (sec Inlrodvdion to Applied Mathematics, pp. 11I-112 1 ). 
When (x, y, z) nre linear functions of (a, b, c ) the elements of J — 

arc constants: hence the elements of = ./*./ and of r; = 

(o, b, c ) 

X(Af — A’ 3 ) are constants. Conversely, if the elements of 77 are con- 
stants (i.e., if the. strain matrix is a constant matrix), the variables 
(x, j/, s) are linear functions of the variables (a, 6, c). In fact, the 
elements of M are constant and so 2x„ 2 and 2x a X6 nre constants 
(2 denoting summation with respect to x, y, and z ). On differentiat- 
ing the first of these with respect to a, b, and c, we obtain 2x a x aa = 
0, 2x n x„t, = 0, 2x a x a( . = 0. On differentiating the second with 
respect to a and using the relation 2 x a x cb - 0, we obtain 2x?,x, 4a = 0. 
Similarly, 2x f x aa ~ 0, and the three relations 2 x a x aa = 0, 2xtx aa = 
0, 2x c x aa «= 0 assure us, since the matrix ./* is non-singular, that 
x ca - 0, y oa - 0, z an ~ 0. Similarly, xi. b = 0, y bb = 0, z bl , - 0 and 

x rr -- 0, pre ~ 0, z ce — 0. We have already seen that 2x„x a i = 0, 

and, on interchanging a and b, this yields 2x^„j, = 0. On differentiat- 
ing the constant 2x { x„ with respect to b, we obtain 2xeX a f> = — 2x a x tc 
and, on differentiating the constant 2 xtx c with respect to a, we obtain 

2x a Xb with 
2x c x af , = 0. 

The three relations 2x a r af) -• 0, 2x/,x„t = 0, 2 x c x nh — 0 assure us, 
since the matrix J* is non-singular, that x nb - 0, y ab — 0, z a t, - 0. 

Similarly, x ir = 0, y hr = 0, z tr - 0, and x f „ = 0, y e(l - 0, z fo - 0. 

Hence all the second derivatives of (x, y, z) with respect to (a, b, c) are 
zero and so the variables (x, y, z) arc linear functions of (a, b, c). We 
say that a strain is homogeneous when the strain matrix jj is constant. 
We have, then, the following result: 

A strain is homogeneous, when, and only when, the Jacobian matrix 

j « Sl-L'a) „ constant matrix. When this is so the variables (x, u, z) 
(o, b, c ) J 

arc linear functions of the variables (a, b, c). 




2 xiXra. On differentiating the constant 
v_ j v. 


respect to c, we obtain 2x a Xf, c + 2 xizr ea — 0, and so 


When the strain is homogeneous 




J * is a constant matrix and 


fO we have the following result : 

When the strain is homogeneous (he applied force per unit mass is 
zero (so (hat when the applied force per unit mass is not zero the strain 
cannot I'e homogeneous). 

The following sections discuss two important examples of homo- 
geneous htrain. 


s By Frw.nv J). Mumsrttxn, John tVjhy A fvm*, IB? 1 :. 
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2, Simple shear 

In simple shear the deformation is described by the equations 
x = a + kb, y = b, z = c, 


thus each particle is displaced, parallel to the x a\is, an amount 
proportional to its distance from the 6-plane We have seen (E\am- 
ple 1, p 33) that the principal axes of strain are obtained by rotating 


the (a, 6) -axes around the c-axis through an angle a 


K , 1 . _ 

■j + jtan 


i » 
2 


The characteristic numbers of M are (I -f- k am 2a -f- k" 1 mo* a, 
1 — k sin 2a + k 8 cos* a, 1), and if we use the relation k — —2 cot 2 a 
these turn out to be (tan* a, cot 2 a, 1) Hence the characteristic 
numbers of are (tan a, cot a, l), and so 


( cos a — sin a 0\ / tan a 0 0\ / cos a sin a 0\ 

sin a cos a Oil 0 cot a 0 1 1 — sin a cos a 0 ] 

0 0 1/ \ 0 0 l/\ 0 0 1/ 


( sir 
— cos 

0 


sin 2a 
cos 2a 


— cos 2a 

2 cosec 2 a — sin 2 a 
0 


3 


■ a + *■’)- 


/2 k 0 \ 

r !i (fc 2 + k ‘ 0 

\0 0 (4 + fcW 


Hence = (4 + k ’)■ 


/2 + i. 2 -k 0 \ 

-» ( -k 2 0 I 

\ 0 0 (4 + ItW 


and so the 


rotation matrix R = JM M 13 given by the formula 


( 2k 0 \ /sin 2a — cos 2a 0\ 

-12 0 } = I cos 2a sin 2a 0 

0 0 (4 + Jfc 1 )’ 5 / \ 0 0 1/ 


Hence if is a rotation through ^ — 2a = —tan 1 - around the c axis 

Tlius the coordinates of the stress tensor T, in the reference frame 
obtained bj rotating the original (x, y, z )- axes around the z-a\is 

through the angle — tan“ l arc the elements of the matnx 
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The coordinates of T in the original (x, y, r)-referenee frame are the 
elements of the matrix 



Since del ./ - 1 . p- = p a (so that no compression is involved or, 

equivalently, volumes are preserved in a simple shear). The modified 

. . m . d<?> 
stress matrix is 7 a = ./ — — 

dr, 

We first consider the situation in which the medium is isotropic, 
and we take the initial stress to he a hydrostatic pressure p. Then 

4 —71/1 4- — — 2p/ ; -r — ~r /i" — 2 ml il = 4* n/3, 

u 

— — —plJz "h (X/ 1 Ez 4* Spry) 4* (II i" — 2 mIn)Ez 4~ 2 ml\T) 4~ n co 77 

drj 

where 


/() t: 0 \ /0 0 0 \ 

77 « 4 ( I: Jr 0 ), co 77 = 0 0 0 )• 

\0 0 0/ \0 0 -\Jr/ 

Since we have stopped the development of <£(17) at terms of the third 
order in the elements of 77, we are entitled only to keep terms not 

#<f> 

involving J: to a higher power than the second in the expression J — 

dt} 

Since / 1 - ik~, 7 : - — $fr. we have, to this order of approximation, 


dr, 



1 


(X 4' in) Jr 


E 3 + 2707 -f 77 co 77, 


r a c. j 


d 6 

d*l 


4* 


-plh 

f^(\ 4 2 p -r w)/- 2 (71 — p)I; 


& 

0 


0 

4 (X 4* 2;i 4- rn)Jr 0 

0 4(X 4- m - hiU'j 


The force on a matrix element of area r!S~ of the boundary of the 

/ 1 0 ON 

deformed medium is given by T a tlS'\ Since (J*) -1 - j — /.• 1 () 

\ 0 0 1 , 

v.e have ■/.'>' rj rti'i'*, (/y* ~ / 7 .v ** — /„• </X’, r/,v r-r tis f , An element of 
an a which lay initially in the n-planc (».<*., for which cAS 5 ~ 0, dS r - 0) 
deformed into an element of area 0) who-v magnitude is 
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the product of the magnitude of dS a by (1 -f fc 2 ) w The force, per 
unit initial area, on an element of area which lay initially in, or parallel 
to, the a plane is found, then, by multiplying the first column of T a by 
(1 -j- JL*)~ W = 1 — £A, 2 Thus the coordinates of the force (per unit 
initial area) that acts upon an element of area which was initially m, 
or parallel to, the o-plane are (— p + i(X + 2 /j + m -f p)fc 2 , nk, 0) 
Similarly, the coordinates of the force (per unit initial area) that acts 
upon an element of area which was originally in, or parallel to, the 
b-plane are ((#» — p)k, —p + -£(X + 2ft + m)A. 2 , 0), and the coordi 
nates of the force (per unit initial area) that acts upon an element of 
area w hich w as originally in or parallel to, the c-planc are (0, 0, — p-f 
^(X + m — £n)k 2 ) If p = 0 the first-order part of this force per 
unit area is furnished by the matrix 

/0 m 0\ 

A ( M 0 0] 

\0 0 0 / 

This is the result of the linear or infinitesimal theory of elasticity The 
second-order correction that must be added to this is furnished by the 
matrix 

( X + 2/* + m 0 0 

0 X + 2 m + m 0 

0 0 X + m - £n 

/0 0 0\ 

= i(X •f' 2fi -f- ni)k 2 Ez — (m ~h ^n)k 2 ( 0 0 ® ) 

\o o 1/ 

Thus in order to produce the simple shear a hydrostatic tension of 
magnitude -£(X + 2ft + m)k 2 must be applied to the medium and, W 
addition, a simple linear compression ( *= negativ e tension) of magni 
tude (ft 4* ^ n )k 2 parallel to the z-axis, these forces being added to the 
shearing force furnished by the linear theory If these surface forces 
which are necessary, m addition to the shear forces furnished by the 
linear theory, to maintain the simple shear are not supplied, themediuffl 
will contract (because of the absence of the hydrostatic tension) , fur 
thermore, it will expand m the direction of the z axis and contract in 
any direction perpendicular to this, because of the absence of thesimple 
linear compression in the direction of the z axis (see the following sec- 
tion, in which simple tension is treated) Since both the hydrostatic 
pressure ( = absence of hydrostatic tension) and the simplelinear tension 
(= absence of simple linear compression) are second-order effects, the 
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magnitude of each having Jfc 2 ns n factor, we may use the linear theory 
in appraising the effects due to them. Thus the deformation due to the 
hydrostatic pressure is furnished by tlie formulas x' — (I ~ $)x, y' 

- (1 - $)y. t* * {I - fl)r; /s * £ 2 ' It will be shown in 

tlie next section that the deformation due to the simple linear tension is 


n 4- , 

“TT 1 


J* , 


2/ 


1 


t 1 + "A 11 T"l , 
a — — /r j y\ 


„// 




where cr 


. „ . , , T , p( 3X + 2/i) . 

is Poisson s ratio and L — — — — ; is Young s 


2(X 4~ y) X + /i 

modulus. On rombining these deformations with the simple shear a* = 
a 4- kb. y >=• h, r « r, we obtain (on neglecting powers of /; higher than 
the second) 

x" - (1 - «/;*)o + /.?>, !/" - (1 - «fc 2 )b, z" = (1 + yk 2 )c 

4* in 

2(3X 4 - 1/0 ’ 

In order to verify the correctness of this solution we have merely to 
start with the deformation 


X -r 2y 4* tri <r / 1 \ y -f- -}-n X + 2y 

where a <&■ 1 1 u -] — n).y — - — — 

2(3\ + 2u) ^ EV 4 / 7 72 2(3X4 


x *» (1 — a/:‘)o 4" kb, 


2/ 


(1 - «/r)6, s « (I + 7 I- ! )c 

where o and y are undetermined constants. We obtain (always 
neglecting all powers of /; higher than the second) 

/l - air k 0 

./ - ( 0 1 - air 0 

\ 0 0 14- yfrj 

h ~ 2ak- k o 

if ( h 1 4- (1 - 2 a)/;" o 

\ 0 0 14- 2>/;-y 

M. 0 \ 

~ *)*•* 0 L h (|- - 2a 4- v)/r f / 2 «= -|/r, 

->/;V 




(I 


0 


<’4 


| x (h “ 2a + ■») + 5 m 


IrK. 


-2a)r k o 

k (1 - 2a) fr 0 
0 0 2 y/r 


fO 0 0\ 


- nlr I 0 0 0 

VO 0 I, 
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as ~ 2a + 7) + + /i(l — 2a)}& 2 /ifc 0 \ 

J T~ ~ ( A |X(£ - 2a + 7) + + A»(l - 2a)}fc 2 0 ) 

av \ 0 0 {X(£ - 2a + 7) + im + 2 my - *n}JL 2 / 

3^ A* 0\ 

In order that J ~ should reduce to ( nk 0 0 ], the undetermined 

a ” \0 0 0 / 

constants a and y must satisfy the two equations 


2(X + f»)a — X7 = -^X + + m 

2Xa — (X + 2p)y = -^X + ^tn — -J-rc, 


and a simple calculation (perform it) shows that a and y have the 
values given above 

The values of a and y given above correspond to an applied force per 
unit area on the planes, into which planes parallel to the a-plane are 
deformed, of amount yk parallel to the y axis and an applied force per 
unit area on the planes, into which planes parallel to the 6-plane are 
deformed, of amount yk parallel to the z-axis The second of these 
applied forces is a shearing force since the orientation of any plane 
parallel to the 6-plane is not changed by the deformation (any such 
plane being deformed into a plane parallel to it), but the first is not 
Any plane element of area of the a plane is deformed into an element of 
area which is perpendicular to the vector v(l + (7 — «)fc 2 , —k, 0), 
for, on setting a = constant in the equations defining the deformation, 
we obtain dx = kdb, dy - ( 1 — ak^db, dz = (1 + yk*)dc, and so 
dS z dS” dS 1 =» 1 + (7 — a)k* —k 0 (powers of k above the second 
being methodically neglected) If we wish, then, the forces on those 
planes, into which planes parallel to the a-plane are deformed, to be 
shearing forces, the first column vector of T a must be perpendicular to 
the vector t>(l + (7 — a)k 2 , —k, 0) We start out with the trial 
deformation 

x =* (1 — afc a )a + kb, y = (1 — 0fc 2 )6, 2 = (I + yk?)c 


where or, 0, and 7 arc undetermined constants We find, as before, 


M 


( 1 - eck 2 k 0 \ 

0 1 - 0/L 2 0 ), 

0 0 1 + yk 9 / 

( 1 - 2 afc 2 k 0 \ 

k 1 + (I - 20) fc 2 0 

0 0 1 + 2yk 2 / 
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1 05 



(—air 

la 

0 \ 


V « 

ifc 

(L ~ fl* 2 

o ,/. 

= (L — a — /? + y)/r, / 2 = 


\ 0 

0 y/;7 


f v, 

" l x G 

\ 1 1 

-a-/J-Fy)+-Wj 

FEn + 





(-2 air k 0 \ 1 /° 0 °\ 

fi ( 1; (1-2 ft)lr 0 j — - nk" ( 0 0 0 , 

\ 0 0 2 ylr) 4 \0 0 1 / 

/{X(i — a — ft + y) "f L>n + p(l — 2a)}k~ yk 0 \ 

I fik {A(L — a — ft + y) 4* Lw -f- p(l — 2/3) }A“ 0 J. 

\ 0 0 !A(L -«-/? + y) + bn + 2yy \lr - J-,)/ 


The throe equations that serve to determine «, ft, and y are 

(X + 2a)« + A/9 - Ay = La + Lm. 

A« + (A -j- 2^)/3 — Ay = LA + Lm + a*, 

A« -f A/3 — (A -f- 2y)y = La + Lm - 

From the. first two of those \vc obtain ft - a + L, and, on substituting 
thi^ value in the second and third, we obtain 


2 (A + y)a — Ay = Lm, 

2A« — (A + 2n)y — Lm — -J-n, 

and so 

m ^ »A __ m 7i (A 4- n) 

2(3 X + 2p ) '* Sg(3X+lhri’ 7 ~ 2(3A + 2p) + 4p(3X + 2,7)' 

With these values of «, /J, y, T„ — J — is 

M* 8 a/.’ 0\ 

3’« « pfr 0 0 
\ 0 0 0/ 


and the shearing force per unit area on planes into which planes 
parallel to the n-planc and to the fc-plane are deformed is The 
d'fotmatitm is 


nk*)n 4- kb, y «{!-(„ + l)!r]b, 


(1 + y Ir)c, 


and the compression ratio — is furnished bv the formula 


tdi-t jy 




dm — n j 
•i(3A -f- 27} j k 
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The stress matrix T is 




»k 0 \ 
0 0 
0 0 / 


If the initial hydrostatic pressure p ^ 0, the values of the elastic con 
stants (X, n, l, m, ») are different since these are functions of p The 
equations determining a, 0, and y are also slightly modified since 

— and, hence, T a — J~ involve p We have (show this) 
dy drj 


-p + {X(£ — a — 0 + y) + (j i — p)fc 

-f* /»(! ” 2a) + pa}k 2 

pk —p + {A(i — a — 0 + y) -f tyn + p(l — 2/3) 0 
+ PI 3)fc 2 

0 0 — p 4- {X(i - a - /3 + 7 ) + fan + 2^7 i 

“ - P7}fc 2 


An element of area of magnitude dS a of the a plane is transformed into 
the matrix element of area 



whose magnitude is dS a {1 + (i + 7 — P)k 2 } Hence the resolved 
part of the force per unit area on dS*, normal to it, is (show this) the 
product of 


~P + {X(£ - a - 0 + y) + fon - 2ap + p(ct + 0 - 7 ))*:* 

by I - (i + 7 - P)k 2 > le , -p + {X(-£ - a ~ 0 + 7) + - 2ati + 

p(i + a) i k 2 , and if this is to be —p (so that the additional applied 
force on the matrix element of area dS z is a shearing force) the unde- 
termined constants a, 0, 7 must satisfy the relation 


(X "t" 2/x p)a + \0 — X 7 = -J-X + + \p 

An element of area of magnitude dS a of the b plane is transformed into 
the matrix element of area 

/ 0 

dS * = dS a I 1 (7 - a) 

\ 0 

and so the resolved part of the force per unit area on dS x , normal to it, 
is (show this) 

— V + |X(i — a — 0 + 7 ) + + m(1 — 2/3) + p(0 + 7 — a) \k* 
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Hcnro. in order that the additional applied force on dS~ be a shearing 
force, the undetermined constants «, 0 t y must satisfy the relation 

(X 4- p)« + (H 2;i — p)0 — (X 4- p)y + hn 4- a- 

Similarly (show tliis), on considering an element of area of the c-plane, 
we obtain the relation 

(X -f- p)rt -f (X 4“ p)0 — (X 4* 2/i — p)y ~ -A-X 4* -A-m — -J-n. 

When a, 0, y are assigned the values that satisfy these equations, we 
have 

f—p -f (a - lp)k' (ft ~ p)k 0 \ 

T„ -■ ( id: - p 4- p(a - y)lr 0 J, 

\ o o -p + r(n + my 


and, since ~ 1 — (y — a — ff)fr. 

Pn 

f-p-r 1 2#t 4- p(y -0- |) 1 /r (p - p)h 0\ 

7’,,/* = (a - p)* -p 0 . 

\ o 0 ~vj 

Thus the initial pressure affects even the first-order approximation, the 
effect being an apparent reduction of the modulus of rigidity p by p. 
.Since, however, p is itself a function of p all we can say is that p — p 
plays, when p ^ 0, the role played by p 0 when p - 0, p n being the 
value of p when p ~ 0. 

On adding together t he second and third of the three equations that 
icrve to determine a, 0, and 7, we obtain 



(X 4* p)n 4- (X 4- f‘)(0 ~ y) ~ -AX 4 - -Aw 4~ ■A/* — •J-n. 

and, on combining this with the first of the three equations, we obtain 

_ (X 4- m)(p — p) 4 - (Ap — t ?i) (X 4 - 2p - p) — A p(X 4- p) 
3Xj* 4* 2^* p(2X -f- ju) 

Amp 4- Ap(\ 4 -ji) 

3Xa 4- 2p' - P (2X+7T* 

‘hi subtracting the third of our three equations from the second, wc 
obtain 

„ , Aa 4- Jt-n 

o ? 

- p 
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and these formulas determine a, 0, and y, when X, p, m, and « are 
known as functions of p Conversely, a knowledge (obtained hj 
experiment) of how the compression ratio, for example, vanes with 
p furnishes information as to the nature of the dependence on p of the 
function (0 — y) -f- a of X, p, m, n, and p 
The theory is more complicated when the medium that is being 
subjected to the simple shear is, like wood, non-isotropic We con 
sider first a medium that is elastically insensitive to any rotation 
around the a-axis (the direction of the shear) There are, then, three 
additional second-order elastic constants and seven additional third 
order elastic constants (the medium having fifteen elastic constants in 
all, five of the second order and ten of the third order) However, 
since / 1 and many of the elements of i? are infinitesimals of the second 
order and not, as is generally the situation, of the first order, not all 
these ten additional elastic constants will affect our second order 
approximation Adopting the notation of Exercise 12, p 95, where 
the three additional second order elastic constants are denoted by 
(5i, 52, 53) and the seven additional third order clastic constants by 

(5 4 , , S10), the addition to the matrix — due to the non-isotropy 

drj 

of the medium is 


/hvi + 5j(7 1 + jji) 0 
I 0 SsVl 4* & 3 V 3 

\ 0 0 


/— 5i« + 5 2 (£ — 2 a — 0 + 7) 
fc 2 0 

\ 0 


: ). 

S2V1 + 53172/ 

0 0 

— 5 2 a + 5 j7 0 

0 — 5 2 a + — 0), 


(Remember that, when evaluating — > <f> must be written symmetrically 
dv 

and so 172173 — T74 2 , for example roust be replaced by VbbVcc — v belch) 
Thus the seven additional third-order elastic constants do not affect 
our second order approximation (only the three additional second- 
order constants (5i, 5 2 , S3) being involved) The addition to T a is 

the same as the addition to (since the addition to — is of the second 
017 drj 

order m k), and the addition to the stress matrix T is also, for the same 
reason, the same as the addition to The equations that serve to 

di? 

determine the constants a, 0, and y are (the initial stress being taken 



SIMPLE SHEAR OF A NON-ISOTROPIC MEDIUM 109 


to be zero) 

(X -f 2 ft -i- Oj 4- 2o : )« + (X 4' 548 — (X + S;)7 = IX + 4 m + 

{X -f 5j)o 4- (X 4- 2/i)3 — (X 4* o z )y - IX 4- bn 4- n. 

(X 4- 5 3 )« 4- fX 4- 5s )0 - (X 4- 2 ^)y = IX 4- lm - Jn 4- Ifo- 

When the values of «. £. and 7 which satisfy these equations are used, 

fuk- »k 0\ /2pjr 0\ 

Tn « jifc 0 0 J, T = ( a/: 0 0 , 

\ 0 0 0/ \ 0 00 / 

and 

- = 1 -f (« + ^ ~ 7 )/>'. 

Pa 

When the medium is insensitive to any rotation around the 6-axis, 
the additional terms in <t> arc of the form 


4' 5»/i^ 2 4- — ijr.') 4- 4®t , '2 3 4* ^bI\Vi~ 

4" 5 C *) 2 ()J3171 — t;i : ) 4- 5 ; /i(??3Jn — I?3*) 4- 5ftTj«(77c" 4~ ipf) 

4- 5 3 fi(ijf,- 4- nF) 4* 6 io{»7<»)c(»)3 ~ ’h) 4* ur,(nc' — ^-T)!, 


and so the additional terms in — are 


1 * 2 

*010*76 


4* 05TJ3 4* 0**7&~ 0 

0 01*72 + o*(I j + *72) 4* 0s*76~ 0 

» t 2 
•r 0s*?r> 

vj^iO^G*’ 0 02*72 4 0 R*? l + 0 r, 17 C- 

J , 


/ 5;(F — d) 4" 5 j 7 4- 4*9 0 i'5io 


/; ; l 


0 


0 


5i(4 — d) 4“ 5j(l — « — 28 4- 7 ) 

4- 45* 4" 4°’ 

45io 0 o 2 (4 ~ d) — 5;rt 4* 45 9 y 


The additional third-order elastic constants oh, 5 9i and 5m nre now 
involved in our second-order approximation. Furthermore, there is 
now a shearing force on planes into which planes parallel to the c-plane 
nre deformed (this shearing force being 4‘5 io6 2 per unit area parallel 
to the x-nxis). 

If the medium i« elastically insensitive to any rotation around the 
i ^ 

r-s\F, the additional terms in - 7 — (over and above the terms that would 

Or 

t«- pne-nt if the medium were isotropic), are furnished by the matrix 

ft- hr '-a thi-A 
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( &2V3 + 3}JJ2 

— S 2 JJ3 + 53»?1 

0 0 


(hy + Wi — (3) 
+ *’( o° 


0 \ /Oh o\ 

0 =-l Jj 0 0 

hV3 + h(I i + V 3 )/ \0 0 0 / 

0 0 ^ 

h7 ~ $3« 0 

0 hi + 5 2 (£ - = - 0 + 2 y)J 


Thus the presence of the additional second-order elastic constant 5* 
affects the first-order or linear approximation, the effect being an appar- 
ent reduction of the modulus oj ngidity n by S 3 We leave to the reader 
the task of writing out the equations that serve to determine the 
constants a, P, and y m this, and the preceding, case and also the task 
of developing the appropriate theory when the medium is under an 
initial hydrostatic pressure 


3. Simple tension; the linear theory 

The problem of simple tension is that of the deformation of a homo- 
geneous cylinder under the application of equal and opposite forces 
applied uniformlj over the two ends of the cylinder We take the 
c-a\is parallel to the generators of the cylinder and consider first a 
cylinder that is elastically isotropic and initially unstressed We 
start off with the trial deformation 

x = (1 — ok)a, y ■= (1 — ok)b, z = (1 + &)c 

where k is an undetermined constant which is so small that we neglect 
any power of it higher than the first o is a second undetermined 
constant that is, as we shall see a function of the elastic constants of 
the medium (being independent of the applied force) Our trial 
deformation is such that J, and hence y, are diagonal, the three diag- 
onal elements of y are 

ill = — ok, 112 — —ok, y 3 = k 

and so 7 X = (1 - 2 a)k ^ = \I 1 E 3 + 2 ny is diagonal its first and 
oy 

second diagonal elements having the common value |A(1 — 2a) — 
2/iaj k Since we are neglecting powers of k higher than the first, the 

ptft dd> , 

matrix T a = J — w the same as the matrix — > and, since there is no 
ay dy 

applied force on the sides of the cjlinder, o must be such that the first 
(and second) diagonal elements of T a are zero (for the matrix clement 
of area has its third coordinate zero on the sides of the cylinder and 
thus at least one of its first two coordinates is different from zero) 
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lienee 


X 

2 (\ + »y 


c i« known ns Poisson's ratio; for many media its value is about 0.3. 
If l is the (initial) length of the cylinder the two ends of the cylinder 
may h<‘ taken to ha given, before the deformation takes place, by the 
equations c ~ 0, c ~ I. If the magnitude of the applied force per 
unit (initial) area on the deformed position of the end c — l is denoted 
by P, we have (since the third diagonal element of T a is \1\ -f- 2gi?3 — 
SX(1 - 2 -t) -f 2, i)/.- = 2g(l + o)k), 2g(l + o)k = P and so k = 
p 

\y e denote 2g(l 4* c) bv E and we term E the Young's 

2p(l 4- tr) 

modulus of the medium: 


E 


Ml 4- <r) — 


n(3 X 4~ 2g) 
X 4- p 


Then k - Thus our assumption that k is so small that its square 

may he neglected means that P is so small in comparison with E that 
p 

the square of ~ may be neglected. Since dx dy — (1 — ok) 2 da db — 

! 1 — 2 ok) da db, any element of area of a cross section of the cylinder is 
decreased in magnitude by the factor 1 — 2 ok; if, then, we denote by 
P' the magnitude of the stress on the end of the cylinder (i.e., the 
magnitude of the applied force per unit, final area), we have P'(l — 
2<?k) — P, and so 

P P % 1 - 2 ok) 

: " E ' " E 


On multiplying both sides of this equation by 1 4- 2 ok and neglecting 
k\ we obtain 


1 : 


r 

E 


In vords, to the degree of approximation furnished by the linear theory 
it makes no difference whether P denotes the applied force per unit 
instnd area or the applied force per unit final area. In either event 
p 

k r * *v and the defonnation L furnished by the formulas 


JT * 





r. 
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The relatne extension ^ C - is equal to ^ The graph of P, as a 

function of the relative extension, is a straight line through the origin 
having slope E 

If the initial stress, instead of being zero, is a hydrostatic pressure p, 

the first diagonal element of — is —p -f- |X(1 — 2a) — 2(ur}k and so 
dn 

d<j> 

the first diagonal element of T„ — J — is — p + {Ml — 2<r) — 2 ya *f 

dt\ 

pa\k An element of area which was originally normal to the a axis 
remains normal to the a-axis after the deformation but its magnitude is 
increased by the factor 1 + Ml — tr) The applied force on such an 
element being a pressure p, wc must have — p + { X(1 — 2 a) — 2ya -}- 
perji. = — p{l + Ml — *)} andsoX(l — 2<r) — 2 n<r + p == 0 Hence 


a 


X + p 

2(X + a) 


Thus the effect of the initial pressure is to change the formula furnish- 
ing Poisson’s ratio to the extent of replacing the numerator of the 
fraction by X -f V Of course, X and y are themselves functions of p 
so that a depends on p not only explicitly (through the term p in the 
numerator X + p of the fraction) but also implicitly through X and y 

The third diagonal element of zr ts ~P + 1*0 — 2<r) + 

dr] 

and so the third diagonal element of T a is — p + {X(l — 2<r) + 2y — 
p\k An element of area which was originally normal to the c-axis 
remains normal to the c-a\is but its magnitude is decreased by the 
factor 1 — 2 ak If, then, the applied force is of magnitude P per unit 
initial area (in addition to a pressure p), we must have 


-p + {X(l - 2<r) + 2y - pi* = P - p( 1 - 2 ak), 

and so 

\X(1 -2a) +2 y-p - 2pa\k = P 
or, equivalently , 
p 

k — — where E — X(1 — 2a) + 2p — p — 2 pa 
L 


As w hen p — 0, it does not matter w hether P denotes the magnitude of 
the applied force per unit initial area or the magnitude of the applied 
force per unit final area (since we are neglecting <l 3 ) Since X(1 — 
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2?) «= 2pcr -- p we have 


E - 2(/i - p )( I 4 a) 


(n — p)(3\ 4 2/i 4 p) 
X 4 M 


Thus the effect of the initial pressure, in the formulas furnishing <r 
and E, is to increase X by p and to decrease /i by p (so that the denom- 
inators of the fractions that furnish a and E are unaffected formally; 
they are, of course, affected by the fact that X 4 u is a function of p). , 
If the initial stress is a non-scalar stress, we know that the medium 
cannot be elastically isotropic (why?). Let us consider an initial 
stre-'s matrix all of whose elements arc zero save the one in the third 
row and third column, and let us denote this non-zero element by P. 
We assume, further, that the medium is elastically insensitive to any 
rotation around the c-axis. Then <j(t?) is of the form 


<^(tj) ~ /hjs -j — ~ 7i* — 2 /i/; - 5)7)3* 4 0;/ iris -f 03(7)1772 — 7j r> ') 


and so, since tj is, bv hvpotbesis, a diagonal matrix, — is a diagonal 

dr) 


matrix whose diagonal elements are 


XI 1 -f 2/071 4 027)3 -f O3772, X/ 1 4 2/17); -f- 0;7)3 -f- 637)1, P 4 X/ 1 

4 2 / 17)3 -f 0 1 f) 3 4 5 ;(/ 1 4 7)3) , 


i.e., 


!X(l — 2 <t) — 2/kt 4 5; — ojir j/;, {X(l — 2<r) — 2/nr -f- 6; — ojer } A*, 

P + |X(1 - 2<r) + 2 /j 4 6j 4 25;(1 - 


Hence, since the first and second diagonal elements of T a 
he “(to (why?), a is determined by the formula 


- J 


• — must 
dr) 


X -j- 6; 

2\ 4 2/i 4 63 

The equation determining /: is 


{XU - 2c) 4 2/i 4 h 4 2o;( I - <r) 4 P\k - A P 

v hen' /* 4 aP j>; the applied force per unit initial area. The coefficient 
of k b Young’s modulus E and, since X(1 — ‘2c) ~ — <5- 4 6 ;{ r, 

we have 

E - 2ix(! 4 <*) 4 h 4 5;(1 - BO 4 5y» 4 P. 
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Since k is the relati\e extension — — — > we have 


z — c _ A P 
c ~ E 

where E is a function of P both explicitly and through the elastic con- 
stants n, a, i), Sgj nnd S 3 On denoting z — c hy Ac, wc obtain, on 

letting A P — * 0 , from the approximate equation — = — the exact 

c E 

equation 

dc _ dP 

7 ~~E' 


and, on integrating this, 
lo{ 




where l is the length of the cylinder when there is no applied force and 
A l is the extension when the applied force per unit final area is P' 
Note P' is the applied force per unit final area, and not the applied 
force per unit initial area, since P is the applied force per unit area, 
the area being evaluated w hen the c> Under is under the tension P 


In order to perform the integration 


r r dp 
Jo E' 


wc must Know how the 


Young's modulus E of the medium vanes with P If w , o make the 
drastic assumption that n and a are independent of P and neglect the 
lack of isotropy introduced by the application of the tension /’, E 
is the linear function 2/<(l + e) + P of P and wc obtain 


lo g (l+y)-log(l +0 

or, equi\ alcntly, 

At _ P^_ 
l E 0 


where Eo = 2/i(l 4- <r) is the Young's modulus of the medium when 
unstressed Thus, under these unrealistic hypotheses, the extended 
linear theory furnishes the same result as the simple linear theory, 
namely 

The relatne extension is proportional to the applied stress (or force 
per unit final are*'' 

This result ike’s law m its simplest form The formula 1 *• 

< 
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j* mav be regarded as the general form of statement of 

\pj <>n 

Hooke's law. 

If we suppose that E may be sufficiently approximated, over the 
interval (0, P’), by the linear function E 0 -f aP, where a is a constant 
Hint is determined by experiment, we obtain 

l0B log (' + a k) 

or, equivalently, 



If n happens to be negative this formula indicates that Af— + as 

p 

P’ ----- In fact, on setting a - -/9, we obtain 
— a 


A! 

/ 


0 


1 


- 0 



- 1 , 


and, as P’ —* — • ( 1 — ft — V — * 0. Thus the extended linear theory 

is capable of explaining the phenomenon of flow of a cylinder under 
tendon; but in order for it to do so one must take into account the 
dependence of the elastic constants on the applied tension. It seems 
hardly likely that the mere dependence of ;i and o on P would be 
sufficient to account for this phenomenon of flow; it is to be expected 
that the elastic constants 5;, 65 (which measure, by their differences 
from rero, the lack of isotropy of the medium) must be considered. 


4. Simple tension; the sccond-ordcr approximation 

We consider first a medium, supposed to be initially unstressed, that is 
hntropie. We assume the trial deformation x — (1 — ok -f- filr'ja, 

P 

V (l — ek -f / 3k‘)h. : ~ (1 + /; 4 o!r)c where /.* *= — and where a, 

J j 

E are the elastic constants (Poisson’s ratio and Young’s modulus) of 
the «lmp!e linear theory; 3 and 0 arc two new undetermined constants. 
V, .• neglect methodically powers of /; above the second, v, is a diagonal 
matrix- xxho-e diagonal elements are 

— ok 4 (5 4 = t;; 
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and so 

lx = (1 - 2 a)k + (20 + 5 + \ + o 2 )k 2 

J s = (o 2 - 2ir)A: 4 

The constant a is such that the coefficient of k m the first diagonal 
element of X/jEj 4- 2/xij is zero and so the first diagonal element of 

7^ is the product of A. 2 by X(2/3 + 5 + H ffJ ) + 2/i(0 4~ -Lr 2 ) -f 
dr) 

1(1 — 2«r) 2 + 2rrur(l a) — no Since we are neglecting powers of 
k above the second, this expression also furnishes the first diagonal 
element of T a (why?) and so the two undetermined constants 0 and 5 
must satisfy the following relation (why?) 

2(X + #i)/3 + xa = -X(i + v 2 ) - HIT 2 - 1(1 - 2cr) 2 

— 2mo(l 4* it) 4* tiff 

The constant E is such that the linear part of the third diagonal ele- 
ment of T a is P, the applied force per unit initial area Since the third 
dd> 

diagonal element of — is {X(l — 2a) 4- 2^]fc 4- {X(2/9 4- 5 4- £ + 

v 2 ) 4- m(2S 4- 1) 4- 1(1 - 2<r) 2 4* 2m(l - o 2 ) 4- rur 2 \k 2 , the third 

diagonal element of T a = J ~ is P 4* {X(2/S -f S 4- $ — 2<r -4 v 2 ) 4* 
o*i 

ji(25 4- 3) 4- 1(1 — 2<r) 2 -f 2m(l — o 2 ) 4- rur 2 \k 2 On setting the 
coefficient of k 2 in this expression equal to zero, we see that the unde- 
termined constants 0 and 5 must satisfy the relation 

2X/S 4- (X 4- 2ji)5 = -X(£ - 2(r 4- it 2 ) - 3h - /(I - 2<r) 2 

— 2m(l — a 2 ) — no 2 , 

and, on solving for 0 and 5 the two linear equations that they satisfy, 
we find 

X(3X 4- 4 M ) X(3X 2 4- GX, U 4 V) 

P 8(X + I*r (3X + 2 U )(X + + 8K3X + 2rf(X + " 

, _ _ 3 (3X + 2ri t 

2 (3X + 2rf(X + ri ! 2(X + ii)‘ m 

3X 1 

M3X + 2rf(X + M ) 

(z — c) P 

The graph of the relative extension — - — plotted against k = 

is a parabola passing through the origin with slope 1 and second 
dem stive 25 
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The dependence of Young’s modulus E on the applied tension, 
recording to the second-order approximation, is obtained ns follows: 

(z — c ) P 

On writing the relative extension in the form —> we have 

c ii 


IL 

E' 


* k + Mr 



and so 



or, equivalently, 



Thus the constant a of the extended linear theory is —5 if the second- 
order approximation is adopted. In order that the second-order 
approximation may furnish (through the method of the extended linear 
theory) an explanation of the phenomenon of flow, it is necessary that 
c he positive; this being the case, the medium will flow before the 


applied force per unit area attains the value 


E a 


To examine the effect of an initial hydrostatic pressure j> (following 
the theory of the second-order approximation), we observe that the 


iir.-t diagonal element of — is 

On 

T I Ml ~ 2<r) — 2fi<r i I: -f- 1X(2/1 -}- 6 -f- b + c~) + g(2/3 

+ tr) + i(l - 2<r)' J + 2m<r(l + c) - r,~\!r, 

*4 A 

mul so the first diagonal element of 1\ = J — is 

On 
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the coefficients of k 2 we obtain the equation 

2(X + m) 0 + (X + p)S - -X(i - <r + 3c r 2 ) - 3 m* 2 - ?(I _ 2cr) 2 

2m<r(l -j~ o) 4" fltr + pa 

connecting the undetermined constants 0 and 5 On using the relation 
p = 2(X + jx)<y — X we may write the right hand side of this equation 
in the form — X(£ + <r) 2 — pa 2 — 1(1 — 2<r) 2 — 2mcr(l -f cr) -f rw 
To obtain a second equation connecting 0 and 5 we observe that the 
d<J> 

third diagonal element of — is 
6tj 

-p + {\(1 - 2j) + 2 flit + {X(20 + 8 + i + „ ! ) + M2S + 1) 

+ ((1 - 2 *) J + 2m(l - „*) + 

The third diagonal element of T a is the product of this expression 
by 1 + k + Sk 2 i e 

-p + {X(l - 2cr) + 2 m - p\k + (X(20 + « + f - 2<r + <r 2 ) 

+ m( 25 + 3) + 1(1 - 2cr) 2 + 2m(l - <r 2 ) + na 2 - pa)k 2 

and this must be equal to the sum of P (the applied force per unit 
initial area) and the product of — p by the ratio of dxdy to da db 
namely 1 — 2<r k + (2/3 + a 2 )k 2 The Youngs modulus E is deter 
mined (in accordance with the linear theory where the initial pressure 
is different from zero) so that the linear parts of the t\\ o sides of the 
equation obtained in this way are equal Equating the coefficients of 
k 2 we obtain the equation 

2(X + p)P + (X + 2m - p)5 - -X(! - 2«r + «r 2 ) - 3 m 

- 1(1 - 2<r) 2 - 2m (1 - a 2 ) - no 2 - p * 2 

connecting the undetermined constants /3 and 5 The determinant 
of the two linear equations that /? and S must satisfy is 

2(m ~ P)(3X + 2 m + P) 

and so 0 and 5 are determined unambiguously by the two equations 
unless p — m The theory breaks down when p = m because the 
effective rigidity of the medium is m ~ V and so the smaller is m V 
the more liquid like (i e less able to support tension) is the medium 
It is important to remember that ft is a function of p it may nell be 
that m “ P is for positive values of p, an increasing function of p 
and so m — P is never zero 
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1. The compression of a spherical shell 


We consider an elastically isotropic spherical shell whose internal and 
external radii, when the shell is free from stress, are R, and R n respec- 
tively. The problem we wish to solve is the following: Wlmt is the 
deformation of the shell under an applied internal pressure p, and an 
applied external pressure p/l From reasons of symmetry each particle 
of the shell undergoes a radial displacement whose magnitude is a 
function of the initial distance r of the particle from the center of the 
shell. If we denote this displacement by ku, the space polar coordi- 
nates (relative to a reference frame whose origin is at the center of 
the shell) of the final position of a particle whose initial coordinates 
had the space polar coordinates (r. 0, <f>) are (r -f* ku, 0, <}>). Here k is 
a constant which is such that in the linear theory we ncglcet powers of 
k above the first whereas in tho second-order approximation we neglect 
powers of k above the second. The space polar coordinates of the 
vector da are (dr, rdO, r sin 0 d$), and the space polar coordinates of 
the vector dx arc {dr -f /: du, (r + ku) dO, (r -f kit) sin 0 d$). Thus 
da is the 3 X 1 matrix 


and dx is the 3 X 1 mat rix 



(dr -j- k du 
(r -f ku) dO 
\(r -r ku) sin 0 dd/ 


Hence the 3 X 3 matrix ./ that transforms, by means of the formula 
dx >■= J da, da into dx i« the diagonal matrix 


J « 


+ ku' 
0 


no 



0 


0 
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where the prime indicates differentiation with respect to r It follows 
that the space polar coordinates of the strain matrix rj are furnished by 
the elements of the 3X3 diagonal matrix 



Neglecting powers of I above the second we have 

/, = k (v.- + 2 “) + {(„')’ + 2^j, 



Since the medium is by hypothesis isotropic and free from initial 
stress 


t 2 - \IiE a + 2m 4 . (f/, 2 - 2mh)E s + 2ml lV + 7ico v 

drj 

(In this equation, <f> is, of course the energy density and not the third 
d<j> r d<l> 

space polar coordinate ) Thus — is diagonal and so T a = J ~~ is 
3>y or) 

diagonal We assume that the shell is free from applied mass forces, 
and we find the conditions that a must satisfy by equating to zero the 
3X1 matrix div B T a * or, equivalently, since T a * = T a the 3X1 
matrix div„ T a Note The axes of the space polar coordinate refer- 
ence frame that is attached to the final position (r -f* hu , 8, 4>) of any 
particle of the medium have the same directions as the axes of the 
space polar coordinate reference frame that is attached to the initial 
position (r, 8, <j > ) of this particle In other words the rotation matrix 
Ri that transforms the latter of these two reference frames into the 
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former F the identify matrix E s . Tims the equation div„ (T n *R x *) - 
15, which expresses that the mass, or body, force? arc zero, reduces to 
dir„ (Tn*) - 0. In evaluating div c T„ we must use the expression for 
she divergence of a 3 X 3 matrix in space polar coordinates. This 
f'\pro--ion, and its derivation, are given in Introduction to Applied 
Mathematic*, p. 112. 1 For a diagonal 3 X 3 matrix whose diagonal 
dements rr, 00, £?> are functions of r alone, the divergence reduces to 
the 3 X 1 matrix 

rr' -f - - rr — - 00 $<}>, (0$ — o\ 

r r r r / 

and (.iriee the second and thiid diagonal elements of T„ arc equal it 
follow." that it must satisfy the single condition furnished by the equa- 
2_ 1 ~ 

tion rr' + - rr (00 -f- $$) — 0. 

r r 

We first obtain the solution furnished by the linear theory by 
neglecting powers of / t above the first. T a reduces to the diagonal 
matrix whose diagonal elements are 

-f- 2 -j- 2 'fill 1 




and so u must satisfy the differential equation 


* ft + 2a) 


( 


u" + 2 • 



- 0 . 


Use general solution of this equation is w =» Ar + — where A and B arc 

r~ 

rotislants which must be determined by the facts that T a dS* has the 
v; Urn dS' when r — lt { and the value —p c dS x when r = R e . The 

'' N 1 matrix dS z is the product of the 3X1 matrix dS a by ^1 + k 

".Imn-eonJ and third elements of each of these matrices being zero). 

A 

Hence .-1 and B are determined by setting rr — — p t - ^1 -f- /: - J when 
P, and rf — — p f -f f: when r = R r . Thu? p , and p r are 

* -b f r.a.:* D. ' ’umsgbnn. John Wiley A- Son-, 10 IS. 


r 
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of the same order of magnitude as the products of X and n by k and, 
to the degree of approximation afforded by the linear theory, we may 

replace — p, + A, ^ and — + by — p, and — p e , 

respectively When we substitute u = Ar + ^ in the expression for 
rr we obtain 

?? - * |(3X + 2 m )A - ^J, 
and so the two equations that serve to determine A and B are 
(3X + 2„ W -^=-(!) 

(3X + 2„)A- 4 -0=- 

Hence 

. P,R. 3 - P.R , 3 (P. - P,)RW 

fc(3X + 2 I 0(R,’ - R, 1 )’ iMR, 3 - R,") 

We denote by Ui this linear theory solution, thus the radial displace- 
ment fcui that is furnished by the linear theory is 

, P,R, 3 - P.R.’ , (i>, - v.)RW 1 

‘ (3X + 2 ri(R, 3 - R, 3 ) MR 3 - R. 3 ) r* 

The solution furnished by the second order approximation is of the 
form u = Ui kw where, in determining the equation that w must 
satisfy, we may replace u by Ui in the coefficient of k 2 in T a (for we 
neglect, in the second order approximation, powers of k above the 
second) T a is the 3X3 diagonal matrix whose diagonal elements are 
furnished by the formulas 

rr = k jx ^u' -f 2^ -f 2/iu'J + k 2 jx Q u' 2 + 

4* 3 ftu' 2 4* l + 2 + 2m ^u' 2 — 4- n Jjj’ 

00 = k jx ^u' 4- 2 ^ 4* 2(i 4- k 2 jx Q«' 2 4-^ + 3 ^ 

+ 3 ,^f(n' + 2f) 2 + 2m(^-f) + nfj=^ 
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where it is understood that, in the coefficients of Zr\ v is to be replaced 
byvu whereas, in the coefficients of /.*, v is to be replaced by «i -j- lew. 
Oft calculating div„ T„ we find that the coefficient of k is zero (pre- 
rjccly because v\ was determined by means of the linear theory), and 
on equating the coefficient of k~ to zero we find that v: must satisfy the 
equation 


(A + 2m) (i v" + 2 y - 2 = 2(X + 3 P + 2m) 




Hence tr rt Cr 4- — + 
r 


I) (A + 3 m + 2m) B- 


1S(A + 3 m + 2 m)B- 


where C and D are two 


(A + 2m)i' j 

muletermined constants. The radial displacement furnished by the 
-mind-order approximation is 


ku hv\ 4- Irv' » (Avl + k"C)r + 


(l:!i + IrD) 


+ 


(A + 3m + 2m)k'B~ 
(A + 2 M)r l 

'Hie constants .4 and B have been determined, by the linear theory, 
so that the coefficient of k in rr — — \^j' "hen r = / f,-, and — 
when r R ( . We now determine C and D so that the coefficient of k" 
in rr ~2/q when r — and — when r ~ R f . For 

this to he the case we must have 

.:>X 4- 2 m)C - (,\ a- - 7 V) 

Ili* k V /?, 7 

, tax -f 10 m)(X 4- 3m 4" 2m) if 


X 4' 2m 


i) b- , /9 0.1 /f , 3/r ; \ 

It? + R?) 


atai n similar equation obtained bv replacing /?, bv /?<. and p, bv p,. 
*2j 
k 


Hn replacing — by (f>,\ 4- 4 m) . 4 — y-r and collecting terms wc 

/It 
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(3X 


+ 2 P )C - 4m^ s = (|x + m - 91 - nj A 2 


+ (6X + 4m 4* 6m - 


-(' 


m( 7X + 10 M - 8m) 
X + 2jx 


\*! 

/ R t 6 


On replacing R t by R, we obtain a second equation connecting C and 
D, and on solving these two equations for C and D we obtain 

A 2 4 - ( p( 7 * + 1( V ~ 8m ) , n 1 S 2 
+ j(X + 2 m )(3X 4* 2fi) ^ 3X + 2mJ R t *R e 3 ' 


(1 _ 9l + n\ 

\2 3X + 2a/ 

■ - - ^3\ + 2,, + 3m 




AB + 


(7X + lO/i - 8m) 
4(X + 2/i) 


+ 




Special Cases 
1 The solid sphere 


pressure simply by p) A 


Here R, = 0 and so (if we denote the applied 
9 l + n\ 


t B = 0, C 


k\ 3X + 2>i) 2 


A(3X + 2m) 

The radial displacement is 

/I _ 9/ + « \ p 2 1 

' \2 3X + 2m/ (3X + 2m) 2 j ? 


/\ 9 l + n\ 

: \2 3X + 2m/ 


1 3A + 2m 

It is easy to verify that this expression agrees with the result furnished 
by the second order approximation in the treatment of the hydrostatic 
compression of an isotropic medium (carry out this verification) 
(Hint If the strain matrix is — eE s the coefficient of r above — —e — 
•£e 2 since, if we denote this coefficient by 6, (1 4* 5) 3 = (1 — 2e) w ) 

2 The spherical canty m an infinite solid Here pz = 0, Rz — 00 
and so (if w e denote the pressure on the surface of the cavity simply 
by p and the initial radius of the cavity simply by R) A — 0, B = 

The radial 


4*m R ’ C " 
displacement is 


0, D = 

[i + 


7A 4~ 10m ~ Bm 
4(X 4- 2m) ^ 

[ 7X 4- 10m ~ 8m 
1 4(A 4- 2m) 


£)(£)' 


3 The thm spherical shell Denoting the thickness of the shell bj 
e, we have R e = R, 4- e On keeping in the expressions for A and B 

only the lowest powers of we obtain 
re. 



Tin: comprkssiox of a circular tube 


12.3 


1 i (7>» - Pi) Hi M} V") fi, _ ^ 

3(3X -f 2/0 7‘ Uf “ 12^ < ' 


<nv; therefore, Avl — - — — It follow from the expressions for 
3X -f 2 jj 

("and I) that 

. |YI 91 + » \ 1<V , m(7X 4- 50m ~ Sm) 

■” r " i\2 3 A 4- 2 m ) (3X 4- 2m)" + (X 4- 2/.) (3X 4- 2m) 


ji ) „ 

4 <?-. 

3X 4- 2m) 


3X 4- 2m V 


3\ 4- 2m 4- 3 m 


1 \ (7X + I0 i? r r jl) 

2 n j ‘ 2(X 4- 2 m)” 


+ 2 m) r ‘ 

Tito radial displacement is found by .substituting these values in the 
c\pro} tion 

ft. , . h ? r . M 4- IrO (X 4- 3m 4- 2ei)/.-7U) 

\ U T 1 1 + + "*lx + 27jr“ ' ) 
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whose magnitude is a function of the initial distance r of the particle 
from the axis of the tube If we denote this displacement by kit, 
the plane polar coordinates (relative to a reference frame whose origin 
is on the axis of the tube) of the final position of a particle whose initial 
position had the plane polar coordinates (r 0) are (r -f- Lu 6) Here 
A, is a constant which is such that m the linear theory we neglect powers 
of k above the first whereas m the second order approximation we 
neglect powers of k above the second The plane polar coordinates 
of the vector da are (dr, r dO), and the plane polar coordinates of the 
vector dx are (dr -f k du, (r + ku) dd) Thus da is the 2X1 matrix 

(rdfl) whereas dx is the 3X1 matrix ^ (r^ki^do) aD ^ ^ X 2 
matrix J that transforms by means of the formula dx — J da, da into 
dx is the diagonal matrix 

fl+ku 0 \ 

0 i+k v 

where the prune denotes differentiation with respect to r It follows 
that the plane polar coordinates of the (plane) strain matrix i; are 
furnished by the elements of the 2X2 diagonal matrix 



Neglecting powers of k above the second we have 

i, . ^ 

r 

Since the medium is by hjpothesis isotropic and free from initial 
stress, 

= X/jZJj -j- 2/tij -f* (i/i* — 2uiI'>)Ez -f- 2/n/iij 
drj 

Note Since we are dealing with a two-dimensional problem so that 
the strain matrix is two-dimensional rather than three dimensional, 
the invariant / j does not appear and we have only the two third order 
elastic constants l and m rather than the three third order elastic 
constants l, m, and n 
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” find T a = ./ — are diagonal. We assume that our tube is free 
bn dij 

from applied mass forces and so we find tlic conditions that w must 
satisfy by setting div„ T a = 0 (win-?). For a diagonal 2X2 matrix 
vdiO’C diagonal elements rr. 00 are functions of r alone, the divergence 
(in plane polar coordinates) reduces to the 2 X 1 matrix 

^rr' + “ (rr — f6), 0^ 

(sec Introduction to Applied Mathematics , p. 98)' and therefore v must 
pati«fy the single condition furnished bj' the equation rr' + - (rr — 00) 

£*- n. 

We first obtain the linear theory solution by neglecting powers of 1: 
above the first. T a reduces, then, to the 2X2 diagonal matrix whose 
diagonal elements are 


rr 


H 


(X -f* 2 p)u' -f- X — 


Ob = k Jxu' + (X + 2p) "J 

and so w must satisfy the differential equation 


ru" + u’ - - - 0. 
r 


B 


T)n> general solution of this equation is u = Ar -} where A and B 

r 

are constants which must be determined by the facts that T„ (!S a has 
dw value — p, r!S~ when r — R t - and the value —p, dS~ when r — /?«.. 
Uic 2 X 1 matrix rh S' is the product of the 2X1 matrix tlS a by ] 4- 

{ thus A ami li are determined by setting rr =■- — p, ^ 1 

v.hcn r It; and rr - — p, -f I: when r — R f . It follows that 

E nud p,. are of the same order of magnitude as the products of X 
and ^ by k, and so (to the degree of approximation afforded by the 
fitirsr theory) we may equate rr to -p,- when r - R; and to ~p f 


v.h<- 


B 


< n r When we substitute v =* Ar 4- in the expression for 


P. Muroscbnn, John Wpcy A: Son*, 10 IS. 
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— T-t- x) D = 0 The radial displacement is 

4 fc 2 (X + u) 2 F 

-i* (i . L . a + p ) 

2(X + a») 1 V T X + m/2(X + m)J 


2 The thm-walled circular tube Denoting the thickness of the walls 
of the tube by t we have R e = R» + * On keeping in the expres- 
sions for A and B only the lowest powers of ^ » we obtain 


IA 


Pi - Vt R* 
4(X + /i) *’ 


kB = (p, ~ y t ) R, 

R* An * 


say, and therefore kA - 
formulas 


X + M 


P C and D are furnished by the 


« = -{( 2 + !T7)(xT->T^)x^ 


k 2 D - - 


f 3X 6/i -f- 4m 2 (m -j- 2m) 1 


-yr . ■ . z . rt & 

X + m X+2 n) 

The radial displacement of any point on the inner wall of the thin tube 


m-M t 


2nP 

.X + n 


3X + Gfi + 4m ; 
+ X + p + ‘ 


m(m + 2m) 

(X + m)(X + 2m) 


+ Gm — X] 

2(X + 2 u) I 


«• 


3. The torsion of a circular cylinder 

We choose the axis of our c\ Imder as the axi« of a system of cylindrical 
coordinates (r, 6, z), and we denote by r 4- ku, 0 + kz, z + kw the 
final coordinates of the particle of the cylinder whose initial coordinates 
were (r, 0, z) Here fcisa constant parameter which is such that its 
square and higher powers may be neglected in the linear theory whereas 
its cube and higher powers may be neglected in the second-order 
approximation We assume that u is a function of r alone, that u> is a 
function of z alone, and that the medium is isotropic and free from 
initial stress, and we try to so determine the unknown functions u =■ 
u(r), to = w(z) that 
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(1) The deformed cylinder is in equilibrium under the action of 
Mirfnn- force.* 1 alone, the mass or body forces being zero. 

(2; There are no applied forces on the sides (or curved portion of 
the bounding surface) of the cylinder, and the forces applied to each 
end ft f the cylinder reduce to a couple about the axis of the cylinder. 

The cylindrical coordinates of the vector da are (dr, r dO, dz) whereas 
the cylindrical coordinates of the vector dx are {(I + hi') dr, (1 + 

(rdO -j- hr dz), (1 km’) dr} where id denotes the derivative of 
r 

u with respect to r and ir' denotes the derivative of tc with respect to 
z. No!'. The cylindrical coordinate** of dx have reference to a 
frame obtained by rotating the reference frame in which the cylindrical 
coordinates of da are (dr, r dO, dz) through an angle kz around the z-nxis. 
The 3 X 3 matrix that transforms, by means of the formula dx — J da, 
tie- 3 X 1 matrix 

/ dr \ 


into the 3 X 1 matrix 


da -- I r dO 
\ dz 


(1 -{- kid) dr 

(l + k r dO -i- k(r 4- hi) dz 
(1 -f hd) dz 


/I + hd () 0 \ 

J - 0 14 -/.-“ k(r + ku) • 

' 0 0 1 + hr’ / 

Hence the cylindrical coordinates of the strain matrix 17 are given by 
the formula 


id 0 () 

„ u 1 


o r u 


fa')" 0 


(1 r'r 4- r- 


tponer* of /. above tlm second b-ing neglected*. It follows that 
/, . h(>dx. , 1. f -) 

V ■ f * “/*' (2 * ~ 2^ r 2 Ucr + 
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l , 
" 2 ™ 


We first determine the linear thcor> solution Here 



u' H h w' I + 2nu' 




+ 2m - 




+ w‘ ) + 2nw l / 


and T 0 = J — is the same as this Since we are using orthogonal 

OJJ 

curvilinear coordinates the matrix whose divergence must bo equated 
to zero (to express the fact that the mass or body, forces are 2 ero) 
is T a *Ri* iv here R x is the rotation matrix that transforms the cylin- 
drical coordinate reference frame at (r, 0, z) into the cylindrical coor- 
dinate reference frame at (r + ku, 0 + kz, z + kw ) 

( cos kz — sin kz 0\ /0 — z 0\ 

sin kz cos kz 0 j = Ei -f- k I z 0 0 J 4- 

0 0 1 / \0 0 0 / 

Since we are neglecting (when obtaining the linear theory solution) 

powers of k above the first the matrix whose divergence must be set 
equal to zero is 





+ 2 „- 


X (u r + + 2 \xv? } 


Using the formula for the divergence of a 3 X 3 matrix in cylindrical 
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coordinates (pee Introduction lo Applied Mathematics, p. 112), 1 we 
obtain the following equations: 

«" + - - “ » 0, w" = 0. 

r r* 


Hence u is a linear combination of r and - and, since u must be bounded 

r 

at r — 0, it follows that v is a multiple of r. Also tr' is a constant and 
m) we may take, without loss of generality (why?), to to be a multiple 
of z: 

u — Ar, w — Hz. 


The fact that there is no applied force on the sides of the cylinder 
yields *2(X -f n)A + X/i ~ 0. The coordinate in the direction of 
the axis of the cylinder of the force, per unit initial area, on the end 
z >- l of the cylinder is constant and equals 2X/1 + (X + 2 y)B. In 
order that the force on the end z — 1 of the cylinder should reduce to a 
couple we must have, then, 2XA + (X -f 2y)B — 0. The two equa- 
tions 2(X -f* fi)A 4- X/i — 0, 2X/1 -{- (X -f- 2y)B = 0 yield A — 0, 
B =“ 0. TThus the solution furnished by the linear theory is Ui = 0, 
ir i •-> 0; thus any cross section z = constant of the cylinder is merely 
rotated through an angle kz. The force, per unit initial area, applied 
to the end z ~ l of the cylinder is perpendicular to the radius to the 
point of application and is of magnitude kyr. Hence the system of 
forces applied to the end z — / of the cylinder is equivalent to a couple 


or torque around the axis of the cylinder of amount 



If o is the angle through which the end z — 



l of the cylinder 


is rotated we have a ~ kl and thus we obtain the formula 


Torque — 


errg/? 4 

~2 r 


v-hich serves to determine the modulus of rigidity a. 

In order to obtain the solution furnished by the second-order 
approximation, we set, since u i and ir\ are zero, n — h\ and ir = !:<■• 
''here c and j are undetermined functions of r and r. respectively; 

Neglecting powers of k above the second, we obtain 


8 / hi. 
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/0 0 0\ 
»j - fc I 0 0 £r J 
\0 0 ) 


'v 0 0 

-h fc 2 1 0 v 0 

fl 0 s + V 


where v denotes the derivative of v with respect to r and $ denotes the 
derivative of s with respect to z Hence 


and a 


h ~ k '{ V + r + ' +1 2 r, ) 
U - -±JbV 

A* 0 0\ 

coij = —\k* I 0 0 0 

\0 0 0 / 


'0 0 
0 0 
i0 n r 


3 




+ 2^t) 


X (v + * + s ^ 

+ (i X + 2 m -l n ) 


'(' + r + 0 

+G‘ + i-) 


+ 2/x 


+ 2 its 


x(. + ^ 4- S ^ 

+ (2 X+ " + 2 M ) 


J 

The matrix whose divergence must be equated to zero (to express the 

dt$> 

fact that the mass or body forces are zero) is — J*Ri* ^\here 

dv 

/0 -z 

Ri = Ez -j-fc( 2 0 0 J 4- ,»! = E$ 4* k 


3- 


/0 0 0 \ 

0 0 r - 

\0 0 0 / 
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Hence 

/0 -z 0\ 

HiJ - Ez + * [z 0 r ) + • 

\0 0 0 / 

/ 0 z 0\ 


J*R |* — 4- /' ( — 2 0 0 J + * • • and so 

V 0 r 0/ 


0 0 0 


10 0 r 
\0 r 0 


X (v' + £ + + 2 

/I 1 1 \ 

+ V2 X + 2 W ”47 


X -j- - 4* + 2/j - 

+ Qx + ). + i ”>) 


X fv' -}- “ + s'j + 2 fts' 

+ Qx + ^ + j».)r ! 


On equating the divergence of this matrix to zero, we obtain the 
following equations: 


(X + 2 ti) + ~ ~ ~ - X - m + | ii^j r, s" ~ 0. 

Since v is bounded at r = 0, we have 

„ (2a - X - 7« + $n)r 3 

•' " Cr + s STSo — 

■‘Vhoro C is an undetermined constant, s is a linear function of z, 
and we may set (why?) 

s »= Dz 


■vhere D is an undetermined constant. Since 7’* - ./ — is given bv 

dr) 


the formula 
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/0 0 0 \ 

, - MO 0 + 

\0 0 / 

X (V + ~ + ,) 4- 2 iaf 

+ Q x + r ‘- i n ) r 


X ^i>' + - + s'^ -f* 2/i 

+ Qx + ^ + i«) r 


X ^t' + - + 6 ^j + 2/ii 

+ G x + '‘ + i m ) r 


the fact that there is no applied force on the sides of the cylinder yields 
jf} 

2(X + (i )C + XD + 10(x + — |20X P + 2V + 4„m + (2X + *)»)- 0 

where R is the initial radius of the cylinder The axial coordinate of 
the force per unit initial area applied to the end z - l of the cylinder is 

the product of X ! by 2XC + (X + 2 „)D + 2J *“ + lc, ' j! ± SjOT+JiXn ^ 

S(X + 2 ft) 

and we choose the constants C and D so that the integral of this over 
the end of the cylinder is zero Thus w c obtain the relation 

2XC + (X + 2„)D + ‘■ HXl. + Ip? 2 + 8„m + 3Xn R , _ Q 
1G(X 4- 2 fi) 

When C and D satisfy the two linear equations just determined, the 
force applied to the end z = l is a couple about the axis of the cylinder 

of amount (prove this) Under this moment (or torque ) not 

only docs the cylinder rotate (as predicted by the linear theory) but 
also its radius increases by the amount 

A*/? (c + — ~ X — ± . fo? /?*] 

I 8(X 4- 2/i ) I 
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and its length increases by the amount JrDl. Since C and D are 
multiples of R' it follows that the increase of the radius of the cylinder 
is a multiple (possibly negative) of JrR 2 whereas the increase of the 
length of the cylinder is a multiple (possibly negative) of IrR*. Accur- 
ate measurements of the changes of the radii and lengths of cylinders 
(of reasonably large radius) would serve to determine the third-order 
elastic constants m and n (it being supposed that the second-order 
clastic constants X and n arc already known). 
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